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Preface 


This  research  was  undertaken  in  an  attempt  to  bring  the  worlds  of  active  controls  and 
structures  a  little  closer  together.  I  hoped  that  by  viewing  damping  from  an  active  controls 
point  of  view  and  from  a  structural  point  of  view  at  the  same  time,  I  could  make  a  little 
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by  my  results,  and  hope  that  you  find  this  dissertation  enlightening.  I  especially  hope  that 
my  approach  to  viscoelastic  damping  will  be  of  use  to  researchers  currently  struggling  with 
structures  containing  viscoelastic  materials. 
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Division,  Flight  Dynamics  Directorate,  Wright  Laboratories,  Major  A1  Janiszewski  (now 
Lieutenant  Colonel),  deserves  the  credit  for  planting  the  seed  which  grew  into  my  research 
topic.  To  his  successor.  Major  Steve  Whitehouse,  I  am  grateful  not  only  for  his  sponsorship, 
but  also  for  his  patience  and  understanding  when  things  took  a  little  longer  than  originally 
planned. 

My  committee  members  have  also  been  very  patient,  and  I  thank  them  for  sticking 
it  out.  I  am  especially  grateful  to  Dr  Liebst  and  Major  Warhola  for  the  many  hours  I 
spent  in  their  offices.  There  is  not  enough  room  to  express  my  appreciation  and  gratitude 
to  my  advisor.  Colonel  Bagley.  From  the  moment  he  first  suggested  I  consider  pursuing  a 
doctorate,  he  has  been  my  staunchest  supporter.  He  had  faith  in  me  and  my  abilities  at 
times  when  I  was  about  to  give  up  hope.  My  husband  Pierre  has  also  been  a  wonderful 
supporter.  He  has  done  more  than  his  share  in  taking  care  of  the  house  and  girls. 

This  dissertation  is  dedicated  to  my  daughters,  Marie  and  Liana,  both  of  whom  were 
born  during  its  development.  They  have  given  my  life  a  whole  new  meaning  and  purpose. 
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Abstract 

A  modified  form  of  the  standard  linear  quadratic  regulator  (LQR)  cost  functional  is 
used  to  optimally  blend  active  vibration  control  and  passive  structural  damping,  whether 
viscous  or  viscoelastic.  Viscoelastic  damping  is  first  modelled  by  a  standard  linear  model, 
and  then  by  a  fractional  derivative  model.  For  the  viscous  damping  case  and  the  classically 
modelled  viscoelastic  damping  case,  a  sub-optimal  closed  form  solution  is  derived  that  is 
independent  of  the  initial  conditions.  An  iterative  technique  that  minimizes  the  average 
value  of  the  cost  functional  and  one  that  minimizes  the  maximum  value  are  developed. 
Both  techniques  are  applicable  to  viscous  and  viscoelastic  damping  and  are  independent 
of  the  initial  state.  The  advantages  and  disadvantages  of  the  different  solution  techniques 
are  given 
examples 


with  respect  to  computation  requirements  and  performance.  Several  numerical 
illustrate  the  similarities  of  and  differences  between  the  various  techniques. 
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SIMULTANEOUS  DESIGN  OF  ACTIVE 
VIBRATION  CONTROL  AND 
PASSIVE  DAMPING 

/.  Introduction 

The  need  for  simultaneous  optimization  of  damping  and  active  vibration  control  is 
driven  primarily  by  the  potential  use  of  large  flexible  structures  in  space.  Simultane¬ 
ous  design  of  active  vibration  control  and  passive  damping  increases  the  effectiveness  and 
efficiency  of  vibration  suppression  systems  for  these  structures.  In  many  methods  of  vi¬ 
bration  suppression,  the  structure  is  made  as  rigid  as  possible,  the  residual  damping  is 
approximated  as  viscous,  and  then  an  optimal  controller  is  designed  with  respect  to  a  pre¬ 
determined  cost  function.  Although  the  controller  effectively  damps  vibrations  at  specific 
frequencies,  the  controller  may  excite  one  mode  while  trying  to  damp  out  another.  To 
stabilize  the  associated  modes,  passive  dampers  are  added.  But  this  changes  the  whole 
system  -  how  “optimum”  is  the  controller  now?  Is  control  effort  still  optimized?  What 
is  the  “best”  size  for  the  passive  dampers?  The  efficiency  as  well  as  the  effectiveness  of 
the  resulting  vibration  system  is  suspect.  By  designing  passive  and  active  control  ele¬ 
ments  simultaneously,  the  efficiency  and  effectiveness  of  the  vibration  suppression  system 
is  increased. 

Thus  a  need  exists  for  a  procedure  to  simultaneously  design  active  control  and  struc¬ 
tural  damping.  The  procedure  should  be  applicable  to  structures  with  viscous  or  viscoelas¬ 
tic  damping.  Mar  alluded  to  this  need  at  the  Vibration  Damping  1984  Workshop  when 
he  encouraged  viewing  damping  as  a  creative  force  in  design  [22].  Others  have  placed 
a  greater  emphasis  on  the  interaction  between  viscoelastic  damping  and  active  vibration 
control. 

One  of  the  efforts  in  this  area  was  the  development  of  numerical  schemes  for  the 
modeling  and  control  of  longitudinal  vibrations  in  a  rod  with  Boltzmann-type  viscoelastic 
damping  by  Burns  and  Fabiano  [7].  Another  effort  was  Hannsgen  and  Wheeler’s  research 
on  the  dynamic  behavior  of  viscoelastic  structures  which  emphasized  the  interaction  be¬ 
tween  passive  viscoelastic  damping  and  active  feedback  damping.  They  determined  the 
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interacting  effects  of  viscoelastic  dissipation  and  feedback  dissipation  in  the  damping  of 
oscillations  in  certain  viscoelastic  rods  and  beams  [15]. 

Simonian,  Major,  and  Gluck  pointed  out  that  incorporating  passive  damping  in  an 
active  control  design  must  be  done  so  that  the  deficiencies  of  one  technology  are  compen¬ 
sated  by  the  strengths  of  the  other  [29].  They  proposed  an  iterative  scheme  that  uses 
modal  strain  energy  analysis  to  determine  the  modal  damping.  A  hybrid  cost  function  was 
used  to  determine  the  tradeoffs.  Fowler  et  al.  implemented  a  variation  of  this  method  in  a 
computer  program  [14]. 

Smith,  et  al.  looked  at  redesigning  a  structure  to  make  it  easier  to  control.  They 
assumed  an  ideal  controller  had  been  designed,  and  adjusted  the  mass,  damping  stiffness, 
and  control  to  give  the  same  performance  with  less  effort  [30]. 

Onoda  and  Watanbe  used  a  direct  numerical  optimization  approach,  which  accounted 
for  the  uncontrolled  residual  modes,  for  the  design  of  an  optimal  controller  composed  of 
a  regulator  and  an  observer.  The  approach  was  incorporated  into  a  structure/controller 
simultaneous  optimization  scheme  [28].  Onoda  and  Watanbe  demonstrated  their  technique 
for  a  beam,  and  optimized  the  structure  by  changing  its  cross-sectional  area. 

This  dissertation  takes  a  different  approach,  attempting  a  closed  form  solution  to 
the  problem  of  how  to  achieve  an  optimal  blending  of  active  control  feedback  gains  and 
passive  structural  damping  parameters  (e.g.,  dashpot  constants,  thickness  of  viscoelastic 
pad,  or  area  of  contact  of  viscoelastic  pad).  A  variety  of  techniques  are  derived,  some  of 
which  are  closed  form  approximate  methods,  and  some  that  require  iteration  to  achieve 
the  optimum. 

'^he  techniques  are  developed  for  structures  in  which  the  equations  of  motion  can  be 
expressed  in  terms  of  finite  elements.  In  addition,  the  damping  matrix  in  the  equations  of 
motion  must  be  symmetric  and  be  able  to  be  expressed  as  a  product  of  three  matrices,  the 
middle  one  being  a  diagonal  matrix  of  the  damping  parameters. 

Chapter  II  contains  the  fundamental  development  of  the  dissertation.  Viscous  damp¬ 
ing  is  assumed;  Chapters  III  and  IV  will  extend  the  development  to  structures  with  vis¬ 
coelastic  damping.  Two  performance  indices,  each  a  modification  of  the  standard  linear 
quadratic  regulator  (LQR)  cost  functional,  are  used.  In  the  first  one,  the  damping  force 
enters  the  cost  functional  in  the  same  form  as  the  control  force.  This  cost  functional  leads 
to  a  closed  form,  approximate,  solution  for  the  damping  parameters  and  active  control 
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feedback  gains.  Next,  an  iterative  technique  for  determining  the  optimal  blend  of  damp¬ 
ing  parameters  and  active  control  feedback  gsuns  that  minimize  the  the  average  value  of 
the  cost  functional  over  a  representative  sample  of  initial  conditions  is  derived.  An  iter¬ 
ative  technique  that  minimizes  the  maximum  value  of  the  cost  functional  is  also  derived. 
Examples  are  included  to  illustrate  and  compare  the  techniques. 

In  the  second  performance  index,  damping  is  treated  as  a  one-time  cost.  Thus  the 
damping  parameters  appear  explicitly  in  the  cost  functional.  Although  a  closed  form 
solution  is  not  obtadned,  two  iterative  techniques  similar  to  those  derived  for  the  first 
modification  of  the  cost  functional  are  presented.  An  example  is  given  to  illustrate  these 
techniques. 

In  Chapter  III,  the  development  of  Chapter  II  for  the  first  performance  index  is 
extended  to  the  case  of  viscoelastic  damping.  Viscoelasticity  is  modelled  using  a  classical 
approach  in  which  the  relationship  between  stress  and  stain  is  expressed  in  terms  of  integer 
order  derivatives  on  stress  and  strain.  As  in  Chapter  II,  a  closed  form,  approximate, 
solution  for  the  damping  parameters  and  active  control  feedback  gains  is  derived.  The 
derivation  holds  for  any  (finite)  number  of  derivatives  on  stress  and  stridn.  For  the  case 
in  which  only  first  derivatives  on  stress  and  strain  are  used  in  the  model,  technlqnes  that 
minimize  the  maximum  value  of  the  cost  functional  and  the  average  value  of  the  cost 
functional  over  a  representative  sample  of  the  initial  conditions  are  derived.  Examples 
similar  to  the  ones  in  Chapter  II  are  presented. 

In  Chapter  IV,  the  development  of  Chapter  11  for  the  first  performance  index  is 
also  extended  to  the  case  of  viscoelastic  damping,  but  in  this  chapter  viscoelasticity  is 
modelled  using  fractional  order  derivatives  on  stress  and  strsdn.  Although  a  closed  form 
solution  is  not  obtained  in  this  ca'  ^.,  two  iterative  solutions  analogous  to  those  in  Chapter  II 
are  derived.  One  minimizes  the  maximum  value  of  the  cost  functional,  and  the  other 
minimizes  the  average  value  of  the  cost  finr^ional  over  a  representative  sample  of  the 
initital  conditions.  Examples  similar  to  the  ones  in  Chapters  II  and  III  are  presented. 

A  variety  of  solution  techniques  for  determining  the  optimum  blending  of  active 
vibration  control  and  passive  structural  damping  is  offered  to  the  design  engineer  in  this 
dissertation.  The  closed  form  solutions  give  him  quick,  adequate  results,  and  provide 
good  starting  points  for  the  iterative  techniques.  Minimizir'^  the  average  value  of  the 
performance  index  tends  to  give  better  pei’i.  •,ian''e  ol  the  two  iterative  techniques. 
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II.  Fundamental  Development 

2.1  Preliminaries 

In  addition  to  developing  techniques  for  simultaneously  designing  active  vibration 
control  and  passive  viscous  damping,  this  chapter  lays  the  groundwork  for  the  following  two 
chapters.  These  three  chapters  rely  heavily  on  a  modified  version  of  the  lineso*  quadratic 
regulator  (LQR)  performance  index,  which  is  commonly  used  in  control  theory  [1].  In 
this  section,  a  short  review  of  the  linear  quadratic  regulator  is  given.  The  next  section 
applies  the  development  of  the  linear  quadratic  regulator  to  the  problem  of  simultaneous 
optimization  of  passive  viscous  damping  and  active  control. 

In  this  thesis,  all  problems  considered  are  in  finite  element  formulation.  The  equations 
of  motion  of  a  structure  modelled  using  finite  elements  is 


Mx(t)  +  Dx(t)  +  Kx(t)  +  bnu(t)  =  0 


(1) 


where  x  is  the  position  vector,  M  is  the  mass  matrix,  D  is  the  damping  matrix,  K  is  the 
stiffness  matrix,  and  is  the  matrix  coefficient  of  the  vector  of  control  forces,  u. 

It  will  be  convenient  to  have  the  equations  of  motion  in  first  order  form: 


x(t) 

x(0 

v 


0  I 


-M-‘K  -M-^D 


+ 


0 

-M-*b„ 


u(0 


(2) 


By  dropping  the  notation  for  dependence  on  time  t  for  brevity,  and  defining  a  state  vector 
y  composed  of  the  position  and  velocity  vectors,  x  and  x 

y  =  ^ 

Elq  (2)  can  be  written  more  compactly 


x(t) 

x(0 


(3) 


y  =  Ay  +  Bu 


(4) 


where  A  is  the  plant  matrix  and  B  is  the  input  matrix.  The  initial  conditions  at  time  to 
are  assumed  known: 


y(<o)  =  yo 


(5) 
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Since  the  goal  of  this  thesis  is  vibration  suppression,  the  final  condition  at  a  specified  time 
t/  is  zero: 

yitf)  =  0  (6) 

The  LQR  performance  index  depends  on  the  state  y  and  on  the  control  u: 

J  =  J*'l(y'^qy  +  n'^Rn)dt  (7) 

The  state  weighting  matrix  Q  is  symmetric  and  positive  semi-definite,  while  the  control 
weighting  matrix  R  is  symmetric  and  positive  definite.  The  relative  magnitude  of  the 
diagonal  elements  of  the  weighting  matrices  Q  and  R  reflect  the  relative  importance  of  the 
state  y  and  the  control  u.  The  optimal  control  is  that  which  minimizes  the  performance 
index,  Eq  (7),  subject  to  the  state  equations,  Eq  (4). 

Using  the  method  of  Lagrange  multipliers,  the  state  equations,  Eq  (4),  can  be  ap¬ 
pended  to  the  performance  index  [12:263]: 

J  =  [^(y^’Qy  +  «"’»•«)  -  ^"■(y  -  Ay  -  bu)]  dt 

To  minimize  the  performance  index,  it  is  necec  sary  to  compute  the  variation  of  the 
performamce  index  due  to  variations  in  the  state  and  the  control. 

6J  =  [(y^Q^y  +  u’^Rfiu)  -  6\^{y  -  Ay-Bu)  -  \^{by  -  A6y  -  B^u)j  dt 

The  independence  of  the  three  unknowns  y,  u,  and  A  leads  to  three  equations: 


y  =  Ay  +  Bu 

(8) 

u  =  -R~*B’’A 

(9) 

0  =  A-l-A’^A-HQy 

(10) 

Assuming  the  Lagrange  multiplier  A  is  linearly  dependent  upon  the  state  y 

A=  Py 
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Eqs  (9)  and  (10)  can  be  written  as 


u=-R-‘B^Py  (11) 

0  =  Py  +  Py  + A^Py  +  Qy  (12) 

Substituting  Eq  (11)  into  Eq  (8),  and  then  Eq  (8)  into  Eq  (12),  one  obtains 

[P  +  PA  -  PBr  ‘b^P  +  A’^P  +  Qly  =  0 

This  equation  holds  true  for  general  y  only  if  the  expression  within  the  brackets  is  zero, 
which  in  turn  implies  that  P  satisfies  the  matrix  equation 

P  = -PA  +  PBR'^B^P  -  A^P  -  Q  (13) 

One  of  the  basic  characteristics  of  the  regulator  is  that  the  final  time  is  very  far 
in  the  future.  So  the  steady  state  solution  of  (13),  i.e.,  P  =  0,  is  of  interest  —  which 
implies  P  satisfies  the  algebraic  Biccati  equation 

PA  -  PBR-‘B’’P  +  A’’?  +  Q  =  0  (14) 

Since  the  algebraic  Riccati  equation  is  symmetric  (i.e.,  taking  the  transpose  of  Eki  (14) 
results  in  the  same  equation),  its  solution  P  is  symmetric  (P  =  P^).  Another  characteristic 
of  the  steady  state  solution  P  is  its  independence  of  initial  conditions.  The  control  u  that 
minimizes  the  performance  index  J  is  given  by  Eq  (11)  where  P  satisfies  Eq  (14). 

In  the  next  section,  viscous  damping  is  considered  analogous  to  active  control,  and 
the  problem  of  optimizing  damping  and  control  simultaneously  is  addressed. 

2.2  Simultaneous  Optimization  of  Viscous  Damping  and  Active  Control 

To  optimize  damping  and  control  simultaneously,  consider  a  modification  of  the  LQR 
cost  functional  to  include  passive  damping  forces.  The  passive  damping  forces  are  treated 
as  co-equals  to  the  active  control  forces,  and  are  weighted  in  a  similar  manner.  Let  v 
be  a  vector  of  passive  damping  forces  in  the  system,  and  S  its  positive  dehnite  weighting 
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matrix.  The  modified  performance  index  is  taken  to  be 

rtf  I 

J=  -(y^Qy  +  u^Ru  +  v’’Sv)d<  (15) 

Jto  2 

There  are  specific  reasons  for  weighting  the  passive  damping  forces  like  the  active 
control  forces.  When  implementing  a  damping  design  using  viscous  fluid  dampers  or  vis¬ 
coelastic  solids,  one  must  take  into  account  the  temperature  sensitivity  of  the  damping 
medium.  Modest  changes  in  temperature  due  to  absorbed  mechanical  energy  can  dra¬ 
matically  alter  the  damping  properties  of  both  fluids  and  solids.  Hence  one  is  motivated 
to  limit  in  some  fashion  the  mechanical  energy  absorbed  by  any  given  damper  as  well  as 
limit  the  peak  value  of  its  damping  force.  The  quadratic  damping  term  appearing  in  the 
performance  index,  Elq  (15),  serves  this  end  rather  nicely.  One  should  note  that  in  the  case 
where  the  passive  forces  v  are  proportional  to  the  velocity  vector,  this  quadratic  damping 
term  is  reminiscent  of  the  Rayleigh  dissipation  function  used  in  conjunction  with  energy 
methods  in  classical  dynamics  [25:89]. 

To  determine  the  form  of  v  and  how  it  fits  into  the  equation  of  motion,  a  closer 
look  at  the  equations  of  motion  is  needed.  In  the  problems  addressed  in  this  chapter,  the 
damping  matrix  is  symmetric  and  is  of  the  form 

D  =  byCbv’’  (16) 


where  C  is  a  diagonal  matrix  of  the  desired  viscous  damping  coefficients.  By  expressing  the 
damping  matrix  D  in  this  form,  the  damping  coefficients  have  been  separated  out  of  the 
damping  matrix.  Since  viscous  damping  coefficients  are  (usually)  constant  and  positive, 
C  is  positive  semi-definite.  The  matrix  b^  is  a  matrix  of  constants  chosen  to  reflect  the 
attachment  pattern  of  the  viscous  dampers.  Then  the  state  equations  can  be  written  as 


x(t) 

x(t) 

\ 


+ 


0  I 
-M-‘K  0 

0 

-M-‘bv 


x(t) 


x(/) 
Cb/x(t) 


0 

-M-^B 


u(0 


The  quantity  Cby^x  represents  the  passive  damping  forces. 


(17) 
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If  stiffness  is  allowed  to  vasy  as  a  design  parameter,  and  the  change  in  stiffness  can 
be  expressed  as 

AK  =  bfcCkb/ 

t 

then  the  third  term  on  the  right  hand  side  of  Eiq  (17)  becomes 


0  0 

-M-‘bk  -M-‘bv 


Ck  0 
0  C 


Hence,  the  passive  structural  forces  v  depend  linearly  upon  the  state  y  according  to 


V  =  C^y 


(18) 


where  C  contains  the  desired  passive  damping  coefficients  (and/or  changes  in  stiffness) 
along  its  diagonal.  The  matrix  ^  is  defined  as  either 


or 


bk^  0 
0  bv^ 


(19) 


depending  on  whether  stiffness  is  allowed  to  vary  or  not.  If  the  number  of  stiffnesses  that 
are  allowed  to  vary  plus  the  number  of  damping  coefficients  is  equal  to  the  length  of  the 
state  vector,  then  $  can  usually  be  chosen  such  that  it  is  invertible. 

By  dcrr.ing  the  matrix  By  in  a  similar  manner. 


By  = 


0 

-M-‘by 


or 


By  = 


0  0 

-M-‘bk  -M-^by 


the  state  equations  can  be  written  in  a  modified,  compact  form: 


(20) 


y  =  Ay  +  Bu  +  By  v 


(21) 


The  optimization  problem  can  now  be  expressed  as  the  minimization  of  the  perfor¬ 
mance  index,  Eq  (15),  subject  to  the  state  equations,  Eq  (21).  Two  additional  constraints 
are  that  the  passive  forces  v  are  of  the  form  given  in  Eq  ( 18)  and  the  matrix  of  damping 
coefficients  C  is  a  constant,  positive  semi-definite,  diagonal  matrix  since  damping  coeffi- 
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dents  are  non-negative.  The  method  of  Lagrange  multipliers  can  be  used  to  append  the 
first  two  constraints  to  the  performance  index: 

J  =  |^^(y^Qy  +  U^R-U  +  v^Sv)  -  A^^(y  -  Ay  -  Bu  -  BvV)  -  Aj  (v  -  C^y)  dt 

Now  the  vectors  y,u,  and  v  are  taken  to  be  independent.  So  Eq  (18)  implies  that  the 
matrix  C  is  independent  of  these  vectors  also.  By  definition,  Ai  and  A3  are  independent  of 
each  other  and  of  y ,  u,  v,  and  C.  Hence,  to  minimize  the  performance  index  J,  one  takes 
its  variation  with  respect  to  y ,  u,  v,  C,  At  and  A2  and  sets  the  variation  equal  to  zero.  The 
variation  of  the  performance  index  J  is 

i  J  =  [(y^Q^y  +  -1-  v^Srfv)  —  6X[ (y  -  Ay-Bu  -  Bvv) 

-Af(rfy  -  ASy  -  B^u  -  B^Sv)  -SA^(v  -  C^y)  -  A^(^v  -  SC^y  -  C^^y)]  dt 

Since  initial  and  final  conditions  sure  spedfied  (Elqs  (5)  and  (6)),  the  term  AfSy  can  be 
integrated  by  parts: 

/  XiSydt  =  A[tfy||^  -  /  XtSydi  =  -  /  XiSy dt 

Jt0  Jto  Jto 

Substituting  this  result  into  the  variation  of  the  performance  index  and  combining  Uke 
terms  yields 

6J=J*'  [(y^Q  +  a[  -1-  A[A  -I-  A^C$)^y  -t-  (u^R  -I-  Af  B)«u 

-l-(v^S  -H  AfBv  -  A2  )^v  -  6Ai(y  -  Ay  -  Bu  -  Byv) 

-«A^(v  -  C^y)  -t-  A[(«C)^y)  dt 

The  independence  of  the  six  unknowns  leads  to  six  simultaneous  equations 


y  =  Ay  -f  Bu  -f  BvV 

(22) 

V  =  C^y 

(23) 

u  =  -R-‘B’’A. 

(24) 

0  =  Ai  +  A’^Ax  +  Qy  +  ^^CAj 

(25) 

A2  =  By^Al  -f-  Sv 

(26) 

A2  =  0 

(27) 
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From  Eqs  (26)  and  (27)  the  passive  force  vector  v  is  proportional  to  the  Lagrange  multiplier 

Ai. 

V  =  (28) 

Substituting  Eqs  (24),  (27),  and  (28)  into  Eqs  (22)  and  (25)  yields  expressions  for  the  time 
derivatives  of  the  state  vector  y  and  the  Lagrange  multiplier  Ai: 

y  =  Ay  -  BR"‘B^Ai-BvS-‘Bv’^Ai 

Ai  =  -Qy  -  A^Ai 

From  Eqs  (23),  (27),  and  (28)  the  relation  between  the  parameter  matrix  C,  the  state 
vector  y,  and  the  Lagrange  multiplier  Ai  is 

C^y  =  -S'‘Bv^Ai  (29) 

Since  the  matrices  C,  S,  and  By  are  constant,  Eq  (29)  implies  that  the  Lagrange 
multiplier  Ai  is  equal  to  a  constant  matrix  times  the  state: 

A,  ~  Py  (30) 

where  the  square  matrix  P  is  constant.  This  is  the  form  of  the  standard  LQR  solution. 
In  fact,  Eqs  (24)  and  (30)  can  be  combined  to  ^ve  an  expression  for  the  active  control  u, 
which  is  the  same  as  that  derived  in  the  last  section: 

u  =  -R-‘B’’Py  (11) 

From  Eqs  (29)  and  (30)  the  relation  between  the  damping  coefficients,  the  state,  and  the 
constant  matrix  P  is 

C^y  =  -S"‘Bv^Py  (31) 

while  Eqs  (25)  and  (30)  yield  the  following  relationship  between  the  state  y,  its  time 
derivative  y,  and  the  matrix  P: 


0  =  Py  +  A^Py  +  Qy 


(32) 
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Substituting  for  the  time  derivative  of  the  state  y  in  Eq  (32)  yields  the  expression 

P(  Ay  +  Bu  +  BvV)  +  A^Py  +  Qy  =  0 

Using  Eqs  (11)  and  (18)  to  substitute  for  the  active  control  vector  u  and  the  passive  force 
vector  V,  one  obtains  an  equation  in  terms  of  the  matrix  P  and  the  state  vector  y: 

PAy  +  P[B(-R'‘B^Py)  +  Bv(-S"‘Bv^Py)]+A^Py  +  Qy  =  0 

This  equation  can  be  written  in  terms  of  a  matrix  times  the  state  vector: 

[PA  +  A^P  -  P(BR“‘B^  +  BvS-*Bv’’)P  +  Q]y  =  0 

For  general  y,  the  above  equation  leads  to  the  algebraic  Riccati  equation 

PA  +  A^P  -  P(BR-^B^  +  BvS-‘Bv’’)P  +  Q  =  0  (33) 

After  determining  the  solution  P  to  the  above  equation,  the  active  control  u  is  determined 
from  Eq  (11).  To  determine  the  viscous  damping  coefficients,  one  notes  from  Eq  (31)  that 
for  general  y, 

=  -S'^bJP  (34) 

If  ^  is  invertible,  it  would  seem  that  C  could  be  determined  from  Eq  (34) 

C  =  -S-‘Bv’’P^"‘  (35) 

where  P  satisfies  the  above  algebraic  Riccati  equation,  Eq  (33).  However,  in  general,  C, 
as  computed  from  Eq  (35),  is  fully  populated,  while  C  must  be  a  positive  semi-definite 
diagonal  matrix  as  the  problem  was  originally  posed.  This  requirement  was  not  specif¬ 
ically  addressed  in  the  above  derivation,  and  is  normally  not  encountered  in  an  LQR 
problem.  The  assumption  that  the  solution  has  the  form  of  the  standard  regulator  (i.e. 
Ai  =  Py)  is  what  led  to  a  non-diagonal  C.  So  the  optimal  solution  derived  in  this  section 
(Eqs  (11),  (33)),  and  (34))  is  not  attainable  since  it  led  to  a  fully  populated  C.  In  the  next 
section,  approximate  solutions  to  the  optimization  problem  of  Eq  (15)  are  proposed. 

Unlike  the  standard  LQR  problem,  the  solution  to  the  optimization  problem  posed 
by  Eq  (15)  and  C  constrained  positive  semi-definite  diagonal  is  not  given  by  a  Riccati 
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equation,  even  when  tf  —*  oo.  The  solution  to  the  present  problem  does,  however,  require 
an  iteration  involving  a  Riccati  equation  similar  to  Eq  (33)  and  is  dependent  upon  initial 
conditions.  It  is  presented  briefly  in  the  following  section. 

2.2.1  Solution  Technique  for  Minimizing  the  Cost  Functional  Error.  The  purpose 
of  this  section  is  to  find  good  approximate  solutions  to  the  problem  of  optimizing  active 
control  and  passive  viscous  damping  simultaneously.  Two  approaches  that  are  independent 
of  initial  conditions  are  given. 

One  approach  is  to  minimize  the  distance  between  the  optimal  solution  derived  in 
the  previous  section  and  the  set  of  allowable  solutions.  The  optimal  solution  is  the  matrix 
product  that  satisfies 

C4f  =  -S-‘Bv^P  (34) 

where  P  is  the  solution  of  Eq  (33).  The  set  of  allowable  solutions  is  the  set  of  matrices 
of  the  form  C#  where  C  is  a  positive  semi-definite  diagonal  matrix  and  ^  is  defined 
by  Eq  (19).  This  set  is  a  convex  subset  of  the  space  of  all  matrices  that  have  the  same 
dimensions  as  C^.  The  distance  will  be  defined  as  the  Frobenius  norm  of  the  difference: 

|1C#-|-S''B^P||f  (36) 

Minimizing  this  norm  is  equivalent  to  a  least  squares  fit.  This  norm  will  be  minimized 
by  computing  the  appropriate  viscous  damping  coefficients,  which  are  the  elements  of  the 
diagonal  matrix  C. 

The  Frobenius  norm  of  an  m  X  n  matrix  A  is  defined  as  the  square  root  of  the  sura 
of  the  squares  of  the  elements  of  the  matrix: 


I|a|1f  = 


n  1/2 


EE  4 

i=i  j=i 


Minimizing  Eq  (36)  over  the  set  of  positive  semi-definite  diagonal  matrices  is  a  minimization 
problem  that  can  be  solved  in  closed  form  (see  Appendix  A).  If 


W  =  -S-‘b/P 
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the  diagonal  elements  of  the  diagonal  matrix  C  are 

If  m  is  the  number  of  damper  coeilicients,  and  n  is  the  length  of  y,  ^  and  W  are  m  x  n 
matrices. 

The  second  approach  involves  examining  the  minimum  performance  index  (Eq  (15)) 
when  the  damping  coefficients  C  are  set  before  computing  the  active  control  gains  rather 
than  computing  them  simultaneously.  The  performance  index  becomes 

J  =  J*'  ^  [y^(Q  +  #^CSC#)y  +  u’’Ru]  di  (38) 

subject  to  the  equations  of  motion: 

y  =  (A  +  BvC^)y  +  Bu  (39) 

With  the  damping  coefficients  C  fixed  in  advance,  the  minimlzaton  is  a  standard  LQR 
problem  with  modified  state  weighting  matrix  Q  amd  modified  plant  matrix  A.  The  mini¬ 
mum  performance  index  J  is  achieved  by  Blq  (11): 

u  =  -R-^B^Py  (11) 

where  P  (in  the  limit  as  tf  -*  oo)  is  a  solution  to  the  algebraic  Riccati  equation 

P(A  4-  BvC^)  +  (A  +  BvC^)’’P  -  PBR-'B^P-I-Q  +  ^’’CSC^  =  0  (40) 

This  algebrmc  Riccati  equation  expresses  the  matrix  P  :n  terms  of  the  matrix  C. 
In  effect,  the  active  control  u  is  now  in  terms  of  the  damping  coefficients  C.  The  next 
step  is  to  determine  how  to  select  the  damping  coefficients  C  such  that  the  cost  functional 
(Eq  (38))  is  minimized. 

One  approach  is  to  express  the  value  of  the  performance  index  as  ->  oo  in  terms 
of  P  and  the  initial  conditions  (see  Appendix  B): 

J  =  J{tj)  =  ^yo  Pyo  (41) 
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Thus,  if  the  initial  conditions  yo  are  known,  one  can  iterate  with  respect  to  the  damping 
coefficients  (i.e.,  positive  semi-definite,  diagonal  C),  using  Eq  (40)  to  determine  P  during 
each  iteration,  to  minimize  the  performance  index  J  given  by  Eq  (41).  However,  opti¬ 
mizing  for  fixed  yo  may  result  in  poor  performance  if  initial  conditions  change  (due  to  an 
unexpected  disturbance,  for  instance).  In  contrast,  the  standard  LQR  problem  does  not 
require  this  iteration  because  the  solution  to  the  Riccati  equation  (Eq  (14)  is  independent 
of  the  initial  conditions  yo,  and  hence  results  in  inherently  robust  performance  with  respect 
to  the  initial  conditions.  To  achieve  robustness  with  respect  to  initial  conditions,  solutions 
independent  of  yo  will  be  sought. 

It  will  be  helpful  to  consider  the  difference  between  the  value  of  the  performance  index 
at  the  (unattainable)  solution  of  Section  2.2  and  that  achieved  by  enforcing  the  constraint 
that  the  damping  coefficient  matrix  C  must  be  diagonal.  Letting  P  represent  the  uncon¬ 
strained  optimal  solution  and  P  a  constrained  solution,  the  change  in  the  performance 
index  is 

AJ  =  yo  Pyo  -  yo  Pyo  =  y?(P  -  P)yo  =  (42) 

In  solving  for  AP,  one  first  recalls  that  P  satisfies  Eq  (33) 

PA  +  A^P  -  P(BR‘^B’’  +  BvS-‘Bv’’)P  -|-  Q  =  0  (33) 

while  P  satisfies  Eq  (40): 

P(A  +  BvC#)  +  (A  -I-  ByC^fP  -  PBR  'b’^P-I-Q  +  =  0  (43) 

Subtracting  Eq  (33)  from  Eq  (43)  leads  to  an  expression  for  AP 

APA  -1-  A^AP  -t-  APBvC#  -1-  PB^C$-l-(BvC^  -I-  PBvC#)^AP 
+  (BvC^  +  FByC^fP  -  AP(BR~‘B^P)  -  PBR“^B’’P 

-  AP(BR“^B’^)AP-l-^^CSC^-bPBvS-'Bv^P  =  0 

which  can  also  be  written  in  the  form  of  an  algebrmc  Riccati  equation 

AP(A  -  BR~‘B^P  -I-  B,C#)  -1-  (A  -  BR"‘B^P-|-BvC$)^AP 

-  AP(BR“^B^)AP-t-E  =  0  (44) 
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where 

E  =  #’’CSC^  +  ^^CBv’’P  +  PBvC^  +  PBvS-‘Bv^P  (45) 

This  Riccati  equation  determines  AP,  which  in  turn  determines  AJ.  Standard  Ric- 
cati  solvers  find  the  real,  symmetric  solution  to  the  algebraic  Riccati  equation  that  stabilizes 
the  plant  matrix.  Now  the  plant  matrix  (A  —  BR~^B’’P  +  ByC^)  is  stable.  Therefore, 
if  E  =  0,  then  AP  =  0  is  the  only  real  symmetric  solution  to  Eq  (44)  [36].  So,  driving  E 
as  small  as  possible  yields  a  small  AP  and  thus  drives  P  to  P.  This  is  highly  desirable 
because  the  unconstrsdned  standard  LQR  solution  P  is  highly  robust  with  respect  to  initial 
conditons. 

So  this  approach  to  an  approximate  solution  is  based  on  finding  those  damping 
coefficients  that  make  the  matrix  E  as  close  to  the  zero  matrix  as  possible.  Based  on 
Eq  (45),  E  can  be  written  as 

E  =  #’’C(SC^  +  By^P)  +  PByS"^(SC#  +  By^’P) 

=  (^^C  +  PByS-^)(SC*  +  By^P) 

=  (*^CS  +  PBy)S“^(SC^  +  By^'P) 

=  [S~^/’(SC*  +  By’'P)f  [S'‘'’(SC^  +  By^P)] 

=  H^H  (46) 


where 

H  =  S~^/*(SC#  +  By^P) 

So  choosing  the  damping  coefficients  such  that  the  diagonal  matrix  C  makes  the  matrix 
product  SC#  as  close  to  the  matrix  product  (— By’^P)  as  one  can  will  make  E  small. 
Since  ||H’’H||2  =  HHH*,  choosing  the  damping  coefficients  that  minimizes  the  two-norm  of 
H  minimizes  the  two-norm  of  E  and  drives  AJ  small.  This  minimization  can  be  accom¬ 
plished  using  the  MATLAB  routine  EMINS,  which  uses  the  Nelder-Mead  simplex  search 
algorithm  [23]. 

It  is  interesting  to  note  how  similar  the  matrix  H  is  to  the  error  term  in  Eq  (36). 
In  fact,  minimizing  the  Frobenius  norm  of  Eki  (36)  gives  the  same  result  as  minimizing 
mw.  Because  ||H^H||p  <  ||H|||.,  choosing  the  damping  coefficients  that  minimize  the 
Frobeniuo  norm  of  H  will  minimize  an  upper  bound  on  the  Frobenius  norm  of  E  and  thus 
drive  AJ  — ►  0.  The  advantages  of  the  Frobenius  minimization  over  minimizing  |1H)|2  are 
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that  it  is  guaranteed  to  return  positive  damping  coefficients  and  it  can  be  solved  in  closed 
form  (see  Appendix  A).  As  will  be  shown  later,  the  latter  advantage  madies  the  Frobenius 
minimization  computationally  more  efficient. 

In  this  section,  two  methods  of  obtaining  an  approximate  solution  to  the  optimization 
problem  posed  in  Section  2.2  have  been  derived.  One  is  a  minimum  singular  value  method 
(minimize  ||H||3),  and  the  other  is  a  Frobenius  norm  minimization  (minimize  ||H||f ).  How¬ 
ever,  these  are  approximate  solutions.  It  is  desirable  to  achieve  am  actual  optimization.  In 
the  next  two  subsections,  the  performance  index  is  redefined  in  two  different  ways,  and  the 
minimums  of  ihe  resulting  cost  functionals  sought. 


2.2.2  Solution  Technique  for  Minimizing  the  Mean  Performance  lex.  In  each 
of  the  next  two  subsections,  a  cost  functional  is  defined  based  on  the  form  of  the  perfor¬ 
mance  index  given  in  Eq  (41).  The  cost  functional  of  the  next  section  is  the  maximum 
value  of  the  performance  index.  In  this  subsection,  a  cost  functional  is  the  average  of  the 
performance  index  over  the  initial  conditions.  Since  the  magnitude  of  the  initial  condi¬ 
tions  is  irrelevant  to  the  value  of  viscous  damping  coefficients  C  and  the  corresponding 
solution  to  the  modified  Riccati  equation,  Eq  (40),  initial  conditions  used  in  the  next  two 
subsections  are  those  whose  magnitude  is  unity,  in  other  words,  the  unit  ball  ||yo||  =  1. 

Minimizing  the  average  value  of  the  performance  index,  denoted  by  J,  over  the  unit 
ball  ||yo||  =  1  does  have  the  potential  to  lead  to  large  values  of  the  performance  index  for 
some  initial  conditions,  but  the  large  values  will  be  offset  by  low  values  of  the  performance 
index  at  other  initial  conditions.  Since  the  damping  coefficients  and  active  control  gains 
that  are  calculated  will  give  the  lowest  average  value  of  the  performance  index,  one  would 
expect  that  this  solution  will  give  better  performance  over  more  than  half  of  the  initial 
condition  compared  to  any  other  combination  of  damping  coefficients  and  active  control 
gains. 

The  average  value  of  the  i>er.^ormance  index  will  be  calculated  by  integrating  the 
performance  index  over  the  surface  of  the  unit  ball,  then  dividing  by  the  surface  area  of 
the  unit  ball.  First  the  performance  index  will  be  expressed  in  terms  of  the  elements  of 
the  initial  state  vector  and  the  solution  to  the  Riccati  equation.  Let  yo  =  [j/i  ^2  •  •  *  VnY  • 
Since  P  is  symmetric,  the  performance  index  is 


J  =  ^yo  Pyo  =  ^  E  E  PayiVi  =  ^  E  P»y^i  +  S  S 


i=i 


(47) 


i  i=i+i 
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The  expression  for  the  average  value  of  the  performance  index  is 


j  4.1.,  (i  Sr.,  ‘‘am)  •I* 

Since  the  sums  are  finite,  the  integrals  and  summations  can  be  interchanged: 


J  = 


(/||yo||=l  ^  ^  (ii|yo||=l 

t  i=i+l 


i=l 


•^llyoM=i 


(48) 


To  simplify  the  evaluation  of  the  integrals  in  Eq  (48),  Gauss’  divergence  theorem  is  used 
to  express  the  surface  integals  as  volume  integrals  [19:355]: 


J  n'udA^  J  divudV 


(49) 


Let  e*  represent  the  unit  vector  in  the  direction  of  the  kth  coordinate.  Then  the  unit  vector 
normal  to  the  surface  of  the  unit  ball  is  n  =  Vk^k-  The  vector  u  must  be  determined 
for  each  of  the  three  integrals  in  Eq  (48).  (Please  note  that  in  this  discussion,  u  is  some 
continuous  vector  function,  not  the  active  control  vector.)  If  n  •  u  is  taken  to  be  yf,  one 
can  infer  that  u  =  52fc=i  =  Vi®**  Therefore  the  divergence  of  the  vector  u  is  unity: 


divu 


dVi 

dyi 


1 


(50) 


So  evaluation  of  the  first  integral  in  the  numerator  in  Ek}  (48)  is  the  volume  of  the  unit 
ball: 

/  y?dA=  /  dV^V  (51) 

•'llyo||=i  ''l|yo||=i 

In  the  second  integral  in  the  numerator  of  Eq  (48),  n  •  u  is  taken  to  be  y^y,- .  It  follows 
that  u  =  2^yi^j  +  yj^i)-  Then  the  divergence  of  u  is  zero: 


div  u 


^Vk  dyi 


+ 


dyj 


=  0 


(52) 
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Hence  the  second  integral  in  the  numerator  of  Exi  (48)  integrates  to  zero: 


I 


|y«ll=» 


Vit/j  dA  =  0 


(53) 


The  integral  in  the  denominator  of  Eki  (48)  is  just  the  surface  area  of  the  unit  ball,  but  it 
can  be  expressed  in  terms  of  the  volume.  In  the  integral,  the  quantity  n  •  u  is  unity: 

n  •  u  =  1 


It  is  sufficient  that  u  =  n.  Thus  the  divergence  of  u  is  equal  to  the  number  of  terms  in  the 
state  vector: 

JL.  fiu.  JL  Bt,. 

(54) 


ifc  9yk  dyt 


Thus  the  surface  area  of  the  unit  ball  is  n  times  the  volume: 


/  dA=  [  ndV  = 


nV 


Combining  Eqs  (48),  (51),  (53),  and  (55), 


(55) 


<=1  i  i=«+i 


=  — trace  P 
2n 


(56) 

(57) 


Hence,  choosing  the  viscous  parameters  C  that  minimize  the  trace  of  the  matrix  P  mini¬ 
mizes  the  average  value  of  the  performance  index  J  over  the  unit  bail  ||yo||  =  1,  regardless 
of  system  order. 

This  subsection  has  derived  a  technique  for  minimizing  the  average  value  of  the  cost 
functional  over  a  representative  sample  of  initial  conditions.  The  next  subsection  will 
present  another  technique  for  optimizing  viscous  damping  and  active  control  simultane¬ 
ously. 


2.2.S  Solution  for  Minimizing  the  Maximum  Performance  Index.  Another  ap¬ 
proach  to  optimizing  viscous  damping  and  active  control  simultaneously  is  to  minimize 
the  maximum  value  of  the  performance  index  J  over  the  unit  ball  ||yo||  =  1.  One  way  to 
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accomplish  this  is  to  minimize  the  two-norm  of  P,  since  the  performance  index  is  less  than 
or  equal  to  the  two-norm  of  P: 

•/  =  jlyfPy.l  <  jllyof  l|P|l.  =  jIIPII.  m 

The  quantity  on  the  right  hand  side  of  the  above  equation  will  be  referred  to  as 

=  i||pii. 

In  this  method,  iteration  with  respect  to  the  viscous  parameters  C  is  carried  out  until 
||P||3  is  minimized  where  P  satisfies  Eki  (40).  This  approach  might  be  very  conservative  in 
general,  since  the  initial  conditions  yo  that  maximize  the  performance  index  may  not  be 
encountered  very  often. 

The  last  three  subsections  have  produced  four  solution  techniques  that  are  indepen¬ 
dent  of  initial  conditions.  In  summary,  the  four  solution  methods  are 

1.  min  llH||f :  Closed  form  solution  based  on  finding  a  diagonal  C  that  minimizes  the 
error  in  the  performance  index  with  respect  to  the  FVobenius  norm.  See  Eiqs  (40)-(46). 

2.  min  ||H||3:  Iterative  solution  technique  based  on  finding  a  diagonal  C  that  minimizes 
the  error  in  the  performance  index  with  respect  to  the  two-norm.  See  EJqs  (40)-(46). 

3.  min  ||P||3:  Iterative  solution  technique  based  on  finding  a  diagonal  C  that  minimizes 
the  maximum  value  of  the  performance  index.  See  Ekjs  (40)  and  (58). 

4.  min  (trace  P):  Iterative  solution  technique  based  on  finding  a  diagonal  C  that  min¬ 
imizes  the  average  value  of  the  performance  index.  See  Eqs  (44)  and  (47)-(57). 

The  minimization  of  HHHf  and  the  minimization  of  ||H||3  do  not  guarantee  a  positive 
solution  for  the  damping  parameters.  The  same  applies  to  minimizing  trace(P)  and  ||P||j. 
However,  as  will  be  seen  later,  positive  damping  parameters  are  obtained  in  all  the  example 
problems.  Following  are  example  problems  which  illustrate  the  above  techniques.  The 
example  problems  are  also  used  to  compare  the  different  solution  techniques. 

2.2.4  Example  Problems. 

2.2.4. 1  Example  #1  -  Single  DOF  System.  In  this  example  problem, 
single  degree  of  freedom  spring-mass-damper  system,  illustrated  in  Figure  1,  is  used  to 
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u 

Figure  1.  Simple  Spring-Mass- Damper  System 

demonstrate  and  compare  the  four  solution  techniques  derived  in  this  section.  The  equation 
of  motion  is  given  by  the  scalar  equation 

X  cx-\-  (0.5  -I-  Ak)x  =  —u  (59) 

« 

where  AA;  and  c  ate  the  passive  design  parameters,  representing  change  in  stiffness  and 
viscous  damping  respectively.  Letting  yj  =  *,  ya  =  x,  then 


=  Ay  -I-  B«  +  ByC^y 


The  weighting  matrices  are  chosen  to  give  equal  weighting  on  position,  active  control,  and 
viscous  damping: 

Q  = 

The  weight  on  the  change  in  stiffness,  Afc,  is  chosen  large  to  assure  negligible  change  in 
stiffness  {Ak  «  0).  This  isolates  damping  as  the  design  parameter  of  interest. 
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Figure  2.  Performance  Index  vs.  Initial  Conditions 


The  optimal  unconstrained  solution  of  Section  2.2  (i.e.,  the  solution  of  £q  (35)  )  for 
this  problem  yields  a  non-diagonal  parameter  matrix: 


C  = 


0.000  0.000 

.500  .707 


(61) 


For  the  single  degree  of  freedom  case,  ||yo||  =  1  can  be  represented  by  a  single 
parameter  6  where  yo  =  [cost*  sin^]^.  The  value  of  the  performance  index  for  the  optimal 
unconstrained  parameter  matrix  is  shown  as  a  solid  line  with  circles  in  Figure  2.  Also 
included  for  comparison  is  the  result  for  purely  active  vibration  control  (i.e.,  C=0). 

Table  1  summarizes  various  results  for  the  single  degree  of  freedom  system.  For  this 
example,  the  results  for  the  alternative  suboptimal  solutions,  the  min  ||H||2  amd  min  ||H||f 
cases,  are  the  same  to  two  decimal  places,  but  in  general  the  results  will  be  different. 
Although  the  damping  coefficients  for  the  min  ||H||2  and  min  ||H||f  cases  are  nearly  twice 
those  for  for  the  min  ||P||2  and  min  (trace  P)  cases,  they  nonetheless  serve  as  reasonable 
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The  result  is  plotted  in  Figure  2.  The  graph  shows  this  solution  to  be  more  robust 
than  the  purely  active  solution,  r'or  at  least  one  initial  condition,  any  other  value  of 
damping  coefficient  (or  change  iu  stiffness)  will  produce  a  larger  value  of  the  performance 
index  than  the  maximum  uf  this  solution. 

Minimizing  the  trace  of  P  for  Example  #1  gave  the  following  diagonal  parameter 
matrix: 

0.000  0 

C  = 

0  .455 

These  results  are  also  shown  in  Figure  2.  Although  this  graph  does  have  points  above  the 
one  for  min  ||P||3,  the  average  value  of  the  performance  index  is  lower  for  this  solution 
than  for  any  other.  The  min  ||H||2  and  the  min  ||H||f  solutions  are  equal  to  two  decimal 
places,  so  only  the  min  ||H||f  solution  is  shown  on  Figure  2. 

In  this  example,  the  different  solution  techniques  were  compared  using  the  original 
performance  index  (Eq  (15)).  In  the  following  example,  three  solution  techniques,  min 
IIHII;.,  min  IIPII2,  and  min  (trace  P)  are  compared  using  the  modified  performance  indices 
introduced  in  Section  2.2.2. 

2.2.4. S  Example  #2  -  Single  DOF  System  xvith  No  Additional  Stiffness.  In 
the  example  just  presented,  the  change  in  stiffness  {Ak)  was  made  small  by  choosing  the 
corresponding  weight  large.  The  following  example  constrains  the  change  in  stiffness  to  be 
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zero  from  the  outset  (Ak  =  0).  The  scalar  equation  of  motion  is 


i  +  ci  +  0.5a:  =  — u 


which  lesuls  to  a  first  order  system  in  which  the  parameter  matrix  is  a  scsdax  and  the  matrix 
^  is  nonsquare: 


=  Ay  +  Btt  +  BvC^y  (64) 


Before  giving  the  weightings  on  the  state,  control,  and  passive  damping,  a  short 
discussion  on  how  to  choose  the  weightings  is  appropriate.  The  weightings  determine  the 
relative  importance  of  driving  the  state  to  zero,  the  expenditure  of  control  energy,  and 
the  absorption  of  damping  energy  [26:782].  In  the  the  previous  example,  no  weighting  was 
placed  on  velocity,  implying  that  driving  velocity  to  zero  was  not  important.  Choosing  the 
weightings  on  the  state,  control,  and  passive  damping  is  “more  of  an  art  than  a  science” 
[18:219].  In  this  example,  the  weighting  matrices  are  chosen  to  minimize  energy.  Hence, 
velocity  is  weighted  more  heavily  than  position  since  the  magnitude  of  the  mass  is  greater 
than  the  magnitude  of  the  stiffness.  Active  control  and  passive  damping  are  weighted 
equally,  and  less  than  position: 


.5  0 

Q=  ,  R  =  .l,  S  =  .l 

0  1 


(65) 


In  practice,  the  design  engineer  may  have  to  experiment  with  different  weightings 
to  get  a  design  that  is  appropriate  to  his  needs.  For  example,  if  the  weightings  on  active 
control  and  passive  damping  are  too  large,  the  damping  ratios  may  be  unacceptably  low.  If 
the  weightings  on  active  control  and  passive  damping  are  too  small,  the  amount  of  control 
and  damping  material  required  may  be  unacceptably  high.  For  a  more  detailed  discussion 
of  the  choice  of  weighting  matrices,  see  Reference  [2]. 
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Table  2.  Active  Control  Feedback  Gains  for  SDOF  System 


c 

feedback  matrices 

m 

(1.72  2.61] 
(1.72  3.00] 
(1.72  3.27] 

Table  3.  Comparison  of  Solutions  for  SDOF  System 


c 

2.52 

itlAl 

Esa 

trace  P 

1.31 

ng 

Ba 

wsm 

1.09 

m 

In  this  case,  the  optimal  unconstrained  solution  of  Section  2.2  (i.e.,  the  solution  of 
Eq  (34))  for  this  problem  yields  a  1  x  2  matrix: 

c*  =  [  1.35  2.52 

The  E^obenius  norm  solution  yields  c  =  2.52,  while  minimizing  trace  P  gives  c  =  0.89,  and 
minimizing  ||P||3  g^ves  c  =  0.47.  The  active  control  feedback  gains  are  given  in  Table  2. 

The  disparity  of  the  values  for  the  damping  parameter  c  indicate  that  there  is  some 
insensitivity  to  the  magnitude  of  the  damping  parameter  c.  This  is  supported  by  the 
fact  that  there  is  little  difference  in  trace  P  and  ||P||3  for  the  three  solution  techniques 
(Table  3).  Figure  3  shows  that  there  is  little  change  in  ||P||3  due  to  the  value  of  c  between 
c  =  0  and  c  =  1.  The  graph  of  the  trace  of  P  is  varies  by  less  than  0.05  between  c  =  0.5 
and  c  =  1.5. 

2.2. 4 ’3  Example  #3  -  Fourteen  DOF  System.  In  this  example,  a  more 
realistic  and  more  complex  structure  is  considered,  the  plane  aluminum  truss  in  Figure  4. 
The  finite  element  model  equations  of  motion  can  be  written  in  the  state  space  form: 


y  =  Ay  +  Bu  +  ByC^y  (66) 

Since  the  fourteen  physical  degrees  of  freedom  lead  to  twenty-eight  mathematical  degrees 
of  freedom  in  a  dynamical  system,  the  dimension  of  state  vector  y  is  28x1,  There  are  four 
control  actuators  and  five  dashpots  on  the  truss,  so  the  dimension  of  the  active  control 
vector  u  is  4x1  and  the  dimension  of  the  damping  parameter  matrix  C  is  5x5.  The 
dimensions  of  the  remaining  matrices  are  given  as  follows:  plant  matrix  A,  28x28;  active 
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c 


Figure  3.  Example  Problem  #2  —  Trace  P  and  ||P||2  vs.  c 

control  input  matrix  B,  28x4;  passive  force  input  matrix  By,  28x5;  and  5x28.  With  E 
representing  the  elastic  modulus  of  aluminum,  A  the  cross  sectional  area  of  the  aluminum 
elements,  and  p  the  density  of  the  aluminum,  the  parameters  of  the  system  are 

EA  =  1.5179  xlO«  N  (3.4125  xlO®  lb) 

pA  =  8.3564  N-s*/m^  (1.212  Xl0“®  lb-s*/in*) 

The  length  of  the  horizontal  elements  is  4.572  m  (15  ft)  while  the  vertical  elements 
are  3.048  m  (10  ft).  The  length  of  the  two  elements  with  dampers  and  Cs 
is  4.819  m  (15.81  ft).  Thus,  the  diagonal  elements  with  dampers  Ci,  C2,  and  C3 
are  5.495  m  (18.03  ft)  in  length.  In  addition  to  the  consistent  mass  matrix  due  to  con¬ 
sideration  of  the  mass  of  the  aluminum  elements,  there  are  0.454  kg  (1  Ibm)  point  m<isses 
at  every  node,  with  an  0.907  kg  (2  Ibm)  point  mass  at  the  tip.  The  undamped  natural 
frequencies  range  from  7.7  rad/sec  to  319.3  rad/sec  (1.2  hz  to  50.8  hz). 

The  state  weighting  matrix  Q  is  chosen  such  that  y^Qy  equals  the  total  mechanical 
energy  in  the  system.  Thus  Q  is  formed  using  the  mass  and  stiffness  matrices,  M  and  K. 
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«3 

Figure  4.  Example  Problem  #3  Truss  Structure 


The  active  control  weighting  matrix  is  chosen  such  that  its  diagonal  elements  are  smaller 
thaa  the  smallest  diagonal  value  of  Q. 

In  the  previous  two  examples,  the  weighting  on  active  control  and  passive  damping 
were  equal.  In  Example  #1,  it  was  shown  that  using  simultaneously  designed  active 
feedback  control  and  passive  damping  gave  better  performance  than  using  active  control 
alone.  In  general,  one  would  expect  that  passive  damping  is  cheaper  than  active  control 
-  it  doesn’t  need  an  actuator  or  power  for  the  actuator.  Hence  for  this  example,  it  is 
assumed  that  active  control  is  more  costly  than  passive  damping,  so  the  passive  damping 
weighting  matrix  S  is  chosen  to  be  10“^  times  identity,  and  the  active  control  weighting 
matrix  R  is  chosen  to  be  10“*  times  identity: 

K  0  1 

Q  =  ,  R  =  10-^1,  S  =  10-®I  (67) 

0  M 

Using  the  weightings  given  above  gives  the  resulting  damping  parameters  for  the 
different  solution  techniques  in  Table  4.  At  the  bottom  of  Table  4  are  the  approximate 
computation  times  on  a  VAX  6420.  The  solution  techniques  can  be  run  with  simple  MAT- 
LAB  routines.  The  Frobenius  norm  minimization  uses  the  algorithm  given  in  Appendix  A, 
while  the  other  minimizations  utilized  the  MATLAB  routine  FMINS,  using  default  toler¬ 
ances  in  the  minimization  routines.  The  minimization  algorithm  used  in  FMINS  is  the 
Nelder-Mead  simplex  search  algorithm  [23], 
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Table  4.  Damping  Parameters  and  CPU  Times  for  14  DOF  Truss 


min  H  3 

mino  ||H||j 

■BM 

60.0169 

59.6515 

59.6514 

dashpot 

WM 

41.7915 

46.2314 

46.2314 

constants 

■■ 

41.9507 

14.7370 

14.7372 

mm 

57.8418 

59.2816 

59.2816 

^9 

55.2145 

45.9636 

45.9636 

computation  time 

17  s 

34  s 

61  s 

min  trace  P 

mino  trace  P 

msmm 

Cl 

54.3826 

54.3826 

76.1842 

76.1839 

dashpot 

Cj 

35.4502 

35.4503 

17.5549 

17.5542 

constants 

C3 

38.4443 

38.4443 

28.5870 

28.5873 

C4 

42.1846 

42.1846 

34.2785 

34.2786 

C6 

46.8419 

46.8420 

45.5261 

45.5255 

computation  time 

22  min 

52  min 

38  min 

78  min 

The  zero  subscript  on  the  minimization  operator  (e.g.,  mino)  indicates  that  an  initial 
guess  of  C  =  0  was  used.  The  other  cases  used  the  damping  coefficients  of  the  min  ||H||f 
solution  as  the  initial  guess.  Using  the  min  ||H||f  solution  as  an  initial  guess  significantly 
decreased  the  computation  time,  by  more  than  one  half  in  some  cases.  The  active  gain 
feedback  matrices  are  given  in  Appendix  F.  The  results  in  Table  4  clearly  demonstrate 
that  if  one  chooses  not  to  spend  the  additional  computer  time  to  optimize  the  dampers 
that  the  inexpensive  min  solution  is  a  reasonably  close  estimate.  For  very  large 

systems  this  solution  may  be  the  only  affordable  alternative. 

In  fact,  as  shown  in  Table  5,  the  min  ||H||f  solution  gave  almost  identical  closed 
loop  damping  ratios  and  natural  frequencies  as  the  min  (trace  P)  and  min  ||P||3  solutions 
for  this  example.  Also  shown  in  Table  5  are  the  open  loop  poles  of  the  system  (no  active 
or  passive  control).  The  average  chsuige  in  frequency  due  to  adding  viscous  damping  and 
active  feedback  damping  was  about  3  rad /sec. 

Since  it  is  not  possible  to  show  a  two-dimensional  plot  of  J  over  the  unit  ball  as 
in  Example  1,  another  method  of  comparing  the  solutions  is  needed.  It  would  be  useful 
to  have  a  way  to  compare  the  solutions  with  respect  to  the  original  performance  index, 
Eq  (15),  regardless  of  system  order.  To  compare  the  solution  given  in  Table  4,  recall  from 
Eq  (42)  that  the  difference  in  J  between  two  different  solutions  can  be  given  by 


AJ  =  yo(Pi  -P2)yo 


(68) 
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Table  5.  Damping  Ratios  and  Natural  Frequencies  for  14  DOF  Truss 


undamped 

frequency 

(rad/s) 

min 

iihiif 

min  trace  P 

min 

IIPII2 

c 

mssm 

c 

lER 

c 

BEISI 

319.3471 

0.1249 

315.7222 

0.1077 

316.8960 

0.0712 

316.9671 

290.2645 

281.1909 

0.0763 

283.5066 

0.0534 

288.0130 

281.1347 

279.7743 

0.0620 

279.4916 

0.0729 

275.2599 

271.9703 

265.4633 

0.0406 

267.0294 

0.0377 

267.9009 

264.1012 

258.9797 

0.0297 

260.8798 

0.0327 

262.1410 

236.3806 

0.1376 

245.6173 

0.0965 

243.8296 

0.0857 

243.0665 

223.1620 

224.8248 

0.0728 

224.4934 

0.0693 

224.1446 

165.0056 

0.1353 

167.7850 

0.1166 

167.0317 

EBI 

166.6382 

150.4895 

153.7449 

0.0862 

152.4666 

0.0811 

152.1203 

86.1426 

86.4941 

0.0792 

86.4090 

0.0682 

86.2700 

64.8809 

RRiH 

65.2589 

0.0952 

65.1966 

0.0962 

65.1802 

54.9372 

RRsfl 

55.0353 

0.0331 

55.0221 

0.0337 

55.0725 

31.9937 

32.2309 

0.1087 

32.2221 

nm 

32.2677 

7,7196 

RWl 

7.7581 

0.0986 

7.7580 

0.0986 

7.7582 

where  Pi  and  Pj  are  two  different  solutions.  If  the  matrix  (Pi  -  Pj)  is  positive  definite, 
then  AJ  >  0  for  all  yo-  Hence,  Pj  would  give  a  lower  value  of  J  than  Pi  for  all  yo.  As 
shown  in  Table  6,  the  min  ||H||f  solution  (Pf)*  the  min  (trace  P)  solution  (Ptr)>  ^d 
the  min  ||P||3  solution  (P2)  give  lower  values  of  J  for  all  yo  than  the  solution  using  only 
active  damping  (Pco);  hence,  for  this  exaunple,  they  axe  all  better  solutions  than  using  only 
active  control.  Table  6  also  shows  that  Pjr  gives  a  lower  value  of  J  than  Pj  in  fourteen 
directions,  and  a  higher  value  of  J  in  six  directions.  Since  most  of  the  positive  eigenvalues 
of  (P2  -  P<r)  are  greater  in  magnitude  than  the  negative  eigenvalues,  it  is  argued  that  P^r 
is  a  better  overall  solution  than  P2.  This  is  consistent  with  the  derivation  of  the  min  (trace 
P)  amd  min  ||P||2  solutions  since  min  (trace  P)  minimizes  the  average  of  the  performance 
index  whereas  min  ||P||2  only  minimizes  the  performance  index  for  the  worst  case.  So  one 
would  expect  that  the  min  trace  P  solution  would  be  a  better  solution  for  the  majority  of 
initial  state  vector  directions. 

Another  method  of  comparing  the  solutions  is  comparing  the  relative  differences  be¬ 
tween  the  solutions  and  the  optimal,  unconstrained  solution  (Eqs  (11),  (33)),  and  (34)). 
Table  7  gives  the  relative  differences  between  the  performance  index  attained  with  the  op¬ 
timal,  unconstrained  solution  and  the  performance  indices  attained  with  the  three  solution 
techniques.  All  three  solution  techniques  resulted  in  increases  in  the  performance  index 
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Table  6.  Eigenvalues  of  Delta  P  for  14  DOF  Truss 


umsm 

%jmjm 

50062.43 

50235.92 

50178.14 

592.16 

778.51 

1092.02 

31627.88 

31654.71 

31349.79 

468.76 

171.34 

-501.75 

24620.13 

24776.08 

24465.41 

-324.75 

111.50 

-608.94 

11273.39 

11227.77 

10975.97 

215.36 

-157.09 

-268.53 

9964.37 

10089.19 

9818.28 

185.94 

-116.79 

-237.21 

9061.48 

9043.97 

8936.96 

138.59 

-69.31 

-177.37 

7178.83 

7266.41 

7214.76 

122.97 

31.17 

-168.83 

6473.13 

6536.19 

6495.66 

-39.70 

-48.81 

-108.45 

6114.26 

6326.30 

6265.98 

87.88 

-43.49 

65.78 

1879.97 

1833.14 

1707.37 

79.01 

-31.22 

-79.37 

751.17 

661.39 

683.74 

52.84 

-16.95 

23.78 

449.52 

416.79 

371.62 

25.30 

-5.84 

-24.24 

60.82 

54.98 

60.69 

-5.85 

2.45 

-0.30 

1.81 

1.60 

1.53 

0.07 

-0.21 

0.17 

0.61 

0.61 

0.61 

-0.03 

-0.06 

-0.07 

0.40 

0.39 

0.35 

-0.03 

0.01 

-0.03 

0.35 

0.35 

0.39 

0.03 

-0.01 

-0.02 

0.03 

0.03 

0.02 

0.01 

0.01 

0.01 

0.06 

0.28 

0.28 

-0.01 

-0.01 

0.01 

0.07 

0.05 

0.23 

0.01 

-0.01 

-0.01 

0.28 

0.24 

0.22 

0.00 

0.00 

-0.01 

0.10 

0.23 

0.20 

0.00 

0.00 

0.00 

0.16 

0.08 

0.06 

0.00 

0.00 

0.00 

0.21 

0.09 

0.11 

0.00 

0.00 

-0.01 

0.14 

0.20 

0.13 

0.00 

0.00 

0.00 

0.24 

0.12 

0.08 

0.00 

0.00 

0.00 

0.24 

0.16 

0.16 

0.00 

0.00 

0.00 

0.12 

0.14 

0.08 

0.00 

0.00 

0.00 
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Table  7.  Example  #3  Relative  Increase  in  Performance  Index 


iQia 

trace(P) 

■mm 

0.1120 

0.1534 

Table  8.  Relative  ChanRes  in  Eigenvalues 


c 

■HI 

0.1169 

7.7733 

0.0185 

0.0185 

0.1490 

32.3715 

0.0398 

0.0391 

0.0714 

55.0899 

0.0373 

0.0377 

0.1486 

65.5923 

0.0476 

0.0530 

0.1415 

87.0082 

0.0537 

0.0739 

0.1097 

151.5400 

0.0146 

0.0290 

0.1237 

166.3686 

0.0145 

0.0232 

0.0791 

224.1759 

0.0069 

0.0065 

0.0841 

238.6824 

0.0619 

■29 

0.0184 

0.0396 

263.7692 

0.0214 

wsn 

0.0092 

0.0595 

271.9458 

0.0296 

0.0260 

0.0262 

IB] 

282.0354 

0.0118 

0.0093 

0.0256 

0.0826 

290.0382 

0.0307 

320.6543 

0.0348 

RwSI 

of  less  than  20%;  minimizing  the  trace  of  the  solution  to  the  algebraic  Riccati  equation 
resulted  in  the  lowest  increase,  11%. 

Table  8  gives  a  measure  of  how  dose  the  attainable  system  is  to  the  optimum  system 
produced  by  the  optimal,  unconstrained  solution.  It  does  this  by  computing  the  eigenvalues 
according  to  the  formula 

|A.(A)-MA,pt)| 

|A.(A„p.)| 

The  damping  ratio  and  natural  frequency  listed  in  the  table  are  those  assodated  with  the 
optimum  system.  Minimizing  the  trace  of  P  is  closer  to  the  optimum  solution  than  the 
other  two  solution  techniques  at  higher  frequendes.  At  lower  frequencies,  it  is  not  dear 
whether  the  Probenius  norm  solution  or  the  minimum  two-norm  solution  is  closer  to  the 
optimum  solution;  both  are  closer  than  the  minimum  trace  solution.  Table  9  lists  the 
changes  in  eigenvalues  in  increasing  order.  Looking  at  this  table,  it  is  clear  that  the  largest 
eigenvalue  change  is  7%  and  occurs  in  the  minimum  two-norm  solution.  Tables  7-9  indicate 
that  the  solution  techniques  are  good  approximations  to  the  the  optimal,  unconstrained 
solution. 
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Table  9.  Example  #3  Relative  Changes  in  Eigenvalues  in  Increasing  Order 


trace(P) 


0.0065 

0.0082 

0.0093 

0.0147 

0.0183 

0.0185 

0.0234 

0.0244 

0.0250 

0.0260 

0.0383 

0.0408 

0.0539 

0.0629 


In  this  example,  the  different  solution  techniques  were  demonstrated  for  a  fourteen 
degree  of  freedom  truss.  A  method  of  comparing  the  different  solutions  was  developed 
that  is  independent  of  system  order  and  initial  conditions.  Before  proceeding  to  the  case 
in  which  the  damping  is  provided  by  viscoelastic  materials,  an  alternate  form  of  the  cost 
functional  will  be  touched  on  in  the  next  section. 


2.3  Alternate  Form  of  The  Performance  Index 

The  cost  functional  of  Eq  (15)  treats  the  passive  damping  forces  as  if  they  were 
similar  to  the  active  damping  forces,  as  do  the  following  chapters  on  viscoelastic  damping. 
This  section  takes  a  short  side  trip  to  address  the  argument  that  passive  damping  is  an 
initial  one-time  cost  item  and  ought  to  be  weighted  as  such.  Hence,  the  weighting  on 
damping  is  accomplished  by  weighting  the  damping  constants  directly: 

J  =  ^c^Sc  -{■  ^  (y^Qy  +  u^Ru)  dt  (69) 

The  equations  of  motion  axe  of  the  same  form  as  Eq  (39): 

y  =  ( A  -I-  BvC^)y  -i-  Bu  (39) 

The  vector  c  in  Eq  (69)  is  a  vector  of  the  passive  damping  parameters  and  the  matrix  C 
again  is  a  dis^onal  matrix  of  the  damping  parameters  (i.e.,  the  elements  of  c).  Regardless 
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of  the  values  of  the  damping  parameters,  the  minimum  of  the  integral  term  is  given  by 
aYo  Pyo  where  the  matrix  P  satisfies 


P(A  +  BvC^)  +  (A  +  P  -  PBR“‘B’’P  +  Q  =  0  (70) 

Hence  the  performance  index  can  be  written  as  the  sum  of  two  terms: 

^  =  ^(yoPyo  +  c’’Sc)  (71) 

To  obtain  solutions  independent  of  initial  conditions,  the  same  approach  used  in 
Section  2.2.2  will  be  used  here.  Since 

j  =  ||y.Py.  +  c’’Sc|  <  i(||y.ri|P||,  +  e’-Sc)  =  i(|lP||,  +  t-'Se) 

the  maximum  value  of  the  performance  index  is  minimized  for  ||yo||  =  1  by  minimizing 

y..„-i(l|P|b  +  c’'Sc)  (72) 

where  P  is  found  by  solving  Eq  (70).  This  is  done  by  iterating  with  respect  to  the  viscous 
damping  parameters  c  in  the  same  manner  as  in  Section  2.2.2. 

To  minimize  the  average  value  of  the  performsmce  index  minimize 

=  trace  P  +  (74) 

2n  2 

where  P  is  again  given  by  Eq  (70).  Here  also,  one  needs  to  iterate  with  respect  to  the 
viscous  damping  parameters  c  in  the  same  manner  as  in  Section  2.2.2. 

For  both  of  the  above  solution  techniques,  the  form  of  the  active  control  is  the  same 
as  with  the  original  performance  index  Eq  (15): 

u=-R-‘B’’Py  (11) 

An  example  demonstrating  both  of  the  above  solution  techniques  follows. 
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Table  10.  Dashpot  Constants  and  Performance  Index  Values  for  Example  #4 


Solution  Technique 

min  J 

min  J„„ 

Cl 

130.1174 

132.6323 

dashpot  C2 

76.3213 

290.2810 

constants  C3 

84.9125 

7.1011 

C4 

63.8805 

33.1698 

C5 

77.5260 

60.7109 

J 

484 

648 

dfnax 

4,587 

4,180 

2.3.1  Example  #4  -  Fourteen  DOF  wiUi  Alternate  Version  of  the  Performance 
Index.  This  example  uses  the  same  structure  as  Example  #3,  the  aluminum  truss 
shown  in  Figure  4.  The  same  weighting  matrix  for  control  was  used.  The  magnitude  of 
the  weighting  matrix  S  for  the  damping  coefficients  is  the  same  as  in  Example  #3,  but  the 
dimensions  are  different  since  the  damping  coefficients  are  being  weighted  directly.  The 
two  solutions  cannot  be  compared  in  the  same  manner  as  the  solutions  in  Elxample  #3  due 
to  the  form  of  the  performance  index  (refer  to  Eq  (71)),  so  the  values  of  J  and  will 
be  compared,  as  was  done  in  Example  #2. 

The  solution  techniques  were  run  with  modified  versions  of  the  MATLAB  routines 
used  in  Example  #3.  An  initial  guess  of  c  s  0  was  used  since  a  closed  form  solution  was 
not  obtmned  for  this  form  of  the  performance  index  (Eq  (69)). 

The  resulting  damping  parameters  for  the  min  J  and  min  solutions  are  given 
in  Table  10.  The  dashpot  constants  are  generally  much  greater  in  magnitude  than  those 
found  in  Example  #3,  even  though  the  magnitudes  of  the  weighting  matrices  were  the 
same.  This  is  due  to  the  difference  in  the  way  damping  was  weighted  and  the  different 
dimensions  on  the  damping  weighting  matrix.  Increasing  the  wdghting  on  the  damping 
parameters  would  decrease  the  magnitude  of  the  dashpot  constants  in  Table  10. 

The  active  control  gains  were  determined  using  Eq  (11),  where  the  matrix  P  is  the 
solution  to  the  algebraic  Riccati  equation  ^ven  by  (70).  The  active  control  gains  for 
the  two  solution  techniques  are  listed  in  Appendix  F. 

Choosing  the  dashpot  constants  and  active  control  gains  such  that  the  average  value 
of  the  performance  index  was  minimized  (min  J)  gave  a  23%  reduction  in  the  value  of  the 
performance  index  over  choosing  the  dashpot  constants  and  active  control  gains  such  that 
the  maximum  value  of  the  performance  index  was  minimized  (min  Jmai)-  The  solution  for 
the  minimization  of  the  maximum  value  of  the  performance  index  gave  a  12%  reduction  in 


33 


Table  11.  Damping  Ratios  and  Natural  Frequencies  for  Example  #4 


undamped 

min  J 

min 

Jmax 

frequency 

(rad/s) 

c 

(rad/s) 

c 

(rad/s) 

319.3471 

0.2747 

305.9835 

1.0000 

1.0000 

372.6796 

144.3131 

290.2645 

0.0421 

264.6086 

0.0342 

297.7766 

281.1347 

0.0911 

261.0290 

0.0322 

266.2804 

271.9703 

0.3560 

258.1865 

0.0790 

265.7667 

264.1012 

0.0236 

256.3547 

0.1430 

256.2697 

236.3806 

0.2604 

245.0388 

0.0500 

254.9518 

223.1620 

0.1130 

239.0424 

0.0683 

236.8992 

165.0056 

0.2803 

178.0419 

0.1579 

181.8822 

150.4895 

0.1637 

157.9386 

0.0934 

153.5848 

86.1426 

0.1927 

88.3404 

0.1719 

92.5748 

64.8809 

0.1588 

66.3746 

0.1125 

69.4808 

54.9372 

0.0420 

55.3605 

0.0418 

55.2103 

31.9937 

0.1202 

32.5571 

0.1187 

32.4441 

7.7196 

0.0986 

7.7597 

0.0988 

7.7609 

Jmm€  over  the  aolotion  for  the  minimization  of  the  average  value  of  the  performance  index. 
Therefore,  using  the  min  7  solution  produces  only  a  12%  increase  in  the  maximum  value 
of  the  performance  index,  but  using  the  min  Jmas  solution  increases  the  average  value  of 
the  performance  index  23%. 

Shown  in  Table  11  are  the  open  loop  frequencies  (no  active  feedback  control  or  passive 
damping)  and  the  closed  loop  damping  ratios  and  natural  frequencies  for  the  minimization 
of  the  average  value  of  the  performance  index  (min  J)  and  the  minimization  of  the  maxi¬ 
mum  value  of  the  performance  index  (min  Jmax)-  The  (min  Jmax)  solution  has  two  poles 
on  the  negative  real  axis  (^  =  1).  However,  the  other  thirteen  modes  all  have  damping 
ratios  less  than  0.2  while  the  (min  J)  solution  has  four  modes  with  damping  ratios  greater 
than  0.2. 

In  this  section  two  methods  of  obtaining  solutions  to  the  case  in  which  the  dashpot 
constants  are  weighted  as  a  one-time  cost  have  been  derived.  These  solution  techniques 
are  similar  to  those  developed  in  Section  2.2.2  in  that  one  minimizes  the  maximum  value  of 
the  performance  index  (Eq  (72)),  and  the  other  minimizes  its  average  value  (Elq  (74)).  An 
example  has  been  given  to  illustrate  the  methods  and  how  they  differ  from  the  methods 
presented  in  the  earlier  portion  of  the  chapter. 
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2,4  Summary 

This  chapter  presented  several  new  design  techniques  to  determine  optimal  blending 
of  passive  viscous  damping  and  active  control.  The  techniques  are  based  on  modified 
versions  of  the  standard  linear  quadratic  regulator  performance  index  of  optimal  control 
theory.  It  was  shown  that  it  is  not  possible  to  find  the  optimum  by  solving  a  single  algebraic 
Riccati  equation  as  in  the  standard  linear  quadratic  regulator.  However,  several  iterative 
techniques  involving  a  Riccati  equation  were  developed. 

The  first  version  of  the  proposed  performance  index  treated  passive  damping  as  a  sep¬ 
arate  control  force,  which  resulted  in  an  additional  control  energy  term  in  the  performance 
index.  This  formulation  led  to  two  solutions  (one  of  which  is  closed  form)  that  attempt  to 
minimize  the  error  in  the  performance  index  due  to  using  viscous  dampers  rather  than  a 
truly  active  control  force.  The  closed  form  solution  gave  a  reasonable  ‘iirst-cut’’  design, 
especially  for  larger  systems.  It  is  useful  when  a  “quick  and  dirty”  design  of  active  control 
and  passive  damping  is  needed.  The  closed  form  solution  can  also  be  used  as  a  starting 
estimate  for  the  damping  coefficients  in  two  additional  iterative  techniques. 

One  of  the  additional  iterative  techniques  minimized  the  maximum  value  of  the 
proposed  performance  index  with  respect  to  all  initial  conditions  represented  by  the  unit 
ball  ||yo||  =  1.  If  it  is  pertinent  that  the  performance  index  remmn  below  a  certain  level, 
this  is  the  method  the  design  engineer  is  likely  to  prefer.  Starting  the  iteration  at  the  closed 
form  solution  saved  approximately  40  CPU  minutes  in  computation  time  over  starting  at 
zero  for  the  28th  order  example  problem  ^ven  in  Section  2.2.4.3. 

The  other  iterative  technique  minimized  the  average  value  of  the  proposed  perfor¬ 
mance  index  with  respect  to  all  initial  conditions  represented  by  the  unit  ball  |lyo||  =  1. 
This  technique  gives  good  overall  response,  and  is  more  likely  to  give  better  performance 
over  the  lifetime  of  the  system  since  the  avers^e  of  the  performance  index  is  minimized 
for  all  initial  conditions  represented  by  the  unit  ball  ||yo||  =  1.  Using  the  closed  form 
solution  as  an  estimate  to  start  the  iteration  saved  about  30  CPU  minutes  over  starting 
at  zero.  This  method  ran  about  16  CPU  minutes  faster  than  the  method  that  minimized 
the  maximum  value  of  the  proposed  performance  index.  Both  techniques  yielded  similar 
performance  relative  to  each  other,  and  much  improved  performance  relative  to  a  system 
with  no  passive  damping. 

The  second  version  of  the  performance  index  treats  passive  damping  as  a  one-time 
initial  cost.  Although  a  closed  form  solution  was  not  developed  for  this  case,  iterative  tech¬ 
niques  for  minimizing  the  average  value  and  the  maximum  value  of  this  performance  index 
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were  developed.  An  example  problem  illustrated  the  two  solution  techniques.  Minimizing 
the  average  value  of  the  performance  index  gave  higher  damping  ratios  for  the  majority 
of  the  oscillatory  modes  than  minimizing  the  maximum  value  of  the  performance  index. 
Hence,  minimizing  the  average  value  of  the  performance  index  led  to  shorter  settling  times 
for  those  modes.  So  the  overall  structure  would  tend  to  stop  vibrating  sooner. 

Although  the  development  of  the  second  version  of  the  performance  index  was  only 
carried  out  for  the  case  of  viscous  damping,  the  extension  to  structures  with  viscoelastic 
damping  is  straightforward,  if  the  development  for  the  first  version  of  the  performance  index 
(where  damping  is  weighted  in  the  same  manner  as  active  control)  has  been  carried  out. 
Hence,  the  following  two  chapters  on  viscoelastic  damping  only  present  the  development 
for  the  first  version  of  the  performance  index. 
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III.  Viscoelastic  Damping  Using  a  Classical  Model 
3. 1  Viscoelastic  Materials 

Before  proceeding  to  the  problem  of  simultaneous  optimization  of  active  vibration 
control  and  passive  viscoelastic  damping,  it  is  appropriate  to  discuss  the  characteristics 
of  viscoelastic  materials  and  how  the  behavior  of  viscoelastic  materials  is  modeled  in  this 
chapter. 

The  basic  characteristic  of  viscoelastic  materials  at  constant  uniform  temperature  is 
that  the  modulus  varies  with  frequency,  which  leads  to  three  behavioral  regions  for  the  ma¬ 
terial  [3].  In  the  rubbery  region,  which  occurs  at  low  frequencies,  the  viscoelastic  material 
dissipates  energy  well,  which  is  why  viscoelastic  material  is  good  for  damping  out  vibra¬ 
tions  at  low  frequencies.  In  this  region,  the  real  part  of  the  modulus  stays  fairly  constant 
while  the  imaginary  part  increases  with  increasing  frequency.  The  glassy  region  occurs  at 
high  frequencies.  The  real  part  of  the  modulus  is  again  fairly  constant,  but  the  imaginary 
part  decreases  with  increasing  frequency.  In  this  region  it  is  necessary  to  rely  mainly  on 
active  control  for  vibration  suppression.  The  re^on  between  the  rubbery  region  and  the 
glassy  region  is  called  the  transition  region.  At  these  intermediate  frequencies,  both  the 
real  and  imaginary  parts  of  the  modulus  increase  with  frequency,  with  the  rate  of  increase 
of  the  real  part  slowly  overtaking  that  of  the  imaginary  part.  The  delineation  between  the 
regions  is  sometimes  unclear,  and  is  highly  dependent  on  the  type  of  viscoelastic  material 
being  used. 

The  problems  under  consideration  in  this  dissertation  are  those  in  which  the  vis¬ 
coelastic  dampers  have  been  constructed  so  that  one  component  of  strmn  dissipates  the 
energy.  Viscoelastic  damping  pads  and  constrained  layer  damping  are  examples  of  this 
type  of  damping  treatment.  Damping  treatments  that  typically  don’t  fall  in  this  category 
are  viscoelastic  tuned-mass  dampers  and  compression-type  dampers. 

Since  one  component  of  strain  dissipates  the  energy,  the  mechanical  properties  of  the 
damper  can  be  derived  from  a  scalar  relationship  between  stress  smd  strain  in  the  material. 
A  standsurd  linear  viscoelastic  model  relating  stress  o{t)  and  strain  f(t)  is  [8:14] 


dt* 


(75) 


where  the  values  of  the  constants  bt,  Eq,  and  Ei  depend  on  what  viscoelastic  material  is 
used.  Usually  the  number  of  derivatives  on  strain  is  equal  to  the  number  of  derivatives  on 
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stress  {j  —  k).  The  next  section  assumes  the  more  general  case  of  the  number  of  derivatives 
on  strain  being  less  than  or  equal  to  the  number  of  derivatives  on  stress  {j  <  k). 


3.2  Simultaneous  Optimization  of  Viscoelastic  Damping  and  Active  Control 

This  section  is  analogous  to  Section  2.2  of  Chapter  II,  in  which  the  equations  that 
minimize  a  chosen  performance  index  were  derived  and  an  approximate  solution  was  pro¬ 
posed,  since  the  solution  that  satisfies  the  derived  equations  was  not  obtainable.  The  first 
task  in  this  section  is  to  derive  the  first  order  form  of  the  equations  of  motion.  Next  a 
quadratic  performance  index  similar  to  the  one  used  in  Section  2.2  is  proposed.  Along 
the  way,  expressions  for  the  active  control  and  a  representative  damping  term  are  given  in 
terms  of  the  state.  Finally,  the  equations  which  minimize  the  performance  index  will  be 
derived,  and  an  approximate  solution  presented. 

In  deriving  the  first  order  form  of  the  equations  of  motion  the  viscoelastic  material 
and  the  structure  are  assumed  undisturbed  for  t  <  0.  At  t  =  0,  it  is  assumed  that  a 
disturbance  wq  occurs.  Using  the  classical  model  £q  (75)  for  viscoelastic  damping,  the 
equations  of  motion  of  the  structure  in  the  Laplace  domain  are 

s’Mx(s)  -h  G(s)KrX(s)  -I-  Kx(s)  -I-  Bu(s)  +  Wo(s)  =  0  (76) 


where 

Gq  +  J0i=i 

The  use  of  Gj  instead  of  indicates  that  the  stresses  are  shear  stresses  only. 

The  matrix  Kv  is  analogous  to  the  damping  matrix  D  of  Section  2.2.  Like  the  matrix 
D,  the  matrix  Ky  is  symmetric  and  can  be  written  as  byCby^  where  C  is  a  diagonal  matrix 
of  the  unknown  damping  coefficients. 

Defining  6o  =  1  and  multiplying  Eq  (76)  by  the  summation  Yli=o^i^*  clears  the 
denominator  in  the  viscoelastic  modulus  G(s): 


M  ^  5is''''^x(s)  -f  b. 


i=0 


c£Gis‘b/x(s) 


»=0 


+  K'£^biS'x{s) 


i=0 


-I-  B^6,s*u(s)  +  ^M‘wo(s)  =  0 


tsO 


t=0 


(78) 
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With  some  manipulation,  Eq  (78)  can  be  written  in  the  time  domain  as  a  first  order  system. 
This  will  be  helpful  in  applying  the  approach  used  in  Chapter  II.  But  first,  the  time  domain 
representation  of  the  viscoelastic  forces  and  the  disturbance  will  be  addressed. 

The  Laplace  transform  of  the  damping  forces  is  given  by  the  second  term  in  Eiq  (76), 
G(s)KvX.  For  simplicity,  the  damping  will  be  represented  by  the  inverse  Laplace  transform 
of  the  term  in  brackets  in  Eq  (78): 

v(0  =  Cb/^G~x(/)  (79) 

t=0 


A  new  disturbance  is  defined  as  the  inverse  Laplace  transform  of  the  disturbance 
term  in  Eq  (78): 

w(t)  =  I"‘  |^6.s‘wo(s)|  (80) 

Making  use  of  Eiqs  (79)  and  (80),  Eiq  (78)  is  written  in  the  time  domain  with  the 
highest  derivative  on  the  position  written  separately  from  the  summation: 


d*+2 


»+i 


d* 


M  E  ‘.-’^*(0  +  •>'''(') + K  E  + B  E  + w(i)  =  0 


(81) 

The  next  step  is  to  express  the  highest  derivative  on  the  position  as  a  function  of  the  other 
terms  in  Eq  (81)  (the  dependence  on  time  t  has  been  dropped  for  convenience): 


Ok 


Mg6,.,^x  +  KE‘.^* 

k 


-  ^M-'bvV-  ^M-‘w 

f'-  ^  dV 


t=0 


bk 


This  relation  can  also  be  written  as 


d*+2 

dt*+2* 


if  *  d*  d*+‘ 

M-‘Kx  +  M-‘K6ix  +  (bi.2  +  M-'Kh.)  — x  +  bk-i^x 

-fM-‘BV6i^u-^M-‘bvV-^M-‘w  (82) 

bk  df  bk  bk 


To  write  this  equation  as  a  first  order  system  will  require  the  construction  of  a  higher  order 
state  vector  y.  While  the  state  vector  in  the  previous  chapter  contained  only  position  and 
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velocity,  this  state  vector  will  include  derivatives  on  position  up  to  and  including  order 
k  -{■  1.  By  letting 


y  =  [x’’  x^  •••  d‘+‘xVrf<‘+‘ r 

ai  =  -4^M-‘K  t  =  0,l 

Hi  = -i^[6i-aH-M-‘K6<]  i  =  2,...,k 


Eq  (82)  can  be  written  as  the  following  first  order  system; 


y  = 


0 

I 

0 

1 

...  0 

0 

0 

I 

•  • 

0 

; 

•. 

0 

y  + 

• 

0 

0 

... 

0  1 

0 

-j^M-‘B  _ 

ao 

aa 

aife+i 

k  ^ 


0 

0 

+ 

v-b 

• 

0 

0 

*  d* 

=  Ay  +  Bi^6,^u  +  BvV  +  BwW 


Now  the  disturbance  Wo(t)  is  assumed  to  be  relatively  short-lived  compared  to  the 
response  of  the  structure.  Assume  the  disturbance  stops  at  to  >  0.  The  result  of  the 
disturbance  is 

y(<o)  =  yo  (5) 

If  the  disturbance  is  known,  yo  can  be  calculated,  but  since  solutions  independent  of  the 
initial  conditions  yo  are  sought,  the  actual  value  of  yo  is  of  little  or  no  importance. 
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Hence  the  equations  of  motion  will  be  expressed  as 


y  =  Ay +  B,5^6j^u  +  BvV 
y(«o)  =  yo 


(83) 

(5) 


Although  the  equations  of  motion  are  now  in  first  order  form,  the  derivatives  on 
the  active  control  will  cause  a  problem  if  full  state  feedback  is  used.  If  the  active  control 
is  allowed  to  have  full  state  feedback,  then  the  active  control  term  in  Eq  (82)  will  have 
derivatives  on  position  that  are  higher  than  that  on  the  left  hand  side  of  Eq  (82).  Since 
interest  lies  primarily  in  active  feedback  control  that  consists  of  the  constant  gain  feedback 
of  position  x  and  velocity  x,  the  active  control  u  will  be  so  constrained: 


u  =  G,[  x’’  x’’  I’’  =  Grl/y  (84) 

I  0  0  •••  0 
0  I  0  •  ••  0 

The  matrix  G,  is  the  matrix  of  feedback  gains  that  will  be  determined  by  simultaneous 
design.  (Note  that  G,  is  a  matrix  and  is  not  related  to  the  modulus  coefficients  G,-.) 

It  will  be  convenient  to  express  the  damping  term  v  in  terms  of  the  state  also.  The 
resulting  form  of  v  is  identical  to  the  form  of  the  damping  forces  in  Chapter  II: 


where  = 


v  =  c[^Gob/  Gib/  •••  Gyb/  0  •••  0 

V  =  C^y 


y 


(29) 


To  simultaneously  design  the  active  feedback  control  gains  and  the  damping  coeffi¬ 
cients,  consider  the  performance  index  used  in  the  previous  chapter  in  which  the  passive 
damping  term  is  weighted  as  well  as  displacement,  velocity,  and  active  control: 


J  = 


X 

X 


>  -|-  u^Ru  -f-  v^Sv  I  dt 


(85) 
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To  express  the  performance  index  in  terms  of  the  state  vector  y,  a  new  weighting  matrix 
Q  is  defined  such  that 


y"^Qy  = 


Q. 


X 

X 


(86) 


Since  the  state  vector  y  contains  derivatives  on  position  x  up  through  the  fc  +  1  derivative, 
the  matrix  Q  is  semi-definite;  it  will  have  zeros  on  the  diagonal  corresponding  to  all 
derivatives  higher  than  the  first  derivative  on  x.  The  performance  index  can  now  be 
expressed  in  terms  of  the  state  y: 


J  =  /  ^(y^Qy  +  u’’Ru  -t-  v’’Sv)df 

Jta  2 


(87) 


Now  that  the  performance  index,  the  equations  of  motion,  the  active  control,  and  the 
passive  damping  term  are  all  expressed  in  terms  of  the  state  y,  the  next  step  is  solving  the 
simultaneous  design  problem  of  active  vibration  control  and  passive  viscoelastic  damping. 
The  goal  is  to  determine  the  damping  coefficients  (positive  semi-definite,  diagonal  C)  and 
the  active  control  feedback  gains  (6,)  that 


1 

minimize:  J  =  /  -(y^Qy  +  u^Ru  -f  v^Sy)dt 

Jig  2 

(87) 

*  cf 

subject  to:  y  =  Ay  -1-  Bi 

i=0  *** 

(83) 

u  =  G,l/y 

(84) 

V  =  C#y 

(29) 

The  approach  used  in  Section  2.2  will  be  followed  as  closely  as  possible  in  this  section. 
Using  the  method  of  Lagrange  multipliers  to  append  the  constraints  to  the  performance 
index  gives 


J  =  I  ky^Qy  +  u^’Ru  -I-  v^Sv)  -  A]’(y  -  Ay  -  Bi 

Jia  [2  ,=0 

-  A[(u  -  G,l/y)  -  A>  -  C^y)  I 


dt 


Now  the  vectors  y,u,  and  v  are  taken  to  be  independent.  So  Eqs  (29)  and  (84)  imply 
that  the  matrices  G,  and  C  are  independent  of  these  quantities  also.  By  definition,  the 
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Lagrange  multipliers  Ai,  Aj,  and  A3  are  independent  of  each  other  and  of  y,  u,  v,Gr  and 
C.  From  the  calculus  of  variations,  the  necessary  condition  for  minimizing  the  performance 
index  is  that  its  change  due  to  variations  in  the  independent  variables  vanish  for  arbitrary 
values  of  the  variations  [12:589].  Hence,  to  minimize  the  performance  index  J,  compute 
its  change  due  to  the  variation  of  the  independent  variables: 


=  f  (y’^Q^y  +  u^Rrfu  +  v'^S^v)  -  rfA^(y  -  Ay  -  Bj  -  Bvv) 

•'‘o  L  7^0  « 

-  A[(rfy  -  A6y  -  -  BJv)  -  /5A[(u  -  G,l/y) 

-  XliSu  -  iSGr)lJy  -  G^l/Sy)  -^Aj(v  -  C^y)  -  AJ(/5v  -  ^C^y  -  C^^y)]  dt 

As  in  Chapter  II,  the  final  value  of  the  state  is  assumed  zero  (£q  (6))  and  the  initial 
condition  is  specified  (Eq  (5)),  so  integration  by  parts  yields 

XfSydt  =  ~  X^Sydt 

Integration  by  parts  on  the  derivatives  of  the  active  control  yields 

r  “-w  (^'■)  =  [ 

Then  the  variation  of  the  performance  index  6J  becomes 


=  [(y^’Q  +  aI  +  AIA  +  XlGA/  +  >^c^)6y  +  (u’’R  +  ^  ^A[Bi6.-  -  A[)^u 

+  (v^S  +  A[Bv  -  Xl)Sv  +  A^«G,l/y  +  Aj(^C)^y 

-  «AJ  (y  -  Ay  -  -  Bvv)  -  ^A2(u  -  G,l/y) 

i=0 

-  S>^(y  -  C*y)j  +  g  g(-l)'  BAS  ' 
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Since  the  variations  of  the  dependent  variables  are  independent  of  each  other,  the 
following  equations  need  to  be  satisfied  in  addition  to  the  three  constraint  equations: 


0  —  Ai  -f-  A^A|  Qy 

(10) 

i=0 

(88) 

V  =  -S-‘Bv^A, 

(28) 

\l  =  \l  =  o 

(89) 

The  minimization  of  the  performance  index  is  now  expressed  as  a  system  of  six 
equations.  The  task  is  to  find  the  active  control  gains  and  passive  damping  parameters 
that  satisfy  these  equations.  The  first  step  is  to  eliminate  the  multiple  derivatives  on  Xi  in 
Ek}  (88).  Due  to  the  special  form  of  the  problem,  this  can  be  accomplished  rather  nicely. 

The  higher  derivatives  on  Ai  occur  in  matrix  products  (B^  so  concentration 

will  be  on  these  matrix  products,  with  the  goal  being  to  express  the  summation  of  matrix 
products  in  EJq  (88)  as  a  single  matrix  product.  The  form  of  the  input  matrix  Bi,  the  state 
weighting  matrix  Q,  and  the  plant  matrix  A  will  be  important  in  this  development.  One 
should  remember  that  the  sizes  of  these  matrices  is  determined  by  the  value  of  fc  —  the 
order  of  the  highest  derivative  on  stress  (se '  Sq  (75))  —  and  on  the  Lagrange  multiplier 
Ai  in  Eq  (88).  From  Eq  (10),  the  first  derivative  of  the  Lagrsinge  multiplier  is 

Ai  =  -Qy-A^Ai 

so  the  product  of  the  input  matrix  and  Lagrange  multiplier  is 

BfAi  =  B?’(-Qy  -  A’’ A,)  =  -Bj  A’’ A,  (90) 

since,  for  systems  in  which  there  is  at  least  one  derivative  on  stress  {k  >  1), 


BrQ=o 


This  occurs  because  the  one  non-zero  matrix  element  in  the  composite  matrix  Bi  is  mul¬ 
tiplied  by  a  zero  matrix  element  in  the  composite  matrix  Q  (see  Eqs  (83)  and  (86)).  If  the 
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highest  order  derivative  on  stress  is  second  order  or  greater  {k  >  2),  the  second  term  in 
the  summation  in  Eq  (88)  is  given  by 

BfAi  =  -BrA’’Ax  =  B[A^(Qy  +  A’’Ai)  =  (-l)^Bf(A^)%  (91) 

In  this  case  Bf  A^Q  =  0.  It  follows  easily  that  the  general  form  for  the  matrix  product  in 
the  term  of  the  summation  in  Eq  (88)  is 

=  (-1)‘B?-(A’’)‘A,  for  i  <  k  (92) 

Substituting  Eq  (92)  into  Eq  (88)  produces  an  expression  for  the  active  control  that  does 
not  contain  derivatives  on  Ai : 

u  =  -R-‘Bf  6.( A^)*  j  A,  (93) 

This  can  be  written  more  compactly  by  letting 

u  =  -R-^BflxAi  (94) 

Substituting  Eq  (94)  into  Eq  (84),  and  Eq  (28)  into  Eq  (29)  gives  the  following  relationships 
between  the  sxtive  control  gains,  the  damping  parameters,  and  the  Lagrange  multiplier 
Ai: 


G,l/y  =  -R-‘BrixA,  (95) 

C^y  =  -S"*Bv’’Ai  (96) 

Since  G,  and  C  are  constrsuned  to  be  constant,  solutions  of  the  form  Ai  =  Py  where 
P  =  0  will  be  sought.  With  this  constraint  and  y  not  specified,  Eqs  (95)  and  (96)  lead  to 

G,I/  =  -R-^B^I^P  (97) 

=  -S-*B/P  (98) 

Now  Ir^Ir  =  I,  SO  the  active  control  gains  are  given  by 

G,  =  -R-‘B[I^PI,  (99) 
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The  Lagrange  multiplier  has  been  effectively  eliminated  from  the  derivation  and  the 
active  control  gains  and  damping  parameters  expressed  in  terms  of  the  matrix  P.  Before 
solving  for  P,  one  needs  to  address  the  summation  term  involving  the  higher  derivatives 
on  active  control  u  that  appear  in  Eq  (83).  It  will  be  shown  that  this  summation  can  be 
expressed  as  a  single  term  that  varies  linearly  with  the  state  y. 

From  Eki  (84),  the  first  derivative  on  active  control  can  be  written  in  terms  of  the 


u  =  G, 


=  GJ/Iiy 


(100) 


where 


0  10 


0  0  1 


Ii=  i 


(101) 


Similarly,  the  second  derivative  on  active  control  can  also  be  written  in  terms  of  the  state; 


u  =  G, 


=  G,l/Ii  y 


(102) 


Extrapolating  from  Eqs  (100)  and  (102),  the  derivative  on  active  control  can  be  ex¬ 
pressed  in  terms  of  the  state; 


Hence,  the  sum  of  the  derivatives  on  the  control  u  can  be  expressed  in  terms  of  the  full 
state  vector,  y; 

(fi)’y  =  Gri/  (fi)’!  y  =  G,l/Ii2y  (103) 

»=0  i=0  Li=0  J 


where  the  matrix  Ii?  is  defined  as 


.t=0 
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Substituting  (99)  into  Eq  (103),  and  Ek^s  (28)  and  (103)  into  E^  (83)  gives  an  expression 
for  the  first  derivative  of  the  state  that  is  in  terms  of  the  state  and  the  matrix  P 

y  =  Ay  -  B,R-‘BflxPCl,2y-BvS-‘Bv’^Py  (104) 

where 

Ir  =  I,l/  (105) 

Recalling  that  the  Lagrange  multiplier  can  be  expressed  as  Ai  =  Py,  Ek)  (10)  can  be 
written  in  terms  of  the  state  vector: 

0  =  Py  +  A’’Py  +  Qy  ( 106) 

Substituting  for  the  time  derivative  of  the  state  vector  y  (Eki  (104))  in  Ek}  (106)  produces 
an  equation  solely  in  terms  of  the  matrix  P  and  the  state  y: 

P(Ay  -  BjR-‘Bj'l^Pf4ijy-B^S-‘B/Py)  +  A^’Py  +  Qy  =  0  (107) 

For  general  y,  the  above  equation  leads  to  a  non-symmetric  quadratic  matrix  equation: 

PA  +  A^P  -  PBiR-^Bfl4PfrIia  -  PBvS'^Bv^P  +  Q  =  0  (108) 

This  equation  serves  the  same  role  as  the  algebraic  Riccati  equation  Ekj  (40)  derived  in 
Chapter  II.  However,  due  to  the  asymmetry  of  the  third  term,  it  is  not  a  Riccati  equation. 
Because  the  equation  is  non-synunetric,  its  solution  is  also  non-symmetric.  The  solution 
matrix  P  of  Eki  (108)  can  be  found  using  the  technique  outlined  in  Appendix  C. 

Once  the  solution  matrix  P  is  determined,  the  damping  coefficients  C  can  be  deter¬ 
mined  in  the  same  manner  as  in  Section  2.2.1  —  minimize  the  EYobenius  norm 

1|C^  -I-  S-‘BTP||r  (36) 

where  P  is  the  solution  of  Ekj  (108).  With 

w  =  -s-‘b/p 
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the  elements  of  the  damping  coefficient  matrix  C  are 

Cii  =  max  ;  Cij=0  for  i  ^  j  (37) 

Thus,  a  closed  form  solution  for  the  damping  coefficients  has  been  found.  The  active 
control  gains  are  given  by  Eq  (99). 

The  damping  coefficients  and  active  control  gains  just  determined  are  an  approximate 
solution  to  the  optimization  problem  posed  at  the  be^nning  of  this  section.  The  multiple 
derivatives  on  the  control  in  the  first  order  state  equation  Ek]  (83)  make  the  mathematics 
intractable  when  trying  to  refine  the  choice  of  active  feedback  gains.  Specifically,  the 
ability  to  express  the  higher  derivatives  of  Lagrange  multiplier  Ai  in  terms  of  the  state  was 
due  largely  in  part  to  the  state  weighting  matrix  Q  having  (at  most)  only  two  non-zero 
matrix  elements,  one  weighting  position  x  and  one  weighting  velocity  x.  If  the  damping 
parameter  matrix  C  is  chosen  a  priori,  then  weighting  on  the  state  in  the  performance 
index  Eq  (87)  becomes  Q  -|-  (refer  back  to  Eqs  (15)  and  (38)  of  Chapter  II).  In 

general!  does  not  have  the  nice  form  of  Q;  it  may  contain  non-zero  weightings  for 

the  derivatives  of  the  position  vector  contained  in  the  state  vector  up  to  the  fe‘*  derivative. 
So  Eqs  (99)  and  (37)  represent  the  general  (approximate)  solution  to  the  minimization 
problem  posed  in  this  section. 

In  an  adternative  approach  to  optimizing  viscoelastic  damping  amd  active  control 
simultaneously,  the  following  two  subsections  take  advantage  of  the  fact  that  the  matrix 
Q  -I-  has  no  weighting  on  the  fc  +  1  derivative  of  position  (^i^x). 

3.2.1  Solution  Technique  for  Minimizing  the  Mean  Performance  Index.  If  the 
viscoelastic  material  is  modelled  using  only  one  derivative  on  stress  and  only  one  on  strain, 
two  alternative  methods  of  determining  damping  coefficients  and  active  control  gains  can 
be  developed  that  are  independent  of  initial  conditions.  One  method  consists  of  minimizing 
the  average  value  of  the  cost  functional  over  the  unit  ball  ||yo||  =  1.  The  second  method 
minimizes  the  maximum  value  of  the  cost  functional  over  the  unit  ball  ||yo||  =  1- 

Given  the  active  control  gains  G,  and  the  damping  coefficients  C,  the  first  order 
state  equations  are: 

y  =  (A  -1-  B,G,l7l„-hB,C^)y  (109) 


48 


Then,  using  the  same  procedure  outlined  in  Appendix  B,  the  value  of  the  performance 
index,  £q  (87)  can  be  expressed  in  terms  of  the  initial  state: 

•^  =  ^y2’Piyo  (110) 

The  matrix  Pi  satisfies  the  following  Lyapunov  equation  (refer  to  Appendix  B  for  its 
derivation): 

Pi(A  +  BjG,l/lij+B^C^)  +  (A  +  B,G,l/l,j+BvC*)’^Pi 

+  (Q  +  I,G/RG,I/+*^CSC«)  =  0  (111) 

The  subscript  1  on  Pi  is  to  distinguish  it  from  the  solution  of  Eq  (108).  From  Section  2.2.2, 
The  average  value  of  the  cost  functional  over  the  unit  ball  is  given  by 

Hence,  minimizing  the  trace  of  Pi  minimizes  the  average  value  of  J  over  the  unit  ball 
||yo||  =  It  regardless  of  system  order.  In  this  method,  given  the  damping  coefficients  C, 
one  solves  first  for  the  active  gmn  matrix  G,  using  the  same  procedure  as  in  Section  3.2. 
The  active  gain  matrix  is  of  the  same  form  as  that  in  Section  3.2, 

G,  =  -R'^BflxPIr  (99) 

but  the  matrix  has  a  slightly  different  definition: 

=  I  +  6i(A  +  BvC#)^ 

The  matrix  P  satisfies  a  different  non-symmetric  quadratic  matrix  equation  than  the  one 
given  by  Eq  (108): 

P(A  +  BvC#)  +  (A  +  B^C^)^P  -  PBiR-‘Bn>»PfrIi2  +  (Q  +  ^^CSC^)  =  0  (113) 

Like  the  matrix  1^,  the  matrix  I12  also  has  a  slightly  different  definition  in  Eq  (113)  than 
in  Eq  (108): 

I12  =  I  +  &1I1 
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where  fj  is  defined  by  Eq  (101). 

In  summary,  once  an  initial  guess  for  the  damping  coefficients  C  is  made,  the  matrix 
P  is  found  by  solving  E>}  (113),  the  active  control  gains  are  calculated  via  Eq  (99),  the 
matrix  Pi  is  found  by  solving  Eq  (111),  and  finally  the  average  value  of  the  performance 
index  is  calculated  using  Eki  (112).  Iteration  with  respect  to  C  is  carried  out  until  the 
trace  of  Pi  is  minimized. 

This  subsection  has  derived  a  technique  for  minimizing  the  average  value  of  the  cost 
functional  over  a  representative  sample  of  initial  conditions.  The  next  subsection  will 
present  another  technique  for  optimizing  viscoelastic  damping  and  active  control  simulta¬ 
neously  when  there  is  only  a  first  derivative  on  stress  and  strain. 

S.2.2  Solution  for  Minimizing  the  Maximum  Performance  Index.  A  similar  ap¬ 
proach  to  the  one  just  presented  for  minimizing  the  average  value  of  the  performance  index 
is  used  for  minimizing  the  maximum  value  of  the  performance  index  for  all  ||yo||  =  1.  As 
in  Section  2.2.3,  this  approach  requires  determining  the  two  norm  of  Pi,  since 

J  <  |liP.lll  (U4) 

In  this  method,  iteration  with  respect  to  C  is  carried  out  until  ||Pi||3  is  minimized.  In  each 
iteration  step,  the  matrix  P  is  found  by  solving  EJq  (113)  for  the  current  value  C,  the  active 
control  gains  are  calculated  via  Eq  (99),  the  matrix  Pi  is  found  by  solving  Eq  (111),  and 
finally  the  maximum  value  of  the  performence  index  is  calculated  using  Eq  (114).  Based 
on  this  value,  new  damping  parameters  are  selected  for  the  next  iteration.  This  approach 
might  be  very  conservative  in  general,  since  the  initial  condition  yo  that  maximizes  the 
performance  index  may  not  be  encountered  very  often. 

Examples  follow  which  illustrate  the  three  different  solution  techniques  developed  in 
this  section. 

3.3  Example  Problems 

3.3.0. 1  Example  Problem  #1  -  Single  DOF  System.  In  this  example,  the 
solution  techniques  developed  in  this  section  will  be  applied  to  the  single  degree  of  free¬ 
dom  spring-mass-damper  system  of  Chapter  II.  A  hypothetical  viscoelastic  damper  will  be 
modelled  with  one  derivative  on  stress  and  one  on  strain.  The  parameters  for  the  model  are 
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Go  =  6.8948  X 10^  N/m*  (1.0  Ib/in*) 

Gi  =  4.8264  X 10'*  N-s/m*  (7.0  Ib-sec/in^) 
bi  =  0.001  sec 

The  the  equation  of  motion  in  the  Laplace  domain  is 

1  +  7s 


s^x  + 


1  +  0.001s 


cx  +  0.5x  =  -ti 


(115) 


Converting  to  the  time  domain, 


—X  =  -500x  -  0.5x  -  lOOOx  -  1000(u  +  O.OOlu)  -  1000c(x  +  7x) 


(116) 


The  state  space  equation  is 


y=^ 


• 

X 

0 

1 

0 

0 

X 

.  = 

0 

0 

1 

y  + 

0 

^x 

-500 

-0.5 

-1000 

-1000 

(u  +  O.OOlu) 


0 

0 

-1000 


c[i  7  o]y 


The  weightings  on  the  state,  active  control,  and  passive  viscoelastic  damping  are 


Q  = 


0.5  0  0 
0  1  0 
0  0  0 


,  J?  =  0.1,  5  =  0.1 


(117) 


The  solutions  for  the  three  different  solution  techniques  developed  in  this  chapter  are  given 
in  Table  12.  The  active  gain  matrices  G,  are  included  as  well  as  the  viscoelastic  coeffi¬ 
cients  c.  The  Frobenius  norm  solution  (min  ||H||f)  gave  a  higher  value  of  the  viscoelastic 
coefficients  and  lower  feedback  gains  than  either  of  the  other  two  solutions.  Its  perfor¬ 
mance,  as  measured  by  trace  P  and  IjPllj,  was  not  quite  as  good  as  the  other  solutions, 
but  it  was  within  10%.  It  is  not  surprising  that  the  value  of  the  Frobenius  norm  for  the 
solution  minimizing  trace  P  (and  therefore  the  average  value  of  the  performance  index)  is 
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Table  12.  Example  Problem  #1  Results 


Solution  Technique 

min  llHIlf 

min  trace  P 

min  PII2 

c 

0.3795 

0.1813 

0.1122 

G, 

[  1.3568  2.5293 

[  1.6676  2.8054  ] 

1.7132  3.0486  ] 

trace  P 

1.1924 

1.1256 

1.1374 

0.9692 

0.8902 

0.8830 

0.3096 

0.5406 

0.6731 

U} 

t... 


*13 


Figure  5.  Example  Problem  -  Truss  Strurture 


smaller  than  that  for  the  solution  minimizing  ||P||2  since  the  former  solution  is  closer  to 
the  Frobenius  norm  solution. 

S.S.0.2  Example  Problem  #2  -  Fourteen  DOF  System.  In  this  example,  the 
solution  techniques  developed  in  this  chapter  will  be  applied  to  a  viscoeiastically  damped 
version  of  the  planar  aluminum  truss  used  in  Example  #3  in  Chapter  II.  In  Figure  5, 
the  Uj’s  represent  active  control  forces,  and  the  Cf’s  represent  the  unknown  damping  pa¬ 
rameters  to  be  determined.  The  hypothetical  viscoelastic  material  was  modelled  using 
one  derivative  on  stress  and  one  on  strain.  The  material  is  similar  to  the  one  used  in  the 
previous  example  problem,  but  the  parameters  Gi  and  bi  are  chosen  to  be  ten  times  their 
counterparts  in  the  previous  example.  They  are  chosen  in  this  manner  so  that  the  natural 
frequencies  lie  in  the  transition  region.  The  parameters  for  the  model  are 
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Go  —  6.8948  X 10^  N/m*  (l.Olb/in^) 

Gi  =  4.8264  X 10*  N-s/m’  (70  Ib-sec/in’) 

6i  =  0.01  sec 

The  real  part  (G')  and  imaginary  part  (G")  of  the  modulus  G(s)  is  shown  in  Figure  6. 
The  natural  undamped  frequencies  of  the  truss  lie  in  the  transition  region,  with  most  of 
the  natural  frequencies  near  the  glassy  region.  The  two  verticcd  lines  indicate  the  range  of 
the  natural  frequencies. 

Using  the  development  presented  in  Section  3.2,  the  equations  of  motion  can  be 
written  in  the  state  space  form: 

I  jf 

y  =  Ay  +  Bi'^bi—u  +  B^v  (118) 

i=0 

Since  the  state  vector  contains  position,  velocity,  and  acceleration,  the  dimensions  of  the 
state  vector  y  are  42x1.  There  are  four  control  inputs,  so  the  control  vector  u  is  4x1.  The 
five  viscoelastic  dampers  lead  to  a  damping  parameter  matrix  C  that  is  5x5  and  diagonal. 
Therefore,  A  is  42x42,  B  is  42x4,  Bv  is  42x5,  and  #  is  5x42. 

The  state  weighting  matrix  Q*  was  chosen  such  that  y’^Qy  equalled  the  total  me¬ 
chanical  energy  in  the  system.  Thus  Q*  was  formed  using  the  mass  and  stiffness  matrices. 
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Table  13.  Example  Problem  #2  Results 


min  |1H||f 

min  trace  P 

min  ||P| 

2 

Cl 

3.335 

2.019 

0.372 

C2 

2.443 

1.554 

1.092 

C3 

2.441 

1.605 

0.896 

C4 

3.330 

2.012 

0.902 

Cs 

2.973 

1.857 

0.715 

computation  time 

3  min 

10  hr  54  min 

23  hr  32  min 

C  -  Passive  only 

.0066  to  .1137 

.0042  to  .0875 

.0018  to  .0511 

C  -  Complete  Solution 

.0217  to  .2491 

.0193  to  .2591 

.0097  to  .2591 

M  and  K.  The  weighting  matrices  are 


Q*  = 


K  0 
0  M 


R  =  0.1I,  S  =  0.011 


(119) 


The  solution  techniques  were  run  on  a  VAX  6420  using  MATLAB  routines.  Special 
routines  were  written  that  calculate  the  solutions  to  the  non-symmetric  quadratic  ma¬ 
trix  equations  derived  in  this  chapter,  using  the  algorithm  given  in  Appendix  C.  The 
minimization  algorithm  utilized  is  the  Nelder-Mead  simplex  method  [23]. 

The  values  of  the  damping  parameters  for  the  three  solution  techniques,  along  with 
their  approximate  computation  times,  are  given  in  Table  13.  The  active  control  gains  for 
each  set  of  damping  parameters  are  computed  using  Eq  (99)  and  the  corresponding  matrix 
P.  For  the  min  ||H||i?’  solution,  this  is  the  matrix  that  solves  Eq  (108).  For  the  other  two 
solutions,  it  is  the  matrix  that  satisfies  Eq  (113)  for  the  specific  damping  parameter  matrix 
C.  The  active  control  gains  are  listed  in  Appendix  F.  Because  the  active  control  gains  for 
these  two  solutions  are  calculated  for  a  certain  diagonal  damping  parameter  matrices  (in 
fact,  the  active  gains  and  damping  coefficients  are  essentially  calculated  simultaneously), 
the  damping  ratios  for  the  complete  solutions  are  better  than  those  for  the  min  llH||f 
solution,  whose  active  control  gains  are  based  on  a  non-dici,gonal  parameter  matrix  (the 
damping  coefficients  are  those  that  give  a  diagonal  damping  parameter  matrix  that  best 
approximates  this  optimal  solution).  The  range  of  damping  ratios  with  just  the  viscoelastic 
damping  portion  of  the  solution  implemented  (i.e.,  setting  the  active  gains  to  zero)  is  shown 
in  Table  13,  as  well  as  the  range  of  damping  ratios  for  the  complete  solution. 

How  well  the  complete  solutions  do  in  approximating  the  optimal  solution,  Eqs  (98) 
and  (99),  is  indicated  in  Table  14.  It  gives  a  measure  of  how  close  the  attainable  systems 
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Table  14.  Example  #2  Relative  Changes  in  Eigenvalues 


c 

W„ 

IlHIIr 

trace(P) 

IIPIU 

0.0863 

322.7332 

0.1612 

0.1130 

0.0890 

0.0843 

291.9234 

0.1578 

0.1045 

0.0831 

0.0637 

283.4323 

0.1275 

0.0859 

0.0589 

0.0642 

272.6237 

0.1159 

0.0735 

0.0594 

0.0383 

,264.3360 

0.0830 

0.0584 

0.0401 

0.0943 

241.8979 

0.1034 

0.0952 

0.0857 

0.0932 

226.7390 

0.1245 

0.1039 

0.0850 

0.1814 

172.5763 

0.2119 

0.1725 

0.1548 

0.1689 

156.5123 

0.1940 

0.1585 

0.1408 

0.2862 

93.6310 

0.2286 

0.1989 

0.1986 

0.3002 

70.6018 

0.2496 

0.2131 

0.1985 

0.1585 

55.9026 

0.0985 

0.0966 

0.0343 

0.2859 

34.2871 

0.1373 

0.1181 

0.1031 

0.1930 

7.8979 

0.0783 

0.0781 

0.0780 

1.0000 

94.9204 

0.1546 

0.4417 

0.3909 

1.0000 

109.3135 

0.5831 

0.2399 

0.2751 

1.0000 

107.3431 

0.4012 

0.1400 

0.1261 

1.0000 

105.7911 

0.4360 

0.2859 

0.1649 

1.0000 

102.6897 

0.0999 

0.3035 

0.1397 

1.0000 

101.6959 

0.4445 

0.3626 

0.1822 

1.0000 

99.5807 

0.0343 

0.3158 

0.1842 

1.0000 

99.9061 

0.4822 

0.0338 

0.0368 

1.0000 

100.0000 

0.0000 

0.0000 

0.0000 

1.0000 

100.0000 

0.0000 

0.0000 

0.0000 

1.0000 

100.0000 

0.0000 

0.0000 

0.0000 

1.0000 

100.0000 

0.0000 

0.0000 

0.0000 

1.0000 

100.0000 

0.0000 

0.0000 

0.0000 

1.0000 

100.0000 

0.0000 

0.0000 

0.0000 

are  to  the  optimum  system  produced  by  the  optimal  solution.  The  relative  changes  in  the 
eigenvalues  are  computed  according  to  the  formula 

|A,(A)  -  A,(A„pt)| 

|Aj{Aopt)| 

The  damping  ratio  and  natural  frequency  listed  in  the  table  are  those  associated  with  the 
optimum  system. 

Table  15  lists  the  changes  in  the  optimum  system’s  oscillatory  eigenvalues  in  increas¬ 
ing  order.  The  relative  changes  range  from  4%  (minimum  two-norm  solution)  to  25% 
(Frobenius  norm  solution).  Tables  14  and  15  indicate  that  the  solution  techniques  are 
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Table  15.  Example  #2  Relative  Changes  in  EiRenvalues  in  Increasing  Order 


H 

trace(P) 

l|P||2 

0.0783 

0.0584 

0.0401 

0.0830 

0.0735 

0.0589 

0.0985 

0.0781 

0.0594 

0.1034 

0.0859 

0.0780 

0.1159 

0.0952 

0.0831 

0.1245 

0.0966 

0.0850 

0.1275 

0.1039 

0.0857 

0.1373 

0.1045 

0.0890 

0.1578 

0.1130 

0.0943 

0.1612 

0.1181 

0.1031 

0.1940 

0.1585 

0.1408 

0.2119 

0.1725 

0.1548 

0.2286 

0.1989 

0.1985 

0.2496 

0.2131 

0.1986 

good  approximations  to  the  the  optimal  solution.  The  relative  changes  are  higher  than 
in  Chapter  II,  Example  #3,  but  this  is  to  be  expected  since  the  state  vector  is  larger  (42 
states  versus  28  in  Example  #3)  and  there  is  a  constraint  applied  to  active  control  as  well 
as  to  passive  damping. 

This  section  has  developed  three  solution  techniques  for  a  classically  modelled  vis¬ 
coelastic  system  analogous  to  the  techniques  developed  for  a  viscous  system  in  Chapter  II. 
A  closed  form  solution  was  derived,  as  well  as  two  iterative  solutions,  one  of  which  mini¬ 
mized  the  average  value  of  the  performance  index  over  the  unit  ball,  and  one  that  minimized 
its  maximum  over  the  unit  ball. 


S.4  Multiple  Damping  Materials 

Up  to  now,  it  has  been  assumed  that  only  one  type  of  viscoelastic  material  is  used 
for  damping.  This  section  will  expand  the  results  of  Section  3.2  to  the  case  where  more 
than  one  type  of  viscoelastic  material  is  present  in  the  struciure.  As  stated  in  Section  3.1, 
the  viscoelastic  dampers  are  assumed  to  have  been  constructed  so  that  one  component  of 
strain  dissipates  the  energy,  so  a  scalar  relationship  exists  between  stress  and  strain  in  the 
material. 

Only  two  different  damping  materials  will  be  considered,  but  the  derivation  can  be 
adapted  for  more  materials.  Both  materials  will  be  assumed  to  have  a  frequency  dependent 


56 


modulus  of  the  form 


Go  +  GiS 
1  +  6s 


Although  the  assumption  that  there  is  only  one  derivative  on  stress  has  been  made  here, 
the  derivation  can  be  expanded  for  higher  derivatives  on  stress.  Two  hypothetical  materials 
with  simple  moduli  are  used  so  that  it  is  easier  to  follow  the  derivation. 

The  first  step  to  expressing  the  equations  of  motion  as  a  first  order  system  is  to  write 
the  equation  of  motion  in  the  Laplace  domain: 


s^Mx  +  bvi 


tjQl  +  CjuS 
1  +  biS 


CibJjX  +  b 


'v2' 


V*  f  — 1^  — 


1  -f  hjs 


CjbyjX  +  Kx  +  Bu  =  0 


It  will  be  convenient  to  clear  the  denominators  in  Eq  (120): 


[1  +  (6i  +  62)5  +  6i62S^]s*Mx  -f-  bv}(l  +  b'2s)(Goi  +  G'iis)CibyjX 
+  bv2(l  +  6is)((702  +  Gi2s)C2bva^x  +  K[1  +  (61  +  62)5  +  6i52^*]x 

-j-B[l  -f-  (61  +  62)5  +  6162^*]^  ~  ®  (121) 


As  in  Section  3.2,  the  Laplace  transform  of  the  damping  will  be  represented  by  the 
terms  which  contain  the  damping  psu’ameter  matrices.  For  the  first  damping  material,  the 
Laplace  transform  of  this  term  is 


—  [^01  +  (Goibj  +  <Jii)3  Gii62S*]Cibyix(s) 


(122) 


In  the  time  domain,  this  can  be  written  as 


Vi(t)  —  Cib„i  ^  Goil  (Goifta  +  G^ii)I 


(123) 


Similarly,  the  expression  for  the  representative  damping  due  to  the  second  material  can  be 
written  in  the  time  domain: 


^2(0  —  C2bva^ 


G02I  (G0261  +  Gi2)I  G1261I 


X 
-  X 
X 


(124) 
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Rewriting  the  equation  of  motion  Eq  (121)  in  the  time  domain,  and  then  substituting 
Eqs  (123)  and  (124)  into  it  leads  to  the  following  differential  equation: 

X  +  (6i  +  62)M—  +  Mx  +  bvi  Vi  +  bv2V2  +  Kx 
+  {bi  +  62)Kx  +  bib^Kx  +  Bu  +  (6i  +  62)Bu  +  fri62Bu  =  0  (125) 


Hence  the  highest  derivative  on  x,  ^x  can  be  written  as 


d*  1  ,  -{-  62  , 

^x  =  M-‘Kx  -  44-iM-‘Kx  + 


6162 

6,  +  62 


6162 
1 


M-‘K  +  — ix 
01O2J 


X  — 


6162  dt^  6162 


M  ^B  [u  +  (61  +  62)“  +  5i52u] 


‘bviv,  -  ]~M  *bv2V2 
&1O2  0x^2 


(126) 
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Following  the  procedure  of  Section  3.2,  Eqs  (82)  -  (83),  this  equation  can  be  written  as  a 
first  order  system: 


X 

0 

I 

0 

0 

X 

y  =  < 

X 

0 

0 

I 

0 

• 

X 

^x 

' 

0 

0 

0 

I 

< 

X 

•^x 

_  T 

*i*j  *  J 

^x 

0 


+ 


1 

hiht 


0 

0 

M-‘B 


[u  +  (6i  +  hj)*  +  hiiaii] 


0 

0 

0 

0 

+ 

0 

Vl  + 

0 

For  compactness,  the  matrices  can  be  represented  by  a  single  letter: 

y  =  Ay  +  Bi  [u  +  (6i  +  63)*  +  bib2vi\  +  BviVj  +  B^aVa 
Eqs  (123)  and  (124)  can  be  written  in  terms  of  the  state  y: 

vi(<)  =  Ci^iy 

V2(0  =  Cj^ay 


(127) 


(128) 


(129) 

(130) 


Now  that  the  equations  of  motion  have  been  expressed  in  first  order  form,  a  perfor¬ 
mance  index  is  needed.  Following  the  same  approach  as  in  the  single  damping  material 
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case,  the  performance  index  will  be  defined  as 


X 

X 


+  u’^Ru  -I-  Vj  SiVi+v,  S2V2  dt 


(131) 


As  in  Section  3.2,  the  state  weighting  matrix  Q  is  constructed  such  that 


and  the  matrix  I,^  is  constructed  such  that 


(86) 


u  =  G,[x^  x^V  =  Grl/y  (84) 

Appending  the  constraints  Eqs  (128),  (84),  (129),  and  (130),  to  the  cost  functional  Ekj  (131) 
and  substituting  in  E}q  (86), 


^  [^(y^’Qy  +  u^Ru  +  v[s,v,+v?s,v,) 

-  {y  -  Ay  -  Bju  +  (61  +  62)6  +  hihauJ-ByiVi-BvaVa} 


-  Aj  (u  -  Gri/y)  -  A3(vi  -  Cj^iy)  -  Aj(va  -  Ca^ay) 


dt 


Setting  the  variation  of  J  equal  to  zero,  and  using  the  same  procedures  as  those  used 
in  the  single  material  case  (Section  3.2),  results  in  six  simultaneous  equations: 


y  =  Ay  +  Bi[u  4-  (hi  +  ha)*  +  hihauj+BviVi+BvaV?  (132) 

u  =  G,l/y  (133) 

0  =  Ai  +  A^Ai4-Qy  (134) 

0  =  Ru  +  BiIyiAi  [I^  =  I  +  (hi  +  h2)ii  +  hih2li]  (135) 

Ci#iy  =  -Sj-'Bv/Ai  (136) 

Ca^ay  —  “Sa  ^Bya^Ai  (137) 


Since  Cj  and  C3  are  constant,  Eqs  (136)  and  (137)  imply  that  Ai  is  of  the  form 
Ai  =  Py.  As  in  Section  3.2,  solving  the  above  simultaneous  equations  leads  to  a  non- 
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Figure  7.  Two  Degree  of  Freedom  System 
symmetric  quadratic  matrix  equation: 

PA  +  A^P  -  PBiR-^B[lxPfrIu-PBviSr‘Bv/P  -  PBvaS, 'Bva^'P  +  Q  =  0  (138) 

This  equation  is  solved  in  the  same  manner  as  the  one  for  one  material,  Ekj  (108)  (see 
Appendix  C).  The  damping  parameters  are  found  by  using  the  same  closed  form  solution 
technique  as  in  Section  3.2,  which  is  given  by  Eqs  (36)  and  (37).  The  two  norms  which 
need  to  be  minimized  to  determine  the  damping  coefficients  are 

||Ci*,  +  Sr'Bv,’’P||f  (139) 


and 


IlCa^a  +  Sa-‘Bv3^Pl|f 


(140) 


This  section  derived  simultaneous  design  of  active  vibration  control  and  passive  vis¬ 
coelastic  damping  when  more  than  one  type  of  viscoelastic  material  is  used  in  the  structure. 
An  example  problem  follows  to  demonstrate  the  techniques  derived  in  this  section. 


3.4.1  Example  Problem  #3  -  Two  Viscoelastic  Materials.  To  illustrate  the  case 
of  using  two  different  viscoelastic  materials  for  damping  in  the  same  structure,  the  solution 
technique  developed  in  this  section  will  be  applied  to  the  two  degree  of  freedom  system  in 
Figure  7.  The  viscoelastic  behavior  of  the  hypothetical  dampers  Ci  and  C2  will  be  given  by 


^  1  -t-  .001s 


G2(5)  - 


1  -f-  9s 
1  -f-  .005s 


(141) 


respectively. 
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Following  the  procedure  outlined  in  Section  3.4,  Eqs  (120)  through  (140),  and  using 
the  weightings 

20  -10  0  0 

-10  10  0  0 

,  Si  =  Sa  =  0.01,  R  =  O.ll 

0  0  10 

0  0  0  1 

the  damping  parauneters  are  determined  to  be 

Cl  =  1.75 

Ca  =  1.16 

The  control  gains  were 

2.88  0.02  1.56  0.44 
-0.45  2.03  0.44  1.06 

This  example  has  demonstated  the  feasiblility  of  designing  passive  viscoelastic  damp¬ 
ing  and  active  feedback  control  when  more  than  one  viscoelastic  material  is  used  for  damp¬ 
ing. 

3.5  Summary 

In  this  chapter,  severad  new  design  techniques  to  determine  optimal  blending  of  pais- 
sive  viscoelastic  daunping  and  aM;tive  vibration  control  were  presented.  The  techniques  are 
baised  on  a  standard  lineau'  viscoelastic  model  amd  on  a  modified  version  of  the  standaird 
linear  quadratic  regulator  cost  functional  of  optimal  control  theory.  The  proposed  cost 
functional  treats  passive  damping  ais  a  sepairate  control  force,  which  results  in  an  addi¬ 
tional  energy  term  in  the  cost  functional. 

Two  iterative  techniques  were  developed  in  addition  to  a  closed  form  solution.  The 
closed  form  solution  was  a  least  squares  solution.  One  iterative  technique  minimized  the 
maximum  value  of  the  proposed  cost  functional.  The  other  iterative  technique  minimized 
the  average  value  of  the  proposed  cost  functional.  All  three  techniques  yielded  similar 
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performance  relative  to  each  other,  although  the  closed  form  solution  took  considerably 
less  computation  time. 

The  three  solution  techniques  were  applied  in  Example  #1  to  a  single  degree  of 
freedom  spring-mass-damper  system.  These  techniques  were  confirmed  in  Example  #2  as 
being  useful  ev<'n  when  natural  frequencies  lie  in  the  transition  region  of  the  structure.  The 
applicability  of  the  solution  procedures  were  demonstrated  in  Example  #3  for  a  structure 
in  which  damping  is  provided  via  two  different  viscoelastic  materials. 

One  of  the  drawbacks  of  the  techniques  presented  in  this  chapter  is  that  the  standard 
viscoelastic  model  may  not  adequately  model  the  behavior  of  the  viscoelastic  material. 
Another  drawback  is  that  the  formulation  of  the  problem  could  yield  very  large  state 
space  matrices.  During  an  attempt  to  solve  a  fourteen  degree  of  freedom  problem  with  five 
derivatives  on  stress  strain,  the  large  state  space  matrices  caused  computational  difficulties 
in  trying  to  solve  the  non-symmetric  quadratic  matrix  equation,  Eq  (108).  The  approach 
presented  in  the  following  chapter  overcomes  these  drawbacks  by  using  a  more  accurate 
representation  for  the  behavior  of  the  viscoelastic  material,  and  by  keeping  the  system 
order  down  by  using  the  more  traditional  state  vector  which  consists  of  only  position  and 
velocity. 
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IV.  Viscoelastic  Damping  Using  a  Four  Parameter  Fractional  Derivative  Model 

4.1  Brief  Overview  of  Generalized  Derivatives  as  Applied  to  Viscoelastic  Materials  [11] 

In  this  chapter,  fractional  order  derivatives  will  be  used  to  model  viscoelasticity 
as  opposed  to  the  integer  order  derivatives  used  in  the  previous  chapter.  But  before 
applying  fractional  order  derivatives  to  structural  problems,  it  is  necessary  to  understand 
the  properties  of  generalized  derivatives  and  their  use  in  the  theory  of  viscoelasticity. 
As  will  be  shown,  generalized  derivatives  behave  in  much  the  same  way  as  conventional 
derivatives.  When  used  to  model  viscoelastic  materials,  generalized  derivatives  typically 
provide  an  excellent  model  over  a  broad  range  of  frequencies  [5].  To  show  how  generalized 
derivatives  can  be  used  to  model  viscoelastic  materials,  it  is  appropriate  to  first  present 
the  properties  of  generalized  derivatives,  especially  the  Laplace  and  Fourier  transforms. 
The  generalized  derivative  is  defined  as  [27:59] 

=  forO<a<l  (142) 

This  definition  is  only  valid  for  a  <  1.  However,  the  definition  requires  only  a  slight 
modification  for  a  generalized  derivative  of  order  greater  than  one.  Let  m  be  a  nonnegative 
integer,  and  a  de.4ned  as  before.  Then  [27:59] 

Z?"*-*-“[i^f)]=  /  f  dr  forO<a<l  (143) 

Although  imposing  in  the  time  domain,  in  the  Laplace  (or  Fourier)  domain  the 
generalized  derivative  manifests  itself  as  a  fractional  power  of  s  (or  u;).  To  calculate  the 
Laplace  transform,  a  change  of  variables  will  be  useful: 

T  =  t  —  Tj 

This  leads  to 

=  (144) 

Applying  Leibnitz’s  rule, 

^”W‘»=f(I^i'^l"<‘-’’>''’'+r(r^  forO<a<l  (145) 


64 


Noting  that  the  integral  is  a  convolution  integral,  and  that 


L  f  ]  =  J- 

Lr(l-a)j 


the  Laplace  transform  is 


(146) 


i|i)”l«(0ll  =  (»iWi)l  -  *(o))  +  ^ 


(i-j) 


or,  more  simply, 


Lmxm=s^L[xit)] 


where  the  Laplace  transform  is  defined  as 


Z,[i(t)]  =  /  x{t)e  dt 
Jo 


(148) 


(149) 


For  initial  conditions  equal  to  zero,  the  Laplace  transform  of  a  generalized  derivative  of 
order  a  has  the  same  property  as  the  conventional  derivative:  the  transform  is  s°  times 
the  transform  of  the  function.  In  fact,  the  generalized  derivative  satisfies  many  of  the 
same  properties  as  the  conventional  derivative,  particularly  linearity  and  the  composition 
property  [27:69-84]: 

Z?“[y(<)  +  xit)]  =  i?“[y(0]  +  D°[xit)]  (150) 


D^[D^[x{t)]  =  I?‘*+'’[x(f)] 


(151) 


Linearity  holds  for  a  and  greater  than  or  equal  to  zero.  The  compostion  property  holds 
for  functions  bounded  at  t  =  0. 


The  Fourier  transform  is  defined  as 


F[x(0]=  r  x(<)e-“‘dl 
*'—00 


(152) 


If  x{t)  =  0  for  f  <  0,  then  the  Fourier  transform  can  be  written  as 


F[x{t)]=  r  x{t)e-''^*dt 
Jo 


It  is  easily  seen  that  the  Fourier  transform  of  a  generalized  derivative  is 


(153) 


Fmx{t)]]  =  (.■u;)“F[x(f)] 


(154) 
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In  the  preceding  discussion,  the  only  restriction  placed  on  a  was  that  it  be  a  nonnega¬ 
tive  real  number  less  than  one.  However,  for  en^neering  applications,  an  irrational  number 
can  be  approximated  by  a  rational  number.  So  a  will  now  be  restricted  to  be  rational  as 
well.  Using  the  term  ‘fractional  derivative”  will  indicate  this  additional  restriction. 


The  standard  linear  viscoelastic  model  relating  stress  and  strain  that  was  used  in 
Chapter  III  will  be  useful  in  illustrating  the  use  of  fractional  derivatives  in  viscoelastic 
theory: 


(75) 


Recalling  Scott-Blair’s  proposal  that  fractional  derivatives  could  be  used  to  relate  time- 
dependent  stress  and  strain  in  viscoelastic  materials,  replace  the  conventional  derivatives 
in  Eq  (75)  by  derivatives  of  fractional  order.  The  result  is  the  general  form  of  the  fractional 
derivative  viscoelasi'c  model  [4]: 


a(t)  /)^«l<T(t)l  =  Eo€{t)  -I-  f;„J5“-[6(<)l  (155) 

m=:l  nsl 


A  large  number  of  materials  can  be  modelled  by  replacing  each  sum  in  Eq  (155)  by 
a  single  term  involving  a  fractional  derivative: 


<r(t)  -i-  6I>^[<7(0]  =  Eoc{t)  +  EiD^[e{t)] 


(156) 


Invoking  the  Second  Law  of  Thermodynamics  requires  that  the  parameters  satisfy  the 
following  constraints  [6]: 

Eo  >  0  El  >  bEo 

El  >  0  a  =  P  (157) 

6  >  0 

These  constraints  ensure  nonnegative  energy  dissipation  and  nonnegative  work.  The  stress- 
strain  relation  in  the  Laplace  domain  is 

€(s)  "1-1-55“  ^  ^ 

This  is  known  as  the  four  parameter  model,  and  has  been  shown  to  be  very  accurate  over 
several  decades  of  frequency  [5,  32,  33]. 


66 


Unfortunately,  using  the  same  approach  as  in  the  classically  modelled  viscoelastic 
problem  will  lead  to  the  appearance  of  a  fractional  order  derivative  on  the  control  forces 
as  well  as  on  the  position  vector  [34].  Integration  by  parts  after  taking  the  variation  of 
J  leads  to  an  infinite  sum.  To  find  a  way  around  these  difficulties  for  a  structure  mod¬ 
eled  using  fractional  derivatives,  the  fact  that  the  position  vector  is  a  function  comprised 
of  exponentially  decaying  sinusoids  and  algebraically  decaying  relaxation  modes  w'll  be 
important.  The  relaxation  modes  describe  the  non-oscillatory  return  of  the  structure  to 
its  equilibrium  position.  Hannsgen  and  Wheeler  [16]  point  out  the  fallacy  of  trying  to 
control  relaxation  modes:  decay  rates  are  essentially  the  same  as  those  without  feedback. 
They  note  that  separating  motion  into  creep  modes  and  oscillatory  modes  appears  to  be  of 
fundamental  importance  in  stabilization  and  co  trol  problems  involving  viscoelastic  damp¬ 
ing.  They  also  point  out  that  systems  in  which  viscoelasticity  is  modelled  using  fractional 
derivative  models  are  among  the  systems  that  are  the  least  susceptible  to  destabilization 
due  to  delays  in  boundary  feedback  [15]. 

Hence,  the  design  techniques  presented  in  this  chapter  only  address  the  control  of 
the  exponentially  decaying  modes.  The  error  introduced  by  this  approximation  is  given 
in  Appendix  D,  and  is  bounded.  To  control  only  exponentially  decaying  modes,  one  can 
approximate  the  state  space  equations  by  retaining  only  the  exponential  eigenvalues  and 
eigenvectors.  The  next  section  derives  this  approximation. 


4-2  Simultaneous  Optimization  of  Viscoelastic  Damping  and  Active  Control 

In  this  section,  solution  techniques  for  determining  viscoelastic  coefficients  and  active 
feedback  gmns  that  optimize  performance  indices  like  those  in  Section  2.2.2  are  provided. 
A  key  point  will  be  the  derivation  of  the  approximation  of  the  first  order  state  equations 
in  the  time  domain. 

Before  deriving  the  time  domain  approximation  of  the  first  order  state  equations,  the 
equations  of  motion  will  be  expressed  in  the  Laplace  domain.  The  stress-strain  relation  in 
the  Laplace  domain  is  represented  by  G(s): 


G{s)  = 


Go  +  G,s° 
l  +  bs° 


(159) 


The  use  of  Go  and  Gi  instead  of  Eq  and  Ei  indicates  that  the  stresses  are  shear  stresses. 
The  equations  of  motion  in  the  Laplace  domain  for  a  multiple  degree  of  freedom  system 
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with  viscoelastic  damping  modelled  using  the  four  parameter  fractional  derivative  model 
are  given  by  [4] 

s^Mx(s)  +  G(s)KvX(s)  +  Kx(s)  +  bu(s)  =  0  (160) 

Before  expressing  the  equations  of  motion  in  the  time  domain,  it  is  appropriate  to 
define  the  stress  relaxation  modulus.  The  stress  relaxation  modulus  Grtiit)  gives  the  time 
history  of  the  stress  in  a  viscoelastic  material  which  is  subjected  to  a  unit  step  in  strain  at 
time  <  =  0  [35].  Thus  the  relaxation  modulus  is 

(161) 

As  in  the  previous  two  chapters,  the  matrix  Kv  is  symmetric  and  an  be  expressed  as 

Kv  =  bvCbv’’ 

where  C  is  a  diagonal  matrix  of  the  damping  coefficients  to  be  selected  by  the  design 
process. 

To  control  the  exponentially  decaying  modes  of  the  structure,  an  approximation  to 
the  equations  of  motion  represented  by  Eq  (160)  will  be  formulated  in  the  time  domain 
as  a  first  order  system  that  has  identical  eigenvectors  amd  eigenvalues  as  the  harmonic 
exponential  response  of  Eq  (160).  This  will  require  an  initial  guess  as  to  the  values  of  the 
viscoelastic  damping  coefficients. 

Taking  the  inverse  Laplace  transform  of  Eq  (160),  the  equations  of  motion  can  be 
written  in  the  time  domain 

Mx(t)  +  bv  [cbv^^  f  G„i(t  -  r)x(T)dr]  +  Kx(t)  +  bu(t)  =  0  (162) 
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In  first  order  form,  Eq  (162)  becomes 


M-*K  0 


+  u 

-M-*b 


J  Grei{t  -  r)x(r)d7 


(163) 


Defining  a  state  vector  in  terns  of  position  and  velocity, 


it  is  apparent  that  the  term  in  parentheses  in  Elq  (163)  is  a  time  dependent  vector  function 
of  the  state.  It  will  be  represented  by  an  operator  A  acting  on  time  and  state,  A(t,y). 
Hence,  with  the  input  matrix  defined  as 


-M-*b 


the  equations  of  motion  are  expressed  in  compact,  first  order  form: 


y  =  A(<,y)  +  Bu 


(164) 


The  time  dependent  vector  A(<,y)  can  be  approximated  by  Ay  where  A  is  constructed 
using  the  eigenvectors  and  eigenvalues  associated  with  the  exponentially  decaying  modes 
of  the  open  loop  (u  =  0)  equations  of  motion.  Devereaux  [11]  presents  an  algorithm  which 
calculates  all  the  eigenvectors  and  the  eigenvalues  of  systems  in  the  form  of  Eq  (160).  The 
eigenvectors  and  the  eigenvalues  of  interest  here  are  those  that  lie  in  the  upper  half  of  the 
principal  sheet  of  the  s-plane  (Im  (A,-  >  0)).  The  matrix  'J'  will  represent  the  matrix  of 
eigenvectors  associated  with  exponential  decay. 


^  [  V’l  V>2  •  ■  • 
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where  n  is  the  length  of  the  position  vector  x.  A  diagonal  matrix  of  the  corresponding 
eigenvalues  will  be  represented  by  A: 

0  0 
Aj  •  •  •  0 

0  •  A„ 


The  approximation  of  A(t,y)  will  be  constructed  using  the  matrices  of  eigenvalues  A  and 
eigenvectors  and  their  conjugates,  A  aind  Because  the  approximate  plant  matrix 
A  will  be  constructed  using  eigenvectors,  eigenvalues,  and  their  conjugates,  only  those 
systems  that  are  underdamped  will  be  considered. 

The  next  step  is  to  define  a  diagonad  matrix  A  of  the  composite  eigenvalues  and  their 
conjugates,  and  a  modal  matrix 


A 

0 

9 

9 

A  = 

0 

A 

9  = 

WK 

(165) 


The  approximate  state  matrix  can  now  be  written  in  terms  of  the  eigenvalue  matrix  A  and 
the  modal  matrix  'I': 


A  =  'i'A'i' 


(166) 


An  alternate  but  equivalent  derivation  of  the  approximation  of  A(t,y)  is  given  in  Ap¬ 
pendix  E.  This  alternate  derivation  shows  that  the  approximation  A  is  a  real  matrix,  even 
though  the  matrices  ^  and  A  are  complex. 

Although  an  approximation  for  the  state  equations  has  been  found. 


y  =  Ay-|-Bu  (167) 

an  approximation  of  the  viscoelastic  damping  forces  is  still  needed.  A  performance  index 
similar  to  one  used  in  Chapters  II  and  III,  Elq  (15)  will  be  used  in  this  chapter  also,  so  in 
order  to  weight  the  viscoelastic  forces,  an  appropriate  representation  will  be  found. 

As  in  the  previous  two  chapters,  again  let  v(t)  represent  the  viscoelastic  damping 
forces.  Before  attempting  to  approximate  these  forces,  the  position  vector  and  its  fractional 
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derivative  will  be  expressed  in  terms  of  a  vector  of  modal  coordinates  The  position 
can  be  expressed  as  a  product  of  the  eigenvalue  matrix  and  the  vector  of  modal  coordinates: 

=  (168) 

Taking  the  time  derivatives  of  both  sides  yields  an  expression  for  velocity: 

x(0  =  {t}  =  (169) 

The  generalized  derivative  of  the  modal  coordinate  vector  is  just  the  eigenvalue  ro.atrix 
raised  to  the  order  of  the  deri’/ative  (neglecting  the  algebraic  response): 

D°ri(t)  =  'J'A“t7(0  (170) 

In  particular,  the  first  derivative  of  the  modal  vector  is  the  eigenvalue  matrix  times  the 
modal  vector: 

Ti{t)  =  9A.ri{t)  {V  =  (171) 

From  Eqs  (170)  and  (171),  it  is  apparei.t  that 

D^vit)  =  9A°ri{t)  (172) 

Taking  the  a  order  derivative  of  both  sides  of  Eqs  (168)  and  (169)  and  applying 
Eqs  (170)  and  (172)  yields  expressions  for  the  a-order  derivatives  of  position  and  velocity: 

I>“x(t)  =  (173) 

£>‘*x(t)  =  '9A“ii{t)  (174) 

Recalling  Eq  (171),  and  also  Eq  (169),  the  a-order  derivative  of  position  can  be  written  in 
terms  of  velocity: 

D^xit)  =  ’i'A"  (A-*^(0)  (175) 

= (’J'-^x(t))  (176) 

The  relations  given  by  Eqs  (173),  (174),  and  (176)  will  be  used  in  the  time  domain  formu¬ 
lation  of  viscoelastic  forces. 
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In  the  Laplace  domain,  the  viscoeiastic  forces  are  given  by  the  second  term  in 
Eq  (160): 

i  +  t,.  ‘I'') 

Multiplying  both  sides  by  (1  +  65“*)  and  neglecting  the  algebraic  portion  of  the  fractional 
derivative  once  again,  this  equation  can  be  written  in  the  time  domain: 

(1  +  6Z)“)v(0  =  C(Go  +  G,Z)“)bvM0 

=  GoCbv’’x(0  +  G,Cbv^£>“x(t) 


Substituting  in  Eq  (176)  yields  the  following  expression: 

(1  +  6D“)v(0  =  GoCbv^x(l)  +  GiCbv^’*'A“"*«'-*x 


(177) 


It  is  desired  that  the  time  domain  representation  of  the  viscoelastic  forces  v(t)  be  linear 
in  the  state  y  and  the  matrix  of  viscoelastic  parameters  C,  as  was  the  case  in  Chapters  II 
and  III.  So  the  form  of  v(t)  that  will  be  sought  is  the  same  as  that  in  the  previous  two 
chapters, 

{ 

X 


v(<)  =  { 


=  C^y 


(178) 


where  4^  is  a  matrix  to  be  determined.  The  a-derivative  of  Eq  (178)  wiU  be  shown  to  be 
linear  in  the  state  also.  First  the  a-derivative  of  v(t)  wiU  be  expressed  in  terms  of  the 
modal  vector  rf{t)  by  taking  the  a-derivative  of  Eq  (178)  and  using  Eqs  (173)  and  (174) 
to  express  the  a-derivatives  of  position  and  velocity  in  terns  of  the  modal  vector: 


f  D“x 

= 

[  £>“x 

'*^A“v  J 

The  definition  of  the  modal  vector  is  used  to  bring  position  and  velocity  back  into  the 
expression  for  D°v{t): 

ijrA“4f'-*X 


Z?“v(f)  =  { 


’4'A“^-*x 
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So  the  a-derivative  of  v{t)  is  linear  with  respect  to  the  state  y: 


£)“v(t)  = 


f«l 

0 

= 

0 

["J 

0  'i'A“'i'-* 

Multiplying  this  last  equation  by  the  model  coefficient  b  and  adding  Eq  (178)  gives  an 
alternate  expression  for  the  left  hand  side  of  Eq  ( 177): 


(1  +  6D")v(t)  =  C^y  +  6C^ 


Applying  the  distributive  property, 


0  'S'A"'®'-^ 


(179) 


(1  +  6D")v(<)  =  C# 


1  0 
0  I 


+  b 


0 


0 


(180) 


The  addition  of  the  two  composite  matrices  cam  be  written  as  one  composite  matrix: 


(l  +  6D“)v(t)  =  C# 


0  I  +  b'9A°'^~^ 


(181) 


Setting  the  right  hand  sides  of  Eq  (177)  and  (181)  equal  will  help  in  finding  an  expression 
for  the  matrix 


I  +  0 

0  I  + 


y  =  GoCb/x+  GiCb/'*'A“"*'J'-^x 


GoChJ  GiChJ^A°-^9-^ 


Hence,  for  general  y. 


C# 


0 


=  Cbv^ 


Gol  Gi’J'A"-''*'-' 


(182) 
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Post-multiplying  this  equation  by  the  inverse  of  the  last  matrix  on  the  left-hand  side, 


yields  an  expression  for  the  matrix  product  C^; 


(183) 


=  Cbv^ 


Go  (I  -I-  (I  -I-  6'I'A“'J'-i)"' 


(184) 


If  one  pre-multiplies  this  expression  by  the  inverse  of  C,  the  matrix  ^  can  be  expressed 
explicitly.  However,  it  always  appears  in  the  following  derivation  as  part  of  the  product 
C^.  Also,  since  the  eigenvalue  and  eigenvector  matrices  depend  on  the  value  of  the 
viscoelastic  coefficients,  leaving  ^  expressed  as  part  of  the  product  emphasizes  that 
the  two  cannot  be  separated  in  this  case. 

Now  that  an  expression  has  been  found  for  the  matrix  product  C^,  the  damping 
forces  can  be  weighted  in  the  performance  index: 


J 


(y^Qy  +  u’^R.u  +  v’’sv)dt 


(15) 


Using  the  relation  in  £q  (178)  for  the  damping  forces  v(t)  yields  a  performance  index  like 
the  one  in  Section  2.2.1: 


J=  ^[y^(Q -I- #’’CSC^)y -I- u’’Ru]dt  (38) 

Jto  2 

Defining  a  new  state  weighting  matrix. 


Q  =  (Q  +  ^’’CSC#) 


the  performance  index  has  a  form  identical  to  Eq  (7)  of  Section  2.1: 

/‘/It-  t 
J=  -(y  Qy  +  U  Ru)dt 

Jto  ^ 


(185) 
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Now  that  expressions  have  been  found  for  the  performance  index,  the  equations 
of  motion,  and  the  passive  damping,  the  simultaneous  design  problem  of  active  vibration 
control  and  passive  viscoelastic  damping  can  be  solved.  The  problem  is  to  find  the  damping 
coefficients  and  active  feedback  gains  that 


minimize: 


1  T  »  T 

2(y  Qy+"  Ru)dt 


(185) 


subject  to:  y  =  Ay  +  Bu 
V  =  C^y 


(167) 

(178) 


This  resembles  the  standard  LQR  problem,  which  is  solved  by 


u  =  -R-‘B^Py 


(11) 


and  results  in  a  value  of  the  performance  index  given  by 

J  =  ^yo^Pyo  (41) 

where  P  satisfies 

PA  +  A^P  -  PBR-*B^P  +  Q  =  0 

However,  the  performance  index  J  has  not  yet  been  minimized  with  respect  to  the  damping 
coefficients  since  the  matrix  P  is  a  function  of  the  unknown  damping  coefficient  matrix 
C.  To  minimize  with  respect  to  the  damping  coefficients,  one  can  either  specify  the  initial 
conditions  yo  or  minimize  the  average  of  the  performance  index  J  or  the  maximum  value 
of  the  performance  index  J  over  the  unit  ball,  ||yo||  =  1.  Since  solutions  independent 
of  initial  conditions  are  preferred,  only  the  latter  two  techniques  are  considered  in  this 
chapter.  The  techniques  follow  the  same  algorithms  as  those  in  Sections  2.2.2  and  2.2.3. 
However,  in  computing  Q,  the  matrix  is  computed  using  Eq  (184).  The  matrix  A  is 
computed  using  Eq  (166). 

The  next  section  will  demonstrate  that  the  solution  procedure  is  valid  v '  en  the 
order  of  the  fractional  derivative  is  a  =  1.  That  is,  when  it  reduces  to  the  case  of  viscous 
damping. 
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Figure  8.  Two  DOF  Viscous  System 

4-3  Validation  of  Techniques  for  Viscous  Case 

This  section  is  a  verification  of  the  algorithms  and  computer  programs  that  imple¬ 
ment  the  solution  techniques  presented  in  the  previous  section.  A  FORTRAN  program  that 
is  a  modification  of  the  one  written  by  Devereaux  [11]  computes  the  eigenvalues  and  eigen¬ 
vectors  and  passes  them  to  a  MATLAB  routine  which  makes  the  necessary  calculations  to 
minimize  the  performance  index.  If  the  programs  return  the  same  dashpot  constants  and 
feedback  gains  for  viscous  damping  as  those  calculated  using  the  techniques  developed  in 
Chapter  II,  there  is  a  high  degree  of  confidence  that  the  results  obtained  for  viscoelastic 
damping  will  be  valid.  The  fractional  derivative  model  reduces  to  that  of  a  viscous  damper 
if  the  order  of  the  derivative  on  strain  is  one: 

Go  =  0  0=1 

Gi  =  1  6  =  0 

Although  Gi  is  not  required  to  be  equal  to  one  for  the  four  pau'ameter  model  to  reduce 
to  that  of  a  viscous  damper,  it  is  set  equal  to  one  here  so  that  the  resulting  equations  of 
motion  are  identical  to  those  in  Chapter  II.  In  this  case,  the  technique  described  in  the 
previous  section  ought  to  return  the  same  parameters  as  the  technique  for  viscous  dampers 
given  in  Section  2.2.2.  Damping  parameters  (or  dashpot  constants)  were  computed  using 
both  techniques  for  a  single  degree  of  freedom  system,  a  two  degree  of  freedom  system,  and 
a  fourteen  degree  of  freedom  system.  The  single  degree  of  freedom  system  used  is  the  one 
described  by  Eqs  (64)  and  (65).  The  damping  parameters  obtained  using  the  min  ||P||2 
solution  procedures  (C2)  are  identical,  as  shown  in  Table  16.  The  damping  parameters 
obtained  using  the  min  (trace  P)  solution  procedures  (c,,)  are  essentially  equal  also. 
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The  two  degree  of  freedom  system  is  shown  in  Figure  8.  The  first  order  system  of 
equations  is 


0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

y  = 

-20 

10 

0 

0 

y  + 

-1 

0 

10 

-10 

0 

0 

0 

-1 

«i 

«2 


0  0 

0  0 

-1  -1 

0  1 


Cl 

0 

0 

0  1 

0 

0 

C2 

0 

0  1 

-1 

=  Ay  +  Bu  +  BvC^y 
The  weighting  matrices  were  chosen  to  be 


Q  = 


20 

-10 

0 

0 


-10  0  0 

10  0  0 

0  1  0 

0  0  1 


,  R  =  I,  S=0.5I 


(186) 


As  seen  in  Table  16,  the  damping  parameters  for  both  techniques  were  identical  for  both 
the  min  ||P|l2  and  min  (trace  P)  solution  procedures. 

The  fourteen  degree  of  freedom  system  is  the  same  as  in  Section  2.2.4.3.  Table  16 
shows  that  the  damping  parameters  obtained  for  this  system  for  both  techniques  are  in 
excellent  agreement. 

So  by  applying  both  the  minimization  of  the  average  performance  index  (c,r)  and 
the  minimization  of  the  maximum  value  of  the  performance  index  (C2)  techniques  from 
Sections  2.2.2  and  4.2,  it  has  been  shown  that  the  techniques  of  Section  4.2  are  equivalent 
to  those  of  Section  2.2.2  for  the  case  of  viscous  damping. 
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Igvva?  X  yAOltAgs. 


Cl 

C2 

C3 

C4 

Cs 

SDOF 

vis 

fv 

.8869 

.8872 

Cj 

vis 

fv 

2DOF 

D 

vis 

fv 

HW 

C2 

vis 

fv 

14DOF 

vis 

fv 

54.3826 

54.3825 

35.4502 

35.4503 

38.4443 

38.4458 

42.1846 

42.1846 

46.8419 

46.8414 

C2 

vis 

fv 

76.1842 

76.1846 

17.5549 

17.5585 

28.5870 

28.5835 

34.2785 

34.2784 

45.5261 

45.5242 

vis  =  viscous 

fv  =  frau:tional  derivative  with  a  =  1 


4^4  Transition  from  Viscous  to  Viscoelastic 

The  last  section  verified  that  the  solution  techniques  derived  in  this  chapter  for 
fractional  order  systems  are  valid  for  the  case  of  viscous  damping.  This  section  will  show 
that  the  viscoelastic  damping  parameters  (the  dis^onaJ  elements  of  C)  are  smooth  functions 
of  the  viscoelastic  parameters  Go,  a,  and  6.  Hence,  it  can  be  inferred  that  the  solution 
technique  for  fractional  order  systems  is  valid  for  all  values  of  a  from  a  =  0  to  a  =  1. 

First,  a  will  be  varied  incremently  from  a  =  1  to  q  =  0,  while  keeping  Go  =  0  and 
6  =  0.  This  is  equivalent  to  varying  the  viscoelastic  modulus  G(s)  from  pure  viscous  to 
pure  elastic. 

The  two  degree  of  freedom  system  of  the  previous  section  (Figure  8  and  Eqs  (186) 
and  (186))  will  be  used  in  the  calculations.  Figure  9  starts  with  the  damping  parameters 
(cj)  for  the  pure  viscous  model  (q  =  1,  Go  =  0,  6  =  0); 

G(s)  =  7.3s 

The  reasons  for  choosing  Gi  =  7.3  will  become  apparent  later.  As  a  is  decreased  incre¬ 
mentally  down  to  zero,  Figure  9  shows  the  resulting  damping  coefficients  for  both  solution 
techniques,  min  ||P||2  and  min  trace(P).  Figure  9  shows  that  the  transition  from  pure 
viscous  to  pure  elastic  (a  =  0) 

G(s)  =  7.3 
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is  smooth  within  the  resolution  of  the  data  and  monotonically  decreasing. 

Next,  Go  and  b  are  varied.  The  end  result  will  be  the  model  which  will  be  used  in 
the  example  problem  in  the  following  section.  Hence,  the  fractional  order  a  will  be  set 
to  a  =  4/7,  and  Gq  will  be  increased  from  0  to  1.  The  parameter  b  will  be  increased  to 
0.0008.  Thus,  Figure  10  shows  that  the  transition  from  the  viscoelastic  modulus 


G(s)  =  7.35“/^ 


to  the  viscoelastic  modulus 

1  -h  7.3s^/^ 

1  +  .0008sV7 


(187) 


is  also  smooth  within  the  resolution  of  the  data.  Eq  (187)  is  based  on  a  model  for  3M-467 
adhesive  at  75®  F,  a  viscoelastic  material  used  for  damping  [3:126].  The  constants  Go, 
Gi,  and  b  have  the  same  value  as  the  actual  model,  but  4/7  is  an  approximation  for  the 
order  of  the  fractional  derivatives  on  stress  and  strain.  Reference  [3]  gives  the  order  of  the 
derivatives  on  stress  and  strain  as  .51  and  .56  respectively.  Since  4/7  «  .57,  a  =  4/7  is 
a  good  approximation  for  order  of  the  fractional  derivative  on  strain,  whose  coefficient  is 
four  orders  of  magnitude  greater  than  the  fractional  derivative  on  stress.  Using  a  lower 
order  fraction  simplifies  the  calculation  of  the  eigenvectors  and  eigenvalues. 
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Example  Problem  -  Fourteen  DOF  System 


In  this  section,  the  fractional  order  solution  technique  will  be  applied  to  the  viscoelas- 
tically  damped  planar  aluminum  truss  of  the  previous  chapter  (Figure  5).  The  viscoelastic 
material,  3M-467  adhesive,  was  modeled  with  the  four  parameter  model  given  by  Eq  (187). 
The  modulus  plot  of  G(s)  is  shown  in  Figure  11,  where  the  real  part  is  plotted  as  G'  and 
the  imaginary  part  as  G".  The  vertical  lines  in  the  graph  indicate  the  lowest  and  highest 
frequencies,  and  show  that  the  structure’s  natural  frequencies  lie  well  within  the  transition 
region.  The  state  weighting  matrix  Q  was  chosen  such  that  y^Qy  equals  the  total  me¬ 
chanical  energy  in  the  system,  while  the  weightings  on  active  control  and  passive  damping 
were  chosen  to  be 


R  =  0.1I,  S  =  10-®I 

The  solutions  were  computed  using  the  algorithms  given  in  Section  4.2  and  Sec¬ 
tion  2.2.2.  The  eigenvectors  and  eigenvalues  of  the  system  for  a  given  set  of  damping 
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80 


0.03 


Minimum  |(Pj(2  Solution 


Minimum  J(r  Solution 


Figure  10.  Example  -  Transition  to  Full  Model 
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parameters  C  were  computed  by  a  FORTRAN  program  that  used  the  algorithm  devel¬ 
oped  by  Devereaux  (see  Reference  [11]).  The  eigenvectors  and  eigenvalues  were  passed  to 
a  MATLAB  routine  which  computed  the  approximate  state  matrix,  the  appropriate  gains 
for  that  state  matrix,  and  the  resulting  value  of  ||P||a  or  trace(P),  depending  on  which 
one  was  being  minimized.  A  new  set  of  damping  coefficients  was  then  selected  using  the 
MATLAB  program  FMINS,  and  the  new  eigenvalues  and  eigenvectors  computed.  The 
iteration  continued  until  convergence  occurred. 

The  values  of  the  resulting  damping  parameters  (c,)  for  the  two  solution  procedures, 
min  ||P||2  and  min  (trace  P),  are  shown  in  Table  17.  The  system  damping  ratios  with  and 
without  active  control  are  included  as  well.  The  active  control  gains  were  computed  using 
Eq  (11)  and  are  given  in  Appendix  F. 

Since  the  properties  of  viscoelastic  materials  are  temperature  dependent,  it  is  con¬ 
ceivable  that  a  drastic  change  in  temperature  could  result  in  a  rather  drastic  change  in 
properties.  Assume  this  has  occurred,  and  the  modtilus  G(s)  is  now  given  by 


G(s) 


1  + 

1  -I-  .08s</7 


The  real  part  of  the  modulus  is  plotted  as  G'  in  Figure  12,  while  the  imaginary  part  is 
plotted  as  G".  Recomputing  the  viscoelastic  parameters  for  this  case  produces  the  results 
shown  in  Table  18.  Even  though  the  natural  frequencies  are  now  very  near  the  glassy 
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Table  17.  Fractional  Viscoelastic  with  Gi  =  7.3,  6  =  .0008 


min  trace  P 

min  P  2 

Cl 

4.2782 

3.6768 

Cj 

3.9100 

4.4106 

C3 

4.2057 

5.8892 

C4 

5.2116 

6.0427 

Cs 

6.0598 

6.5551 

C  -  Passive  only 

.0003  to  .0080 

.0004  to  .0092 

C  -  Complete  Solution 

.0041  to  .0340 

.0044  to  .0341 

Figure  12.  Viscoelastic  Modulus  Plot  -  Gi  =  730,  6  =  .08 


Table  18.  Fractional  Viscoelastic  with  Gi  =  730,  b  =  .08 


min  trace  P 

min  1  P  I2 

Cl 

0.0563 

0.0085 

C2 

0.0576 

0.0001 

C3 

0.0722 

0.1592 

C4 

0.0661 

0.0497 

C6 

0.1048 

0.1194 

^  -  Passive  only 

.0002  to  .0029 

.0001  to  .0036 

C  -  Complete  Solution 

.0023  to  .0340 

.0025  to  .0340 
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region,  the  solution  techniques  can  still  be  used.  However,  damping  relies  more  heavily  on 
active  control.  The  active  control  gains  are  given  in  Appendix  F. 

The  results  above  indicate  that  the  viscoelastic  material  must  be  carefully  chosen 
such  that  the  natural  frequencies  of  the  problem  structural  modes  lie  in  the  transition 
region  for  significant  viscoelastic  damping  to  occur.  This  is  not  a  concern  with  viscous 
dampers  because  they  have  no  high  and  low  frequency  plateaus. 

4.6  Summary 

This  chapter  presents  two  new  design  techniques  to  determine  optimal  blending  of 
passive  viscoelastic  damping  and  active  vibration  control.  The  techniques  are  iterative  in 
nature,  and  are  based  on  a  fractional  order  derivative  viscoelastic  model  and  on  modified 
versions  of  the  standard  linear  quadratic  regulator  performance  index  of  optimal  control 
theory,  given  by  Eq  (7). 

The  performance  index,  Eq  (15),  treats  passive  damping  as  a  separate  control  force, 
which  results  in  an  additional  energy  term  in  the  performance  index.  By  neglecting  the 
algebraic  modes  due  to  the  viscoelastic  dampers,  two  of  the  solution  techniques  developed 
in  Chapter  II  can  be  applied  to  the  viscoelastically  damped  structure;  the  minimization  of 
the  average  value  of  the  performance  index  amd  the  minimization  of  the  maximum  value  of 
the  performance  index.  These  solution  techniques  are  useful  when  the  viscoelastic  material 
is  well  modelled  by  a  four  parameter  fractional  order  model. 
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V.  Conclusions  and  Recommendations 


The  new  design  techniques  presented  in  this  dissertation  provide  design  engineers 
with  the  ability  to  determine  complementary  feedback  gains  and  passive  damping  param¬ 
eters  independent  of  initial  conditions  and  the  ability  to  handle  viscoelastic  dampers  when 
determining  the  optimum  blending  of  active  vibration  control  and  passive  structural  damp¬ 
ing.  In  deriving  new  design  techniques  to  determine  optimal  blending  of  passive  structural 
damping  and  active  vibration  control,  the  treatment  of  passive  structural  damping  forces 
in  this  dissertation  is  unique.  Usually  passive  structural  damping  forces  are  just  included 
as  part  of  the  state,  but  here  they  are  weighted  separately  from  the  state  in  the  linear 
quadratic  regulator  (LQR)  performance  index.  The  passive  damping  forces  are  treated  as 
co-equals  to  the  active  control  forces,  and  are  weighted  in  a  similar  manner.  The  modifi¬ 
cation  to  the  performance  index  is  minor,  but  allows  structural  damping  to  be  treated  as 
a  control  parameter  in  LQR  theory.  Viewing  structural  damping  as  an  LQR  parameter 
increases  the  usefulness  of  LQR  theory  by  enabling  the  design  of  complementary  feedback 
gains  and  passive  damping  parameters.  This  reduces  the  active  control  effort  needed  to 
meet  specifications,  allowing  the  design  en^neer  to  take  full  advantage  of  each  ounce  of 
damping  material.  It  also  allows  the  design  en^neer  to  avoid  implementing  excess  damping 
that  inadvertently  increases  response  times. 

In  addition,  the  techniques  devdoped  in  this  dissertation  enable  the  design  engineer 
to  select  complementary  feedback  gains  and  passive  damping  parameters  that  are  indepen¬ 
dent  of  initial  conditions.  Usually  in  a  problem  dealing  with  active  and  passive  control, 
initial  conditions  are  set  before  working  the  problem.  This  restriction  was  not  used  in  this 
dissertation,  so  the  solution  techniques  presented  here  lead  to  designs  that  are  more  robust 
with  respect  to  initial  conditions. 

Another  restriction  relaxed  in  this  dissertation  is  related  to  the  structural  damping. 
In  control  problems,  structural  damping  is  often  approximated  by  viscous  damping,  since 
other  types  of  damping  lead  to  complicated  equations  of  motion.  This  dissertation  develops 
the  techniques  for  viscous  damping,  but  extends  the  results  to  the  more  realistic  case  of 
viscoelastic  damping. 

Three  techniques  were  developed  for  the  viscous  case  —  one  closed  form  and  two 
iterative.  The  closed  form  solution  is  a  least  squares  solution.  One  iterative  technique 
minimizes  the  maximum  value  of  the  proposed  cost  functional  with  respect  to  the  ini¬ 
tial  state,  while  the  other  minimizes  the  average  value  of  the  proposed  cost  functional. 
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Both  techniques  yielded  similar  performance  relative  to  each  other,  and  much  improved 
performance  relative  to  a  system  with  no  passive  damping. 

Three  similar  techniques  were  also  developed  for  viscoelasticity  modelled  with  a  clas¬ 
sic  model  using  integer  order  derivatives  on  stress  and  strain.  The  closed  form  solution 
proved  to  be  much  cheaper  computationally  than  the  iterative  techniques,  while  giving 
similar  performance. 

The  classic  viscoelastic  approach  can  be  rather  difficult  to  implement  for  large  struc¬ 
tures  as  the  order  of  the  state  space  equations  increases  dramatically  with  higher  order 
derivatives  on  stress  and  strain.  The  order  of  the  state  space  equations  was  kept  low  in 
the  approach  developed  for  the  fractional  derivative  viscoelastic  model.  Another  benefit  of 
using  the  fractional  derivative  model  is  that  it  is  often  a  more  accurate  description  of  the 
materials  behavior.  Two  iterative  techniques,  similar  to  those  for  the  viscous  case,  were 
developed  for  this  case.  A  closed  form  solution  was  not  possible  due  to  the  nature  of  the 
problem. 

All  of  the  techniques  derived  in  this  dissertation  are  applicable  to  structures  that  can 
be  modelled  with  finite  elements.  The  techniques  are  useful  in  cases  where  the  baseline 
structure  is  well-defined  and  the  design  en^neer  wishes  to  use  a  combination  of  active 
and  passive  vibration  suppression.  If  the  damping  is  adequately  represented  by  a  viscous 
model,  then  the  corresponding  solution  techniques  should  be  used.  The  Frobenius  norm 
solution  gives  a  quick,  rough-cut  answer,  and  is  adequate  for  many  applications.  If  a  better 
performamce  is  required,  the  minimization  of  the  average  value  of  the  performance  index 
is  recommended.  It  will  tend  to  give  better  response  times  and  takes  less  time  to  com¬ 
pute.  The  minimization  of  the  maximum  value  of  the  performance  index  is  recommended 
only  if  the  design  engineer  is  concerned  that  the  minimization  of  the  average  value  of 
the  performance  index  may  yield  high  values  of  the  performance  index  for  certain  initial 
conditions. 

If  the  damping  is  to  be  modelled  with  a  classical  viscoelastic  model,  the  design 
engineer  will  use  the  Frobenius  norm  solution.  If  the  model  only  has  first  derivatives  on 
stress  and  strain,  the  design  engineer  may  want  to  refine  his  design  by  minimizing  the 
maximum  value  of  the  performance  index.  This  method  is  recommended  over  minimizing 
the  maximum  value  of  the  performance  index  for  the  same  reasons  as  in  the  viscous  case. 

Also  for  the  same  reasons,  minimizing  the  average  value  of  the  performance  index 
is  recommended  for  the  case  in  which  viscoelastic  damping  is  modelled  using  a  fractional 
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derivative  model.  A  four  parameter  fractional  derivative  model  is  recommended  over  a 
classical  model  for  the  reasons  cited  above;  it  keeps  the  order  of  the  state  space  equations 
low  and  it  is  often  a  more  accurate  model. 

For  situations  in  which  the  design  engineer  wishes  to  weight  passive  damping  as  a 
one-time  cost,  two  8<dution  techniques  were  developed  for  structures  with  viscous  damping. 
It  would  be  a  straightforward  task  to  extend  these  solution  techniques  for  viscoelastic 
damping.  As  in  the  case  in  which  passive  damping  is  weighted  as  an  additional  control  force, 
minimizing  the  average  value  of  the  performance  index  is  recommended  over  minimizing 
the  maximum  value  since  it  will  give  better  performance. 

The  weighting  matrices  in  both  types  of  performance  indices  are  chosen  by  the  design 
engineer  based  on  the  relative  importance  of  passive  damping,  active  vibration  control, 
and  required  damping  ratios.  The  weighting  matrices  may  need  to  be  adjusted  to  attain 
a  solution  that  is  acceptable  to  the  design  engineer.  Two  general  guidelines  in  choosing 
the  weighting  matrices  are  that  the  diagonal  elements  of  the  weighting  matrices  for  passive 
damping  and  active  control  be  larger  than  the  diagonal  elements  of  the  weighting  matrix 
for  the  state,  and  that  the  diagonal  elements  of  the  weighting  matrix  for  passive  damping 
be  less  than  or  equal  to  those  for  active  control.  However,  these  guidelines  do  not  have  to 
be  followed  if  the  design  engineer  has  a  reason  for  choosing  the  weightings  differently. 

The  solution  techniques  in  this  dissertation  are  limited  to  structures  in  which  the 
structural  damping  is  either  viscous  or  viscodastic.  To  use  the  solution  technique  for 
viscoelastic  damping,  the  viscoelastic  dampers  must  be  constructed  such  that  only  one 
component  of  strain  dissipates  the  energy.  Expainding  the  techniques  to  different  types  of 
dampers,  such  as  viscoelastic  dampers  in  which  two  components  of  strain  dissipate  energy, 
is  an  area  for  further  research. 

Another  limitation  is  in  the  numerical  implementation  of  the  solution  techniques. 
As  noted  in  the  summary  of  Chapter  III,  if  the  order  of  the  finite  element  model  is  too 
large,  computational  difficulties  may  arise.  This  is  especially  critical  in  the  case  where 
viscoelasticity  is  modelled  using  integer  order  derivatives,  since  the  order  of  the  state 
equations  are  at  least  three  times  larger  than  the  order  of  the  finite  element  matrices 
defining  the  structure.  The  numerical  stability  for  problems  of  more  than  fourteen  degrees 
of  freedom  was  not  investigated. 

Applying  the  solution  techniques  to  larger  scale  problems  would  be  an  area  for  further 
research.  Structures  other  than  the  spring-mass-damper  systems  and  trusses  in  the  example 
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problems  could  also  be  considered.  Other  structures  might  be  basic  structures,  such  as 
plates  and  beams,  and  also  more  complex  structures  that  are  combinations  of  these  simpler 
structures.  Another  area  for  further  work  would  be  to  experimentally  verify  the  solution 
techniques  using  some  of  these  structures;  the  work  presented  in  this  dissertation  is  ail 
theoretical. 

Since  this  is  an  engineering  dissertation,  an  engineering  approach  was  taken  in  mini¬ 
mizing  the  cost  functional  subject  to  the  constraints  on  the  passive  damping  force  (and  on 
the  active  force  when  viscoelastic  damping  was  modelled  with  a  classical  model).  Follow¬ 
ing  an  exact  analytic  approach  using  vector  space  methods  would  be  an  opportunity  for 
further  research. 

The  two  greatest  contributions  of  this  dissertation  are  the  ability  to  determine  com¬ 
plementary  feedback  gains  and  passive  damping  parameters  independent  of  initial  condi¬ 
tions  and  the  ability  to  handle  viscoelastic  dampers.  These  two  contributions  combine 
to  show  that  viscoelastic  damping  can  be  optimized  in  conjunction  with  active  feedback 
control  to  provide  better  vibration  reduction. 
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Appendix  A.  Calculation  of  Damping  Coefficients 

The  derivation  of  the  solution  to  the  diagonal  C  that  minimizes  the  Frobenius  norm 
of  Eq  (36)  is  given  in  this  appendix.  Although  the  development  is  for  vectors  and  matrices 
of  order  4,  it  can  be  easily  extended  to  vectors  and  matrices  of  order  n. 

Let  X  be  the  space  of  real  valued  vectors  with  four  elements  and  let  Y  be  the  space 
of  4x4  real  valued  matrices.  Define  an  operator  A  ;  X  -+  Y  by 

0  0 
0  0 
13  0 

0  X4 

where  x  €  X.  It  is  easily  verified  that  A(xi  +  axj)  =  A(xi)  +  aA(x2)  for  any  scalar  a,  so 
A  is  linear. 

Given  y  €  Y,  the  vector  x  is  sought  that  minimizes  the  Frobenius  norm  ||y  -  AjcJIp. 
Hence,  a  solution  to  the  following  expression  is  sought: 


y  =  A(x)  = 


ii  0 

0  x^ 
0  0 
0  0 


nunljy- Ax|1p 


(188) 


Define  a  functional  on  Y  x  Y  by 

(y.*)  =  £  V  y,z  G  Y  (189) 

tj=i 

Since  this  functional  satisfies  the  four  reqtiirements  of  a  real  inner  product: 

(y,z)  =  (z,y) 

(yi+y2,z)  =  (yi,z)  +  (y2,z) 

(Ay,z)=  A(y,z) 

(y»y)>0,  with  equality  iff  y  =  0 
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define  the  functional  in  Eq  (189)  to  be  the  inner  product  on  Y.  Then  the  norm  on  Y  is 
given  by 

l|y||  =  (y.y)"’  = 

Notice  that  this  is  the  Frobenius  norm  for  4x4  matrices. 

The  conjugate  operator  /I*  :  Y  — ►  X  is  found  from  the  relationship 

(x,>i*y)x  =  (>ix,y)Y 

The  subscripts  X  and  Y  denote  the  spaces  in  which  the  inner  product  is  taken.  Since 

4  4 

(/lx,y)Y  =  =  (**^*y)x 

ij=l  <=1 

the  conjugate  operator  is  given  by 


A*y=< 


yn 
yn 
yss 
[  y«  j 


The  solution  which  yields  the  smallest  error  in  the  Frobenius  norm,  EJq  (188),  is  given  by 
{A* A)~^ A*y  if  (A*  A)“^  exists  [21:160].  Find  A* A: 


/ 

0 

0 

0 

\ 

A*  Ax  =  A* 

0 

l2 

0 

0 

=  < 

l2 

0 

0 

*3 

0 

\ 

0 

0 

0 

X4  , 

} 

X4 

\  ^  J 

So  A*  A  is  the  identity  operator  and  its  inverse  is  also  the  identity  operator.  Hence  the 
solution  that  minimizes  the  Frobenius  norm  is 


X  =  A*y 

In  other  words,  x  is  the  diagonal  of  y. 
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Using  the  same  procedure  outlined  above,  it  can  be  shown  that  the  dis^onal  C  that 
minimizes  the  Probenius  norm  ||C^  —  Wjl/*  is  given  by 

c  =  (AM)-U*W  = 

where  c  is  a  vector  containing  the  m  elements  of  the  diagonal  of  C,  and  #  and  W  are 
m  X  n  matrices.  In  the  case  of  Eq  (35), 

W  =  -S-*Bv^P 

The  diagonal  elements  of  the  positive  semi-deiinite  diagonal  matrix  C  that  minimizes 
the  FYobenius  norm  (|C^  -  W((/i-  are 

(37) 

The  last  sentence  of  the  following  theorem  will  be  used  to  prove  that  C  as  defined 
above  minimizes  the  norm  ||C$  -  Wjlf  over  the  set  of  positive  semi-definite  diagonal 
matrices: 

Theorem  [21:69].  Let  x  be  a  vector  in  a  Hilbert  space  H  and  let  K  be  a  closed  convex 
subset  of  H .  Then  there  is  a  unique  vector  ko  ^  K  such  that 

Ik  -  fcoll  <  Ik  -  *11 

for  all  k  G  K.  Furthermore,  a  necessary  and  sufficient  condition  that  ko  be  the  unique 
minimizing  vector  is  that  (x  —  k(,,k  —  ko)  <  0  for  all  k  E  K. 

Let  H  be  the  space  of  all  m  x  n  matrices.  The  following  functional  satisfies  the  four 
requirements  of  an  inner  product  and  will  be  defined  as  the  inner  product  on  H: 


Let  n  =  {E  6  is  positive  semi-definite  diagonal}  and  let 

K  =  =  E^  where  E  €  Ji}.  Since  Q,  is  convex,  K  is  convex  also  (0  <  a  <  1): 

OfEi  "h  (1  —  0!)E2  G  fl 


therefore, 

[aEi  -I-  (1  -  a)E2]^  =  oEi^  +  (1  -  a)E2^  6  A' 

Since  K  is  equal  to  its  closure,  it  is  closed. 

The  inner  product  of  W  -  with  E#  —  C#  is 


(W-C$,E^-C#)  =  5^ 


«=1 


E  H  -  E  ^0^‘j  -  E  ^0  +  E  ^0 

i=l  j=l  ;=1  i=l 


For  each  i  such  that  52"=  i  >  0,  the  corresponding  term  in  the  summation  above  is 

zero; 


V  J  k  V  } 


i=i 


i=i 


-l-E.,E^0W;y  =0 

;=i 


For  each  i  such  that  5I"=i  <  0,  the  corresponding  term  in  the  summation  over  t  in 

the  inner  product  is  negative  for  all  E: 


n  n 

Eii  E  <0  for  E  <  0  and  E  e  if 

j=i  /=i 

Hence,  the  summation  over  t  in  the  inner  product  reduces  to  a  sum  of  zeros  amd  negative 
terms,  and  is  therefore  less  than  (or  equal  to)  zero  for  all  E  G  A”.  This  proves  that  C 
as  defined  in  Eq  (37)  minimizes  the  Frobenius  norm  ||C^  —  Wjl/-  over  the  set  of  positive 
semi-definite  diagonal  matrices. 

The  method  of  calculating  the  damping  coefficients  presented  here  is  useful  in  both 
Section  2.2.1  and  Section  3.2. 
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Appendix  B.  The  Value  of  the  Cost  Functional  in  Terms  of  the  Initial  State 

This  appendix  contains  the  derivation  of  a  form  of  the  cost  functional  in  terms  of  the 
initial  state,  which  is  the  form  given  in  Eq  (41).  Define  a  time  dependent  function  J{t)  as 


•^(0=  /  ^(y^Qy  +  u^R.u  + v^Sv)dt 
Jto  ^ 

Assume  J{t)  is  of  the  form 


(190) 


J(t)  =  c-y^Gy 

where  c  is  an  unknown  constant  and  G  is  an  unknown  matrix.  Then  the  derivative  of  J{t) 
with  respect  to  time  is 

j{t)  =  -y^Gy  -  y^Gy  -  y^Gy  (191) 

But  from  the  definition  of  J(t)  and  Eqs  (23)  and  (11)  in  Chapter  II, 

J(t)=  fy'^Qydr  (192) 

Jto 

where  Q  =  Q  -  PBR“‘B^P  +  Taking  the  time  derivative  of  Elq  (192)  gives 


j{t)  =  y'^Qy 


(193) 


Now  y  =  Ay  where  A  =  A  —  BR“*B’’P  +  BvC^.  The  objective  of  this  dissertation  is 
vibration  suppression,  so  P  and  C  have  been  chosen  such  that  A  is  stable.  Combining  Elqs 
(191)  and  (193)  gives 


y^Qy  =  -y’^AGy  -  y^Gy  -  y^GAy 


Now  G  =  0  since  Q  is  constant,  so  Eq  (194)  becomes  a  Lyapunov  equation  in  G: 


(194) 


A^G  +  GA  +  Q  =  0 


(195) 


In  order  to  express  J{t)  in  terms  of  the  initial  state,  yo,  the  state  y  is  first  expressed  in 
terms  of  the  initial  state: 

y  = 
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Then  the  cost  functional  can  be  written  as 


J{t)  =  c  - 

Since  the  initial  value  of  the  cost  functional  is  zero,  J{to)  =  0,  the  constant  c  is  c  =  yj Gyo- 

The  the  cost  functional  J  given  in  Eq  (15)  can  be  expressed  as  As  tj  approaches 

infinity,  the  cost  functional  can  be  expressed  as 

J  =  lim  J(t) 

t->oo  '  ' 

Since  A  is  stable,  the  limit  exists.  Hence,  the  cost  functional  is 

J=  lim  J(t)  =  yjGyo  (196) 

t—*0O 

where  G  is  the  solution  of  Eq  (195).  Comparing  Eqs  (40)  and  (195),  it  becomes  clear  that 
G  =  P.  Therefore  Eq  (196)  is  equivalent  to  Eq  (41). 

A  similar  development  can  also  be  found  in  Reference  [1]. 
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Appendix  C.  Solving  the  Non-Symmetric  Quadratic  Matrix  Equation 

In  this  appendix,  the  method  for  solving  Eqs  (108)  and  (113)  is  described.  The 
method  will  be  demonstrated  using  Eq  (108): 

Q  +  PA  +  A^P  -  PBiR-iBflxPf,I,2  -  PBvS"‘Bv^P  =  0  (108) 

This  equation  is  quadratic  in  P,  but  is  not  solvable  by  any  of  the  standard  MATLAB 
routines.  Newton-Raphson  is  a  standard  technique  and  it  is  employed  as  follows:  First  the 
left  hand  side  of  Eq  (108)  is  defined  as  a  function  of  P: 

Q(P)  =  Q  +  PA  +  A^P  -  PBiR-‘Bf  lAPIrIi2  -  PBvS-'Bv^P 

The  goal  is  to  find  P  such  that  Q(P)=0.  The  function  Q(P)  is  approximated  by 

Q(P)  =  Q(Po)+Z)QlPo](P  -  Po)+o(||P  -  Poll) 

where  the  operator  2?Q[Po]  represents  the  Frechet  derivative  of  Q  at  Po  and  o(P)  repre¬ 
sents  terms  that  are  at  least  quadratic  in  (||P  —  Poll)*  Let  H  =  P  —  Po* 

Since  Q  is  a  function  of  P, 

ZlQ[Po](H)  =  £q(Po  +  aH)|„=o  (197) 

Carrying  out  the  calculation, 

Z?Q[Po](H)  =  HA  -I-  A^H  -  HBiR-‘BriAPoIrIi2  -  PoBiR-^B^UKUiz 
-HBvS“‘b/Po  -  PoBvS"‘Bv^H 

This  may  not  look  much  better  than  Q(P),  but  at  least  it  is  linear  in  H.  The  equation  to 
solve  for  H  is 

Z?Q[Po](H)  +  Q(Po)  =  0  (198) 

Note  that  this  equation  is  similar  to  a  Lyapunov  equation,  although  some  of  the  terms  do 
not  have  H  on  the  left  or  right,  but  in  the  middle.  To  get  around  this  difficulty,  note  that 
I,  is  an  identity  matrix  with  the  bottom  third  diagonal  all  zeros.  Therefore  approximate 
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HI,  and  HIu  by  H  where  appropriate  in  I?Q[Po](H).  This  results  in  the  simplification, 


DQ(Po](H)  =  HA  +  A^H  -  HB,R-‘B[lAPoIrIi2  -  PoBrUH 
-HBvS'^Bv^’Po  -  PoBvS-‘B/H 

With  this  approximation,  Eq  (198)  becomes  a  Lyapunov  equation,  easily  solved  by  the 
simple  MATLAB  routine  LYAP.  The  initial  guess,  Po,  is  found  by  setting  I^,  1^.,  and  I12 
equal  to  identity  in  Eq  (108)  and  solving  the  resulting  algebraic  Riccati  equation.  The 
iteration  scheme  is 


Z)Q[P.1(H<)  +  Q(P),  =  0 
Pi+i  =  P.  +  H. 

Convergence  occurs  when  Q(P)  is  sufficiently  small.  The  sum  of  the  absolute  values  of  the 
elements  of  Q(P)  was  used  as  a  measure  of  “smallness”. 

The  method  presented  here  for  solving  the  non-symmetric  quadratic  matrix  equation 
can  be  used  to  solve  both  Eq  (108)  of  Section  3.2  and  Eq  (113)  of  Section  3.2.1. 


96 


Appendix  D.  Error  Analysis 

In  this  appendix,  the  error  between  the  approximation  Ay  and  the  operation  A(t,y) 
is  examined  in  detail.  The  operator  A(t,  y )  is  defined  by  Eqs  ( 163)  and  ( 164)  in  Chapter  IV. 
The  approximate  state  matrix  is  defined  by  Eq  (166): 


A  = 


(166) 


The  impulse  response  to  Eq  (160),  represented  by  x(<),  will  be  used  to  illustrate  the 
error  due  to  the  approximation  of  A(t,y).  Prom  Reference  [3],  the  impulse  response  can 
be  expressed  in  terms  of  an  integral  that  decays  algebraically  and  a  sum  of  decaying 
exponentials: 


X(()  =  i/m  [jT  X(re-‘')e-'‘rfr]  + 


^  {1.}(1  +  hA*"") 


mi 


X^t 


(199) 


where 


re 


^  mj{r^/”*e-‘’/”'  -  A,)  ^ 

a  =  order  of  fractional  derivative  (a  =  q/m) 

m  =  denominator  of  a 

Xj  =  eigenvalue  associated  with  the  expanded 
equations  of  motion  (see  Reference  [3]) 

=  eigenvector  associated  with  the  expanded 
equations  of  motion  (see  Reference  [3]) 
n  =  twice  the  length  of  the  vector  of  displacements 

N  =  total  number  of  eigenvalues  associated  with  the  expanded 
equations  of  motion  (see  Reference  [3]) 
b  =  viscoelastic  parameter  (see  Eq  (156)) 
mj  =  modal  constant 
{!.}  =  column  vector  of  ones 


The  approximate  state  matrix  is  calculated  using  the  eigenvalues  and  exponential 
modes  of  the  structure,  so  the  error  due  to  using  the  approximate  state  matrix  is  the 
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integral  term  of  the  impulse  response: 


Xb(<)  =  x(t)  -  x(t)  =  i/m  X(re-*')e-^'drj  (200) 

The  error  will  be  shown  to  be  real  and  continuous,  and  its  asymptotic  behavior  for 
small  and  large  time  examined. 

In  a  first  glance  at  the  behavior  of  the  error,  it  appears  to  be  singular  at  t  =  0. 
But  for  real  structures,  the  singular  components  cancel  each  other.  Churchill’s  theorem 
on  the  existence  of  the  inverse  transform  will  be  used  to  prove  that  the  response  is  real, 
continuous,  and  causal  [3:88].  Then  it  will  follow  that  the  error  is  real  and  continuous  since 
it  is  the  difference  between  two  continuous  functions. 

Before  applying  the  theorem,  it  will  be  shown  that  the  Laplace  transform  does, 
indeed,  satisfy  the  transformed  equations  as  they  appear  in  Eq  (160): 


s2Mx(s)  +  G(s)KvX(s)  +  Kx(5)  +  bu(s)  =  0  ( 160) 


The  Laplace  transform  of  the  impulse  response  x(t)  (Elq  (199))  for  which  -bu(<)  =  d(<)  is 


X(s)  [3:92]: 


N 


x(»)  =  2 


-  A,), 


(201) 


To  show  that  X(s)  does  indeed  satisfy  the  transformed  equations  of  motion,  the  equa¬ 
tions  of  motion  will  be  posed  in  expanded  form.  The  first  step  is  to  clear  the  denominator 
of  <j(s): 


G(s)  = 


Go  "I"  Gjs” 
1-1- 6s® 


Hence  Eq  (160)  becomes 


s*(l  -I-  6s“)Mx(s)  -f-  (Go  +  Gis“)Kvx(s)  -|-  (1  -I-  6s“)Kx(s)  =  -(1  -b  6s")bu(s)  (202) 

Rearranging  terms  yields 

[M(s2  +  6s^+“)  +  (GoKv  +  K)  +  (GiKv  +  6K)s“]  x(s)  =  -(1  -b  6s“)bu(s)  (203) 
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The  left  hand  side  of  this  equation  can  be  written  in  terms  of  summations; 


j  j 

+  x(s)  = -(1  +  6s“)bu(s)  (204) 

i=o  ,=o 

where  J  =  2q  +  m  and  Cj  and  K^-  are  zero  where  appropriate.  Letting  Aj  =  Mc^  -f- 
simplifies  Eq  (204)  even  further: 

J 

53  A,V/'"x(s)  =  -(1  +  bs“)huis)  (205) 

3=0 

Hence,  the  equations  of  motion  can  now  be  posed  in  terms  of  two  real,  square,  symmetric 
matrices: 


This  equation  can  be  expressed  in  more  compact  notation: 


s‘/”‘Mx(s)  +  Kx(s)  =  -u(s)  (206) 

The  orthogonal  transform  matrix  (or  modal  matrix)  #  which  simultaneously  diag¬ 
onalizes  the  expanded  mass  and  stiffness  matrices  (  M  and  K  )  is  constructed  from  the 
eigenvectors  associated  with  the  eigenvalue  problem  for  the  expanded  equations  of  motion: 

-  0  (207) 
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The  expanded  displacement  vector  can  be  expressed  in  terms  of  the  modal  matrix  and  a 
set  of  modal  coordinates  a(s): 

x(s)  =  ♦a(s)  (208) 

Premultiplying  Eq  (206)  by  and  substituting  in  Eq  (208)  produces  the  decoupled 
expanded  equations  of  motion  [24:159]  [3:68]: 

s‘/'"^^M*a(a)  +  4^Kia(s)  =  -4^u(5)  (209) 

To  show  the  decoupled  nature  explicitly,  the  equations  of  motion  are  written  in  terms 
of  the  modal  constants  m„  and  kn  where 


m„  =  (210) 

K  =  (211) 


Hence,  Eq  (209)  becomes 


gl/m 

nin 

»(«)  + 

kn 

a(«)  =  -♦’’u(s) 

(212) 

. 

Premultiplying  this  equation  by  the  matrix  di^l/m„)  yields 


^‘/"•^(s)  -1- 

k„lm„ 

a(s)  =  - 

l/m„ 

•  , 

^^u(s) 


(213) 


Premultiplying  Eq  (207)  by  and  recalling  Eqs  (210)  and  (211)  results  in  an  equa¬ 
tion  relating  the  eigenvalues  to  their  corresponding  modal  constants: 


“b  —  0 


(214) 
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Hence,  Ek]  (213)  can  be  rewritten  in  terms  of  the  eigenvalues: 


* , 

‘  , 

A„ 

a(s)  =  - 

l/m„ 

*  , 

4^u(s) 


(215) 


The  individual  elements  of  a(s)  look  like  modal  participation  factors  —  in  fact,  they 
axe.  The  individual  elements  of  a(s)  are  of  the  form: 

-  = 

m„(s»/"‘  -  A„) 

FVom  the  form  of  x(s)  and 


x(s) 

x(s) 

^  L 

a  « 

• 

II 

e 

; 

s~x(s) 

x(s) 

it  follows  that 

“(a)  =  (1  +  (216) 


The  Laplace  transform  of  the  structural  displacements  is  of  the  form 

*(^)  =  53«n(a)l^„ 

nssl 


or 


x(«)  = 


n=l 


*n„(s»/”*  -  A„) 


Therefore,  X(s)  is  the  transform  of  the  impulse  response  to  Elq  (160). 

Churchill’s  theorem  on  the  existence  of  the  inverse  transform  can  now  be  used  to 
prove  that  the  response  is  real,  continuous,  ard  causal. 
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s  that  is  analytic 
>  1;  also  let  f(x) 
7,  where  7  >  Xo, 

F{t)  =  !-'{/(«)}  (-00  <  <  <  00), 

whose  Laplace  transform  is  the  given  function  f{s): 

Hm)  =  /w  (*  >  *0). 

Furthermore  F{t)  is  0(e*®‘),  it  is  continuous  (—00  <  t  <  00),  and 

F(t)  =  0  when  <  <  0. 


Theorem  [9:178].  Let  f{3)  be  any  function  of  the  complex  variable 
and  of  order  0(«“‘)  for  all  s  {s  =  x-\-  iy)  in  a  half  plane  x  >  xq,  where  k 
be  real  (x  >  Xq).  Then  the  inversion  integral  of  f{s)  along  any  line  x  = 
converges  to  a  real-valued  function  F{t)  that  is  independent  of'), 


The  quantity  is  analytic  in  the  half  plane  x  >  Xo  when  Xq  is  positive  and  the 
branch  cut  of  s^!”*  is  chosen  to  lie  along  the  n^ative,  real  axis  in  the  s-plane  and  the  poles 
of  which  occur  at 

do  not  appear  in  the  right  half  s-plane.  (If  the  poles  appeared  in  the  right  half  s-plane, 
then  the  viscoelastic  material  would  be  generating  energy  instead  of  dissipating  energy.) 

The  quantity  X(x)  is  real  for  x  positive.  The  only  quantities  appearing  in  Eq  (201) 
that  can  be  complex  are  m^,  rjtj,  and  because  (l  +  fcx®)  and  x^!”*  are  real  for  x  positive. 
When  mj,  and  A^-  are  complex,  they  occur  in  conjugate  pairs  [3:89].  Hence,  it  follows 
directly  that  when  the  terms  in  Eq  (201)  sxe  complex,  they  occur  in  complex  conjugate 
pairs  and  X(x)  is  real  for  x  positive. 

X(s)  also  satisfies  the  condition  that  the  transform  be  of  order  0(s~*),  A;  >  1.  The 
direct  proof  that  X(s)  is  order  s"’  for  s  large  in  the  right  half  s-plane  (i.e.,  putting  all 
the  terms  over  a  common  denominator  and  adding  them  up)  is  long  and  involved,  so  an 
indirect  approach  will  be  pursued.  The  tramsformed  equations  of  motion  as  they  appear 
in  Ek}  (160)  will  be  used: 

s*Mx(s)  +  G{s)KyX{s)  +  Kx(s)  +  bu(s)  =  0  (160) 
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With  Kt(«)  =  G(s)Kv  +  K,  the  transformed  equations  of  motion  for  simultaneous,  unit 
impulsive  loading  (-bu(s)  =  {1.})  are 

[s’M  +  Kt(s)]X(s)  =  {1.}  (217) 

The  only  terms  in  Kt(s),  other  than  the  constant  terms,  are  those  terms  proportional  to 
G(s),  which  behaves  like  a  constant  as  s  -»  oo.  As  a  direct  result,  Eq  (217)  reduces  to 

3’MX(s)  ~  {!.}  as  s  oo  (218) 


Since  the  elements  in  the  mass  matrix  are  constant,  X(s)  is  order  s~^  for  s  large  in  the 
right  half  s-plane. 

Since  X(s)  satisfies  the  conditions  of  the  theorem,  the  response  (Ek^  (199))  is  real, 
continuous,  and  zero  at  initial  time  t  =  0.  It  follows  that  the  error  (Eq  (200))  is  also 
continuous,  since  it  equals  the  total  response  minus  a  sum  of  decaying  exponentials.  It 
also  follows  (from  Eq  (199))  that  the  initial  value  of  the  error  is  given  by 


xb(0)  =  -  ^ 


mA^-* 


V>r{i.}(i  +  6A"»°) 


TTii 


Before  examining  the  short  time  behavior  of  the  error,  the  short  time  behavior  of 
the  response  will  be  derived.  The  response  will  be  shown  to  monotone  for  short  time; 
a  Tauberian  theorem  will  then  be  appled  to  determine  the  short  time  behavior  of  the 
response. 

The  monotone  nature  of  the  response  for  short  time  can  be  derived  by  examining  the 
equations  of  motion  (Eq  (162)  in  the  time  domain  wth  the  forcing  term  — bu(t)  replaiced 
by  the  impulse  function  tf(t){l.}: 

Mx(f)  G„/(t  -  T)x(r)dT  +  Kx(t)  =  ^(0{1}  (219) 

Integrating  this  equation  from  0~  to  t  yields  an  expression  in  terms  of  the  velocity: 

Mx(t)  +  Kv  /  G„j(t  -  r)x(r)dr  +  K  /  x(f)  =  {!.}  (220) 

Jo  Jo- 
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Now  G,ei(0  is  a  continuous  function.  This  is  evident  when  the  Laplace  transform  of  Grei(t) 
(see  Hlq  (159)  and  (161))  is  expressed  as  a  sum  of  two  terms: 


I{G„,(0}  =  — + 

s 


Go  -  Gi/fe 

s(H-  6s“) 


The  first  term  is  clearly  the  Laplace  transform  of  a  step  function,  while  the  second  term 
satisfies  the  conditions  of  Chuchill’s  theorem  (where  the  Xo  of  the  theorem  is  equal  to  6*/®). 
Hence,  as  t  ].  0  the  integrals  in  Eq  (220)  disappear  since  the  integrands  are  continuous: 


limMx(t)  =  {!.} 


or 

limx(t)  =  M“‘{1.} 

Hence,  there  exists  a  ^  >  0  for  each  e  >  0  such  that  the  absolute  value  of  each  element  of 
x(t)  minus  each  element  of  M~^{1.}  is  less  than  c  whenever  0  <  t  <  ^.  Written  in  terms 
of  vectors, 

-  €{1.}  <  x(t)  -  <  €{1.}  (221) 

Solving  for  velocity, 

M-‘{1}  -  e{l.}  <  x(t)  <  +  e{l-}  (222) 


So  for  the  smallest  e  such  that  an  element  of  (M“^{1.}  — e{l.})  or  (M“*{l.}4-c{l.})  equals 
zero,  there  exists  a  j  >  0  such  that  the  elements  of  the  velocity  x(t)  do  not  change  sign 
on  0  <  t  <  ^.  Therefore,  each  element  of  the  response  x(t)  is  a  monotone  function  on  the 
interval  0  <t  <  6. 

To  examine  the  behavior  of  the  response  for  short  time  in  more  detail,  the  following 
Tauberian  theorem,  slightly  parapharsed  from  Feller,  will  be  used.  (Note  that  a;(s)  is  the 
Laplace  transform  of  u(t).) 

Theorem  [13:446].  IfU  has  an  ultimately  monotone  derivative  u  then  as  s  oo  and 
t  — ►  0,  respectively, 

w(a)~4  (223) 

s 
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FVom  Eq  (218),  the  Laplace  transform  of  the  response  for  large  s  is 


X(s)  -  {1.} 

It  follows  that  the  response  is  linear  for  short  time: 


x{t)  =  L  ‘{X(s)}  ~  M~*  {1.}  <  when  t  — >  0 

Since  the  error  is  the  response  minus  the  exponential  portion  of  the  response,  which  is 
known  exactly,  the  error  can  be  expressed  as 


Xsit)  ~  M  ‘  {1.}  <  -  ^ when  t  0 


It  is  appropriate  to  expand  the  exponential  terms  in  terms  of  t.  To  make  things  easier, 
the  vector  coefficient  of  will  be  represented  by  Cj.  The  vector  coefficients  appesj'  in 
complex  conjugate  pairs  and  the  complex  conjugate  of  Cj  will  be  assumed  to  be  c„/3^). 
The  coefficient  of  t  in  the  exponential  will  be  expressed  in  terms  of  its  real  and  ima^nary 
parts: 

K  =  -Pi  + 

Adding  the  complex  conjugate  pairs  and  carrying  out  the  expansions  yields  an  expression 
for  the  error  in  powers  of  t.  Retaining  only  the  first  power  of  t  yields 


n/2 


XB(t) - 2y^iZe(cj)  + 


J=1 


n/2 


M  ‘{l.}  +  2y3(p,i2e(c,)  +  7i/m(c,)) 


j=i 


for  small  time.  In  long  form,  the  constant  term  is 


-2^  J!e(c,)  =  -2^  m  (1  +  Me(A”“))  Re 

;=1  i=l  \  ^i  / 


which  is  equal  to  the  initial  value  of  the  error,  Xb(0). 

To  examine  the  behavior  of  the  error  for  large  time,  Watson’s  lemma  will  prove  useful. 

Watson’s  lemma  [10:446].  Let  F{r)  =  f{r“)r*,  a  >  0,  6  >  —  1.  Let  /(x)  have  a 
Maclaurin  expansion  for  jxj  <  6.  Let  |F(r)|  <  Me”’  for  some  constants  M  and  c  as 
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r  — ►  oo.  Let  f(x)  be  continuous  for  all  values  of  x.  Then 


F{r)e-'*dT  ~ 

t=o 


aiT{ka  +  6+1) 

^ia+i+t 


as  t 


<x>  where  = 


k< 


With  f{x)  defined  as 


^  -  Aj) 

the  quantity  X(re“**)  in  Eq  (200)  can  be  written  as 

X(re-")  =  /(r*/") 


The  Maclaurin  expansion  for  f(x)  is 


jti  mjXj 


1-^—, - ® 


+  •••+*«  ^6c-"“  + 


+ 


“T”/ 


-iw(i+3)/m 


¥^r-j 


for  |x|  <  |Ai|,  where  |Ai|  is  the  eigenvalue  with  the  smallest  magnitude.  Let  the  coefficient 
of  I*  be  represented  by  a*: 


hi  m,A,  V  A,  ) 

^  f  he-"*/"* 


for  0  <  <  q 

ioT  k>q 


Then  the  Maclaurin  expansion  can  be  written  more  compactly: 


/(*)  = 


k=0 
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Hence,  the  asymptotic  behavior  of  the  integral  in  the  error  is  given  by  a  summation  in 
terms  of  time  t: 


•'o  k=a  * 


a*r(fc/m  +  1) 


for  t  — >  oo 


(224) 


Expanding  the  first  few  terms  of  the  summation  in  Eq  (224)  yields 
r  ao  air(l  +  l/m)  a2r(l  +  2/m) 


for  t  — »  oo 


Now  Oo  can  be  expressed  in  terms  of  the  eigenvectors  and  modal  constants  m^  : 

<*0  =  1^ -  I - 

j=t  mjXj 

Since  the  eigenvectors,  eigenvalues,  and  modal  constants  occur  in  conjugate  pairs,  oo  is 
real.  Hence,  the  error  does  not  have  a  l/<  term: 

,,,  1  ,  /air(l  +  1/m)  ,  a^ril  +  2/m)  .  ^ 

To  show  that  the  imaginary  part  of  ci  is  non-zero,  oj  will  be  computed.  From  the  form  of 


m^A,  V  Xj  J 


Applying  Euler’s  identity  yields 

oi  =  >  - ^  (=08  ■Kim  -  t  sm  x/m) 

'njiXjY 

Hence,  by  the  same  rationale  that  qq  is  real,  Oi  has  an  imaginary  part: 

J!!, 

/m(ai)  =  X;  sin  x/m 

m,(A^)2 

Since  the  error  is  asymptotic  to  for  large  t,  it  does  not  decrease  exponentially. 

The  rationale  for  ignoring  the  error  due  to  the  approximation  is  given  at  the  end 
of  Section  4.1.  Although  the  solution  procedure  outlined  in  Section  4.2  does  not  take 
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into  account  the  relaxation  modes  (i.  e.,  the  error)  due  to  the  presence  of  viscoelastic 
dampers,  the  calculation  of  the  oscillation  modes  does  take  into  account  the  effect  of  the 
viscoelasticity  of  the  damping  material. 
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Appendix  E.  Alternate  Method  for  Calculating  the  State  Matrix 

This  appendix  contains  an  alternate  method  of  calculating  the  approximate  state 
matrix  to  the  one  given  in  Section  4.2.  The  derivation  of  this  method  makes  it  clear  that 
the  approximate  state  matrix  is  real-valued. 

In  Section  4.2  A  is  constructed  using  the  principal  value  eigenvectors  and  eigenvsdues 

^  =  _ 

where  A  is  a  diagonal  matrix  of  the  eigenvalues  and  'i'  is  the  matrix  of  corresponding 
eigenvectors.  Then  A'i'  =  'i'A.  Writing  A  in  terms  of  four  submatrices  A,-^,  this  becomes 

All  Ai2  ^  '*'1' 

Aji  A22  ’*'A  WK  9A  WA 

Carrying  out  the  matrix  multiplication 

Aii^' -f  Aia’i'A  Aii¥-1-Ai2^ 

Aji^'  -I- A22’*'A  A21^-1-A22^ 

results  in  four  equations.  The  top  two  will  be  considered  first: 


Aii^'  -I-  Aia^'A  =  9 A  (227) 

Aii^  -I-  Ai2^  =  ^  (228) 

Adding  and  subtracting  these  two  equations  gives 

Aii(^-|-¥)-1- Ai2(’*^A-1-^)  =  (^A-H^)  (229) 

Aii(«^  -  ^)  -I-  Ai2(^A  -  ■^)  =  (^A  -  ^)  (230) 


A  =  4'A4' 
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The  quantities  in  brackets  are  either  two  times  the  real  part  of  ^  ('i'A)  or  2i  times  the 
imaginary  part  of  'J'  Hence,  these  two  equations  can  be  written  as 


Aiii2e(^)  +  A,2/te('J'A)  =  Re(^A) 
Anlm(9)  +  Ai2/m('i'A)  =  Im(9A) 


Solving  for  An  in  each  of  the  equations  yields 


An  =  [Re(^)  -  Ai,Re(^A)]  [Re(^)]-^  =  [/m('J')  -  A,2/m(«'A)]  [/m('t')r‘  (231) 


Combining  the  terms  containing  An  gives 


Ai2  {iieC'J'A)  [Re{9)]~^  +  /mC^'A)  [/mC’J')]"*} 


(232) 


From  this  equation,  it  is  apparent  that 


A,2  =  I 

which,  from  Eq  (227),  implies 

All  =  0 


The  second  two  equations. 


A21’®'  +  A22«'A  =  ^A^ 
A21^  4-  A22^A  =  ^A^ 


(233) 


(234) 


(235) 

(236) 


result  in  the  two  equations. 


AiiiRc(^)  +  Ai2fle('®'A)  =  Re{9A^)  (237) 

Anlm{^)  +  Ai2/m('i'A)  =  /m('*'A^)  (238) 

Solving  for  A21  in  each  of  the  equations  yields 

A21  =  [Re{iif^)  -  A22M'®'A)]  [f2e('*')r‘  =  [/7n('®'")  -  A22/m(^A)]  [/m( «')]"'  (239) 
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Hence, 


A22{i?c(’*'A)[iZc(1')r'  +  /m(4'A)  [/m(4')r‘} 

(240) 


Finally,  solving  for  A22  yields  an  expression  which  is  real. 


A22  =  {f2e(«'A")(f?e('®')r‘  +/m('J'A")[/m('I')r'}  (241) 

X  {fle'J'A)  [fle('t)r‘  +  Im{9A)  [/m('i')r‘}~‘  (242) 


Substituting  this  expression  into  Eq  (239)  above  yields  A21.  Therefore,  A21  is  real  also. 
Since  An  and  An  are  real  also,  the  state  matrix  A  is  also  real.  Since  this  matrix  is 
equivalent  to  the  approximate  state  matrix  derived  in  Section  4.2,  the  approximate  state 
matrix  is  also  real. 
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Appendix  F.  Feedback  Matrices  for  the  Truss  Examples 

This  appendix  contains  the  active  gain  feedback  matrices  for  all  the  fourteen  degree 
of  freedom  examples  in  this  dissertation.  E^h  feedback  matrix  has  dimension  4  x  28.  The 
gains  for  each  chapter’s  examples  are  presented  in  separate  sections. 

F.l  Chapter  II  Gains 

This  section  contains  the  feedback  matrices  as  computed  by  £q  (11)  for  the  four 
different  solution  techniques  developed  in  Sections  2.2.1  and  2.2.2  for  Example  #3.  It  also 
includes  the  feedback  matrices  for  the  two  solution  techniques  presented  in  Section  2.3  for 
Example  #4.  The  matrices  aure  on  the  foUowing  pages.  Table  19  is  the  key  to  the  solutions. 


Table  19.  Key  to  Matrix  Names  for  Section  F.l 


Example 

Matrix  Name 

Solution  Technique 

Example  #3 

Gf 

min  ||H||r 

Gih 

min  lIHlla 

Gtr 

min  trace  P 

Ga 

min  yPlla 

Example  #4 

G,.. 

min  ^  trace  P  +  Jc^Sc 

Gnio* 

niin  JdlPlIa  +  c^Sc) 
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G/  = 


Colunns  1  through  7 


60.8854 

137.2134 

-75.7822 

-173.6511 

26.2906 

-34.7551 

-8.8621 

69.6113 

135.6511 

-52.4885 

14.8513 

-147.8120 

-106.1664 

168.0595 

78.5921 

104.6444 

18.1517 

-62.9460 

-22.1549 

-298.9220 

103.0132 

235.7793 

136.3770 

112.7833 

50.3280 

72.5824 

-59.3887 

-143.9000 

Colunns  8 

through  14 

48.0909 

-8.4912 

33.3340 

-2.3179 

-32.7392 

5.9136 

7.3087 

19.5168 

132.7994 

30.9326 

-134.7725 

-65.6421 

-21.4194 

41.8575 

312.1403 

-23.4110 

130.7269 

-71.1841 

-175.1829 

42.5333 

20.4618 

-4.5445  - 

-170.0058 

-15.7503 

108.4971 

24.7631 

24.4720 

-10.5425 

Colunns  15  through  21 

6.1961 

-1.4179 

0.1507 

-0.1824 

1.0709 

-0.2874 

1.4822 

-1.4142 

11.4926 

0.7460 

3.4615 

0.3602 

2.7212 

-0.7177 

0.1534 

0.7517 

7.6018 

-0.8572 

2.7187 

1.2781 

1.9533 

-0.1758 

3.4552 

-0.8656 

11.9417 

-0.0934 

6.7203 

0.9148 

Colunns  22  through 

28 

-1.6310 

1.5495 

-1.7135 

1.0639 

-1.3217 

2.6922 

-3.8026 

5.8867 

0.1494 

3.8181 

0.2867 

3.0427 

0.8136 

6.3550 

1.4761 

1.7575 

0.6129 

2.9212 

1.4884 

6.1589 

2.5689 

3.1828 

0.7333 

3.3875 

0.3460 

4.1678 

1.4643 

6.3543 
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Colims  1  through  7 


62.9054 

134.8581 

-66.3055 

-176.0679 

18.3592 

-26.0768 

-16.5618 

79.3484 

144.3491 

-37.6207 

21.9677 

-166.3707 

-113.2804 

163.3008 

68.4678 

95.8321 

16.8779 

-68.9060 

-18.3349 

-266.4717 

67.5165 

243.8168 

145.9967 

122.4900 

58.4604 

71.4145 

-89.2929 

-136.6324 

Colunns  8  through  14 

38.2794 

-4.4152 

42.6332 

-1.6232 

-39.9349 

3.4473 

6.7362 

-12.4108 

144.4190 

82.1159 

-131.7938 

-70.3460 

-12.1800 

23.5068 

271.8900 

0.4167 

180.2037 

-69.4202 

-221.3717 

38.4422 

23.2290 

-10.0895 

-167.2828 

-17.2681 

118.7133 

69.9598 

23.8387 

-30.7180 

Coluams  16  through 

21 

6.2410 

-1.4277 

0.1415 

-0.1732 

1.0468 

-0.3986 

1.6337 

-1.4234 

11.7814 

0.7319 

3.7110 

0.4766 

2.5659 

-0.6713 

0.1439 

0.7369 

7.6211 

-0.9163 

2.6754 

1.2443 

1.9206 

-0.1663 

3.7063 

-0.9241 

12.2309 

-0.0261 

6.5101 

1.0204 

ColuBus  22  through 

28 

-1.6967 

1.5318 

-1.7559 

1.0256 

-1.2416 

2.7222 

-3.8331 

5.6853 

0.1806 

3.9373 

0.3396 

3.1580 

0.9155 

6.2907 

1.5366 

1.7642 

0.6783 

3.0228 

1.3866 

6.0855 

2.6713 

3.0823 

0.7998 

3.4006 

0.3364 

4.3195 

1.5819 

6.3279 

G.,= 


Coluans  1  through  7 


62.5198 

136.9156 

-84.9311 

-173.4920 

19.8862 

-36.8601 

5.6607 

70.7117 

138.3706 

-63.3180 

17.3930 

-158.8348 

-102.9273 

174.5300 

88.3765 

113.3171 

19.2687 

-63.9681 

-29.4842 

-302.0402 

109.8547 

230.0153 

139.5345 

114.5842 

54.4877 

79.6198 

-75.2200 

-152.8309 

Columns  8 

through  14 

53.4789 

5.0294 

35.8599 

-9.0953 

-40.5677 

-0.1541 

10.6292 

5.0429 

144.8681 

20.4038 

-144.3306 

-55.2011 

-34.8869 

46.1979 

318.4706 

-25.8206 

128.6146 

-63.7649 

-179.7290 

34.9756 

21.2320 

-0.1763  - 

■180.1297 

-4.7210 

111.2053 

15.1257 

14.5805 

-9.0842 

Coluans  15  through 

21 

6.1509 

-1.3413 

0.1804 

-0.3025 

1.1789 

-0.4110 

1.5747 

-1.3369 

11.7083 

0.7295 

3.7868 

0.3732 

2.7948 

-0.7119 

0.1831 

0.7350 

7.7083 

-0.8737 

2.8665 

1.2882 

2.0999 

-0.2958 

3.7803 

-0.8823 

12.1434 

0.0067 

6.6483 

0.9537 

Coluans  22  through 

28 

-1.4605 

1.6688 

-1.6170 

1 . 1401 

-1.3985 

2.8010 

-3.8226 

5.8102 

0.1104 

3.7480 

0.3557 

3.0772 

0.9515 

6.3343 

1.6487 

1.8831 

0.7459 

3.1165 

1.4338 

6.5630 

2.4971 

3.2901 

0.8715 

3.4161 

0.3831 

4.0682 

1.5877 

6.3704 
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Gj  = 


Coluans  1  through  7 


59.4241 

111.0129 

-66.6250 

-148.7496 

11.8285 

-44.6229 

10.5544 

129.2189 

166.2954 

-86.3966 

9.1834 

-162.3672 

-72.3976 

144.7425 

85.9633 

125.1300 

20.7201 

-50.8535 

-34.4722 

-309.3312 

142.5468 

244.6480 

165.7833 

95.4493 

54.5290 

101 . 1949 

-136.9265 

-144.0311 

Coluans  8 

through  14 

65.7956 

0.8674 

22.3652 

-15.4089 

-36.6680 

2.0923 

13.4014 

-59.8052 

148.0188 

9.7843 

-137.8634 

-30.7686 

-80.7885 

61.7748 

323.0306 

-52.1648 

133.3272 

-62.2221 

-194.8828 

40.4834 

20.4625 

26.2492  - 

■201.3158 

13.9502 

82.8149 

-1.6140 

-13.0974 

1.0109 

Columns  15  throiigh 

21 

6.4417 

-1.3761 

0.1911 

-0.3406 

1.1791 

-0.4620 

1.5320 

-1.3712 

11.9444 

0.6470 

3.8516 

0.2777 

2.5384 

-0.7196 

0.1939 

0.5529 

8.1370 

-0.8070 

2.7663 

1.2321 

2.0935 

-0.3339 

3.8460 

-0.8164 

12.4176 

-0.0346 

6.1569 

0.8482 

Columns  22  through 

28 

-1.3896 

1.6949 

-1.5766 

1.0828 

-1.3600 

2.7290 

-3.8698 

5.3474 

-0.0081 

3.8167 

0.2607 

3.1503 

0.6429 

6.4370 

1.9055 

1.7932 

0.8369 

3.0411 

1.4902 

6.6023 

2.3926 

3.0848 

0.6363 

3.4532 

0.2701 

4.1371 

1.3288 

6.5021 
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Gawe  = 


Columns  1  through  7 


64.1695 

144.5988 

-50.0861 

-183.6278 

17.4595 

-16.5797 

-57.6956 

177.3899 

163.9929 

-69.1734 

-14.5471 

-200.8170 

-104.2373 

213.3795 

74.6711 

131.9043 

-7.9591 

-94.1001 

-13.6606 

-307.3552 

127.4211 

328.3747 

171.8565 

134.3263 

6.8416 

145.0328 

-40.5235 

-181.5181 

Columns  8 

through  14 

23.6747 

-34.7961 

22.5202 

29.2920 

-17.9947 

10.8837 

3.2705 

41.7881 

168.7589 

0.6942 

-159.4808 

-55.8531 

-57.0812 

51.9287 

323.8085 

-3.0948 

126.3508 

-92.6014 

-179.7254 

45.7161 

30.1199 

13.7485  • 

-230.9777 

-6.5328 

124.6472 

6.2781 

18.9491 

-10.7131 

Columns  15  through 

21 

2.6100 

-2.0565 

0.0589 

0.3003 

0.9235 

-0.2244 

1.2754 

>2.0540 

7.7074 

0.8255 

2.3542 

0.2910 

2.5540 

-1.0319 

0.0631 

0.8334 

4.5988 

-0.9877 

2.4842 

1.3050 

1.4064 

0.3068 

2.3489 

-0.9968 

8.3288 

-0.6151 

6.8957 

0.6598 

Columns  22  through 

28 

-1.8606 

1.3283 

-1.7513 

0.9407 

-1.3099 

2.3204 

-3.8801 

6.0478 

-0.0930 

3.9044 

-0.1169 

3.0212 

-0.0418 

6.6334 

1.0905 

1.2959 

0.5579 

2.4848 

1.6559 

5.2348 

2.6133 

2.9706 

0.1518 

3.4567 

0.0377 

4.3167 

0.5836 

6.5274 
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! 


G 


mac  — 


Coliums  1  through  7 


92.3936 

154.8387 

-58.1994 

-199.3819 

-67.3366 

20.0614 

195.5806 

171.6741 

-94.6911 

-43.1889 

-181.9405 

-79.6736 

105.0575 

141.6966 

-45.9438 

-166.4785 

8.7368 

-164.3084 

361.4868 

183.6536 

127.2080 

-27.1159 

225.2402 

-36.0190 

ColuianB  8 

through  14 

-25.4892 

17.9443 

5.2245 

33.1036 

-0.0822 

-27.6276 

0.4522 

272.0096 

100.9069 

-140.9888 

-58.9851 

-35.6951 

238.3149 

25.3230 

188.2920 

-82.7019 

-268.2625 

41.4951 

22.7196  ' 

-239.5681 

-9.4620 

126.9080 

88.0696 

73.3432 

Columns  15  through 

21 

2.8996 

-2.0661 

0.0321 

0.4696 

0.9380 

-0.7617 

-2.0612 

7.8104 

1.0140 

2.1628 

0.5056 

3.0963 

0.0378 

1.0808 

4.8366 

-1.3009 

2.9930 

1.5393 

0.4781 

2.1757 

-1.3070 

8.2418 

-0.0999 

6.5909 

Columns  22  through 

28 

-1.5725 

1.2578 

-1.5825 

0.8400 

-1.3211 

1.6844 

5.5949 

-0.0739 

3.4076 

0.1453 

2.8400 

0.8144 

0.5598 

0.9202 

0.9424 

2.3418 

1.3532 

4.3504 

3.5836 

0.4557 

3.0788 

0.2729 

3.7776 

1 . 1063 

-58.6918 

212.1390 

74.2271 

-144.7645 


8.3576 

10.2324 

47.1680 

-59.2934 


1.2350 

-0.8654 

0.8631 

0.8618 


-3.8714 

6.8125 

2.5067 

6.8835 
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F.2  Chapter  III  Gains 


This  section  contains  the  feedback  matrices  as  computed  by  Eq  (99)  for  the  three 
different  solution  techniques  developed  in  Sections  3.2  and  3.2.1  for  Example  #3.  The 
matrices  are  on  the  following  pages.  Table  20  is  the  key  to  the  solutions. 

Table  20.  Key  to  Matrix  Names  for  Section  F.2 


Matrix  Name 

Solution  Technique 

G/ 

min  IjHllf 

min  trace  P 

Gj 

min  IIPII2 
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G/  = 


Columns  1  through  7 


489.0862 

28.6248 

-71.3414 

-1.5727 

-72.1660 

37.4579 

37.6215 

17.7314 

515.4182 

109.0362 

-297.4294 

31.6795 

-146.2687 

-20.8456 

-44.9506 

87.4180 

604.7917 

-50.9229 

52.4029 

-85.1045 

-90.0218 

105.2035  • 

■271.3641 

-18.1287 

468.7974 

-19.0622 

97.2782 

64.8550 

Columns  8 

through  14 

-26.5444 

7.8820 

0.3209 

-3.7119 

0.8848 

-2.5209 

9.5871 

130.7301 

25.1807 

21.3275 

8.7155 

-38.4993 

33.8261 

1.2374 

59.3063 

-56.4086 

-14.8269 

-10.3330 

-0.5436 

24.4887 

1.9621 

-120.3438 

1.8507 

-26.0600 

24.9368 

24.5943 

32.8938 

-12.1145 

Columns  15  through  21 

5.4058 

-0.1816 

-0.4012 

0.2324 

0.0592 

0.2144 

0.8973 

-0.2802 

8.3526 

0.9596 

-0.1919 

0.1924 

0.7148 

-0.9664 

-0.4366 

0.9121 

6.7926 

-0.7942 

1.2161 

-0.3428 

-0.0866 

0.4123 

-0.1276 

-0.7620 

7.8868 

-0.7648 

3.7372 

0.5457 

Columns  22  through 

28 

-0.7086 

0.7830 

-0.6640 

0.3901 

-0.6744 

1.1016 

-1.6770 

3.7738 

-0.3749 

2.2375 

-0.1275 

1.6387 

-0.3797 

3.8566 

0.4802 

-0.1450 

0.0443 

0.8866 

0.5386 

1.9161 

1.4046 

0.8781 

-0.1590 

1.8312 

-0.3545 

2.3890 

-0.4807 

3.6853 
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G.,  =  lO^x 


Columns  1  through  7 


5.6539 

0.6351 

0.7863 

-0.0965 

-0.9973 

-0.2201 

0.4552 

4.2154 

-0.0477 

-2.5036 

0.0183 

-0.7224 

0.3678 

0.5871 

6.7266 

0.3822 

-0.6905 

-0.5118 

1.4910 

-1.8828 

0.8736 

3.7913 

-0.0626 

-0.2264 

Columns  8  through  14 

0.1717 

-0.2739 

-0.3260 

-0.7621 

0.2627 

0.2695 

0.2298 

-0.0481 

0.3316 

0.1400 

-0.3537 

-0.1322 

0.0057 

-0.5344 

0.2580 

-0.3523 

-0.4444 

0.0257 

-0.5208 

-0.5672 

-0.4576 

-0.1986 

0.3467 

0.3019 

Columns  15 

through 

21 

0.0861 

0.0089 

0.0119 

-0.0025 

0.0011 

-0.0074 

-0.0040 

0.0758 

0.0062 

-0.0166 

-0.0043 

0.0052 

0.0031 

0.0065 

0.1020 

0.0069 

-0.0021 

0.0018 

0.0135 

-0.0135 

0.0097 

0.0766 

-0.0042 

0.0198 

Columns  22 

through 

28 

-0.0025 

0.0018 

-0.0101 

-0.0018 

-0.0040 

0.0144 

0.0223 

0.0021 

0.0156 

0.0000 

0.0081 

-0.0037 

0.0004 

0.0049 

0.0027 

0.0078 

0.0037 

0.0215 

0.0100 

-0.0030 

0.0093 

-0.0014 

0.0164 

0.0056 

-0.3016 

0.5SS3 

-1.1454 

0.1935 


0.2142 

-0.0269 

-0.0366 

0.0697 


0.0018 

0.0017 

0.0045 

0.0003 


-0.0127 

0.0208 

0.0113 

0.0218 
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Gz  =  lO^x 


Columns  1  through  7 


6.0104 

0.2557 

0.9050 

0.3489 

-1,5184 

-0.7742 

-0.0029 

0.6289 

4.1351 

0.2488 

-2.2366 

-0.2825 

-1.0163 

0.4448 

0.3753 

0.6546 

6.8706 

0.5112 

-0.9782 

-0.7147 

-0.9912 

1.5668 

-1.8007 

0.8274 

3.8946 

-0.2072 

-0,6900 

0.2627 

Columns  8  through  14 

0.7349 

-0.1738 

-0.5530 

-0.9896 

0.3881 

0.3245 

0.2369 

0.0403 

0.1154 

0.5935 

-0.0895 

-0.3598 

-0.1564 

-0,0388 

0.1271 

-0.5036 

0.3764 

-0.3193 

-0.6004 

-0.1108 

-0.0595 

-0.4776 

-0,4433 

-0.5820 

-0.2490 

0.7063 

0.0476 

0.0356 

Columns  15 

through 

21 

0.0910 

0.0048 

0.0129 

0.0016 

-0.0032 

-0.0128 

0.0009 

-0.0041 

0.0765 

0.0063 

-0.0137 

-0.0048 

0.0037 

0.0004 

0.0053 

0.0065 

0.1017 

0.0063 

-0.0029 

0.0024 

0.0079 

0.0146 

-0.0123 

0.0054 

0.0773 

-0.0032 

0.0175 

0.0017 

Coliunns  22 

through 

28 

0.0024 

0.0037 

-0.0121 

-0.0049 

-0.0021 

0.0174 

-0.0127 

0.0202 

0.0018 

0.0175 

-0.0014 

0.0069 

-0,0051 

0.0203 

0.0046 

0.0065 

0.0036 

0.0093 

0.0039 

0.0217 

0.0108 

0.0117 

-0.0036 

0.0059 

-0.0025 

0.0195 

0.0024 

0.0213 
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F.3  Chapter  IV  Gains 


This  section  contains  the  feedback  matrices  as  computed  by  Eq  (11)  for  the  two 
different  solution  techniques  developed  in  Section  4.2  for  both  cases  of  the  fourteen  degree  of 
freedom  example  given  in  that  chapter.  The  matrices  are  on  the  following  pages.  Table  21 
is  the  key  to  the  solutions. 


Table  21.  Key  to  Matrix  Names  for  Section  F.3 


Example 

Matrix  Name 

Solution  Technique 

Gi  =  7.3 

G,, 

min  trace  P 

G2 

min  IIPII2 

Gi  =  730 

Gave 

min  trace  P 

Gma* 

min  IIPII2 
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Columns  1  through  7 


-11.7488 

3.2884 

13.4986 

-0.5558 

17.9921 

-13.8203 

-16.7188 

-9.9056 

-16.9476 

29.3157 

0.7307 

9.9056 

-19.5694 

-14.2028 

-16.6610 

-34.5369 

-12.3350 

26.5726 

7.9564 

26.8844 

16.1763 

-3.4407 

-13.1607 

-32.9885 

-16.3686 

0.2130 

47.2643 

12.5992 

Columns  8 

through  14 

14.0544 

-15.3505 

-16.4734 

15.3605 

18.9762 

-2.2705 

-3.6532 

39.0706 

-24.1762 

-4.5332 

8.6402 

-4.2683 

10.3667 

-2.2846 

-29.0592 

-3.8376 

15.8942 

-18.6049 

-11.6161 

10.3083 

-0.9486 

-23.9520 

23.5508 

-1.9967 

-7.1669 

-4.2658 

-2.3863 

3.4267 

Coliunns  15 

>  through 

21 

-5.1343 

0.1716 

-0.3069 

-0.0404 

-0.1499 

0.3782 

-0.4648 

0.1692 

-6.0402 

-0.2692 

-0.7329 

-0.2864 

-0.8601 

0.0000 

-0.3071 

-0.2697 

-5.4743 

-0.0852 

-0.6116 

-0.3406 

-0.5974 

-0.0422 

-0.7302 

-0.0836 

-6.2146 

0.0723 

-1.2644 

-0.6249 

Columns  22 

!  through 

28 

0.0626 

-0.3449 

0.4453 

-0.2715 

0.5341 

-0.9002 

1.1157 

-1.2409 

-0.1138 

-1.0375 

-0.4907 

-1.4666 

-0.6478 

-1.8196 

-0.5316 

-0.7306 

-0.2372 

-0.9703 

-0.1914 

-1.9094 

-0.7640 

-1.0689 

-0.4573 

-1.5007 

-0.0909 

-1.0637 

-0.9986 

-1.9607 
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Gj  = 


Columns  1  through  7 


-11.7936 

3.6822 

13.4141 

-1.1972 

19.0378 

-14.1206 

-15.4661 

-9.7794 

-16.1182 

29.9986 

0.9010 

10.6674 

-20.1413 

-15.9227 

-16.9676 

-36.4136 

-12.7730 

27.7664 

8.3364 

26.3197 

18.2846 

-2.4169 

-12.0716 

-34.9031 

-16.2340 

1.2867 

46.2696 

15.8351 

Columns  8 

through  14 

16  2768 

-18.1244 

-16.9983 

14.6906 

18.8738 

-0.9257 

-3.6643 

36  J842 

-23.8281 

-3.3703 

9.6093 

-2.9364 

8.9960 

-3.0381 

-26.6124 

-6.6469 

14.0736 

-20.2280 

-11.1046 

11.0074 

-0.6140 

-23.8286 

21.6736 

-1.6667 

-9.0406 

-3.0817 

-2.3342 

2.8039 

Columns  16  through 

21 

-6.1660 

0.1681 

-0.3420 

-0.0606 

-0.0630 

0.3829 

-0.4401 

0.1667 

-6.9862 

-0.2674 

-0.6662 

-0.2846 

-0.9341 

0.0109 

-0.3426 

-0.2679 

-6.6387 

-0.1240 

-0.4237 

-0.3109 

-0.4992 

-0.0624 

-0.6623 

-0.1222 

-6.1774 

0.1294 

-1.3122 

-0.6896 

Columns  22  through 

28 

0.0732 

-0.3376 

0.4428 

-0.3043 

0.6322 

-0.9076 

1.1218 

-1.3018 

-0.1148 

-1.0139 

-0.4733 

-1.3978 

-0.6668 

-1.8391 

-0.6020 

-0.7731 

-0.2466 

-0.9706 

-0.2047 

-1.9324 

-0.7673 

-1.1380 

-0.4672 

-1.4613 

-0.1094 

-1.0390 

-0.9881 

-1.9717 
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Gatie  — 


Columns  1  through  7 


-22.8709 

22.8106 

-6.5041 

-28.3029 

54.2380 

0.8188 

-20.0364 

-16.5469 

-20.1663 

65.4339 

6.3215 

4.6520 

-44.3262 

-21.4800 

-6.0726 

-69.4766 

-23.9072 

54.6779 

-5.3017 

1.4459 

69.3602 

16.8109 

-11.9736 

-61.8080 

-22.9949 

4.4964 

75.6195 

28.6704 

Columns  8 

through  14 

6.5656 

-20.1967 

-24.4718 

21.3485 

28.5355 

-8.7677 

-3,4032 

62.5311 

-31.8067 

21.0302 

3.1264 

-31.6115 

23.7250 

-2.5897 

-0.6412 

-24.9763 

48.2131 

-66.3328 

-52.2274 

32.1660 

0.1134 

-41.2860 

28.8172 

-23.4676 

-7.8490 

14.2654 

-8.1036 

4.3680 

Columns  15  through  21 

-7.4396 

0.5249 

-0.8452 

-0.7396 

0.1190 

0.9985 

-0.4188 

0.5212 

-7.0872 

-0.2929 

-0,4262 

-0.6486 

-0.9977 

-0.0066 

-0.8442 

-0,2941 

-7.6642 

-0.4416 

0.2043 

-0.4186 

-0,4776 

-0.7415 

-0.4217 

-0.4407 

-7.5236 

0.6443 

-0.5406 

-0.9764 

Columns  22  through 

28 

-0.6363 

-0.6846 

0.7457 

0.2731 

0.3761 

-0.9087 

0.9068 

-0.8074 

-0.1145 

-0.6607 

-0.9885 

-2.0818 

-0.2452 

-1.6534 

-0.3066 

-1.1602 

-0.2901 

-1.3372 

-0.0700 

-1.2076 

-0.7120 

-1.3012 

-0.5369 

-2.0333 

0.3319 

-0.6358 

-1.3158 

-1.9474 
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Gfna*  — 


Columns  1  through  7 


-22.8709 

22.8106 

-6.5041 

-28.3029 

54.2380 

0.8188 

-20.0364 

-16.5459 

-20.1653 

55.4339 

6.3215 

4.6520 

-44.3262 

-21.4800 

-6.0725 

-59.4756 

-23.9072 

54.6779 

-5.3017 

1.4459 

69.3602 

16.8109 

-11.9735 

-61.8080 

-22.9949 

4.4964 

75.5195 

28.6704 

Columns  8 

through  14 

6.5656 

-20.1957 

-24.4718 

21.3485 

28.5355 

-8.7677 

-3.4032 

62.5311 

-31.8067 

21.0302 

3.1264 

-31.6115 

23.7250 

-2.5897 

-0.5412 

-24.9753 

48.2131 

-56.3328 

-52.2274 

32.1660 

0.1134 

-41.2860 

28.8172 

-23.4676 

-7.8490 

14.2654 

-8.1035 

4.3680 

Columns  15  through 

21 

-7.4396 

0.5249 

-0.8452 

-0.7395 

0.1190 

0.9985 

-0.4188 

0.5212 

-7.0872 

-0.2929 

-0.4252 

-0.5485 

-0.9977 

-0.0066 

-0.8442 

-0.2941 

-7.6642 

-0.4416 

0.2043 

-0.4186 

-0.4775 

-0.7415 

-0.4217 

-0.4407 

-7.5235 

0.6443 

-0.5406 

-0.9764 

Columns  22  through 

28 

-0.5363 

-0.5846 

0.7457 

0.2731 

0.3761 

-0.9087 

0.9068 

-0.8074 

-0.1145 

-0.6507 

-0.9885 

-2.0818 

-0.2452 

-1.6534 

-0.3056 

-1.1602 

-0.2901 

-1.3372 

-0.0700 

-1.2076 

-0.7120 

-1.3012 

-0.5369 

-2.0333 

0.3319 

-0.6358 

-1.3158 

-1.9474 
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